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(54) Title: SIGNAL PROCESSING APPARATUS 
(57) Abstract 

The present invention involves 
method and apparatus for analyzing 
two measured signals that are modeled 
as containing primary and secondary 
portions. Coefficients relate the two 
signals according to a model defined in 
accordance with the present invention. 
In one embodiment, the present invention 
involves utilizing a transformation which 
evaluates a plurality of possible signal 
coefficients in order to find appropriate 
coefficients. Alternatively, the present 

invention involves using statistical , . . 

functions or Fourier transform and windowing techniques to determine the coefficients relating to two measured signals, 
invention is described in particular detail with respect to blood oximetry measurements. 
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SIGNAL PROCESSING APPARATUS 

PurtniTTH "f tiff '"«"tion 

p«M nf thP Invention 

The present invention relates t. the field of signal processing. More speafically. the present 
invention relates to the processing of measured signals, containing . primary signal portion and a secondary 
signal portion, for the removal or derivation of either the primary or secondary signal portion when little 
is known about either of these components. More particularly, the present invention relates to modebng 
the measured signals in a novel way which MM* minimrnng the correlation between the primary s,gnal 
portion and the secondary signal portion in order to produce a primary and/or secondary signal. The present 
invention is especially useful for physiological monitoring systems including blood oxygen saturation systems. 

np«rri ptinn of t *ff Belated Art / 

Signal processors are typically employed to remove or derive either the primary or secondary s.gnal 
portion from a composite measured signal including a primary signal portion and a secondary signal porfon. 
For example, a composite signal may contain noise and desirable portions. If the secondary signal portion 
occupies . different frequency spectrum than the primary signal portion, than conventional ffltenng 
techniques such as .ow pass, band pass, and high pass filtering are available to remove or derive «ther the 
primary or the secondary signal portion from the total signal Fixed single or multiple notch filters could 
also be employed if the primary and/or secondary signal portion® exist at a fixed frequency®. 

„ is often the case that an overlap in frequency spectrum between the primary and secondary 
signal portions exists. Complicating matters further, the statistical properties of one or both of the prunary 
and secondary signal portions change with *» In such case, conventional filter*, techmques are 
ineffectWe in extracting either the primary or secondary signal. If. however. . description of ether the 
primary or secondary signal portion can be derived, correlation canceling, such as adaptive noise cancelm* 
can be employed to remove either the prunary or secondary sign* portion of the signal .solafng the other 
portion, m other word, given sufficient information about one of the signal portion,that signal porfon 
can be extracted. . ' ... 

^ «^« 1. - ~ »^ -«"«•"'■ Howmw. correlation cancatas 

rao*. orU»r . —ft* <*— •< • I— » «*■—— '° * """^ 

Fw instinct, lor a Matured signal contammg noise and 

•» .««, «ha e«* AiU-v* th. «H of «- »<«™- «P* « "« -—-»» — •• 
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amplitude of the corresponding primary or secondary signal portions, they have a frequency spectrum which 

is similar to that of the primary or secondary signal portions. 

In many cases, nothing or very little is known about the secondary and/or primary signal portions. 

One area where measured signals comprising a primary signal portion and a secondary signal portion about 
5 which no information can easily be determined is physiological monitoring. Physiological monitoring generally 

involves measured signals derived from a physiological system, such as the human body. Measurements 

which are typically taken with physiological monitoring systems include electrocardiographs. Wood pressure. 

blood gas saturation (such as oxygen saturation), capnographs. other blood constituent monitoring, heart 

rate, respiration rate, electro^mcephalograph (EEG) and depth of anesthesia, for example. Other types of 
10 measurements include those which measure the pressure and quantity of a substance within the body such 

as cardiac output, venous oxygen saturation, arterial oxygen saturation, bffiruhin. total hemoglobin. 

breathalyzer testing, drug testing, cholesterol testing, glucose testing, extra vasation, and carbon dioxide 

testing, protein testing, carbon monoxide testing, and other in-vivo measurements, fx example. 

Complications arising in these measurements ere often due to motion of the patient, both external and 
15 internal (muscle movement, vessel movement end probe movement for example), during the measurement 

process. 

Many types of physiological measurements can be made by using the known properties of energy 
attenuation as a selected form of energy passes through a medium. 

A blood gas monitor is one example of a physiological monitoring system which is based upon the 

20 measurement of energy attenuated by biological tissues or substances. Blood gas monitors transmit light 
into the test medium and measure the attenuation of the Dght as a function of time. The output signal 
of a blood gas monitor which b sensitive to the arterial blood flow contains a component which is e 
waveform representative of the patient's erterial pulse. This type of signal, which contains a component 
related to the patient's pulse, is called • plethysmography wave, and is shown in Figure 1 as curve s. 

25 Plethysmography waveforms are used in blood gas saturation measurements. As the heart beats, the 
amount of blood in the arteries increases and decreases, causing increases and decreases in energy 
attenuation, illustrated by the cyclic wave s in Figure 1. 

Typically, a digit such as a finger, en ear lobe, or other portion of the body where blood flows 
dose to the sKn. is employed as the medium through which light energy is transmitted for blood gas 

30 attenuation measurements. The finger comprises skin, fat bone, muscle, etc. shown schematically in 
Figure 2, each of which attenuates energy incident on the finger in a generally predictable and constant 
manner. However, when fleshy portions of the finger are compressed erraticaly. for example by motion 
of the finger, energy attenuation becomes erratic 
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An example of a more realistic measured waveform S is shown in Figure 3. illustrating the effect 
of motion. The primary plethysmograprdc waveform portion of the signal s is the waveform representative 
of the pulse, corresponding to the sawtooth-Tike pattern wave in Rgure 1. The large, secondary motion- 
induced excursions m signal amplitude obscure the primary plethysmography signal s. Even small variations 
5 in amplitude make it difficult to distinguish the primary signal component s in the presence of a secondary 
signal component n. 

A pulse oximeter is a type of Wood gas monitor which non-invashiely measures the arterial 
saturation of oxygen in the blood. The pumping of the heart forces freshly oxygenated blood into the 
arteries causing greater energy attenuation. As weB understood in the art, the arterial saturation of 
10 oxygenated blood may be determined from the depth of the valleys relative to the peaks of two 
plethysmograprdc waveforms measured at separate wavelengths. Patient movement introduces motion 
artifacts to the composite signal as illustrated in the plethysmography waveform illustrated in Figure 3. 
These motion artifacts distort the measured signal. 

Summary of t he invention 

15 This invention provides improvements upon the methods and apparatus disclosed in U.S. Patent 

Application No. 08/132.812, fHed October 8, 1993, entitled Signal Processing Apparatus, which earDer 
application has been assigned to the assignee of the instant application. The present invention involves 
several different embodiments using the novel signal model in accordance with the present invention to 
isolate either e primary signal portion or a secondary signal portion of a composite measured signal. In one 

20 embodiment, a signal processor acquires a first measured signal and a second measured signal that is 
correlated to the first measured signal. The first signal comprises a first primary signal portion and a first 
secondary signal portion. The second signal comprises a second primary signal portion and a second 
secondary signal portion. Tha signals may be acquired by propagating energy through a medium and 
measuring an attenuated signal after transmission or reflection. AJteraatively, the signals may be acquired 

25 by measuring energy generated by the medhim. 

In one embodiment, the first and second measured signals are processed to generate a secondary 
reference which does not contain the primary signal portions from either of the first or second measured 
signals. This secondary reference is correlated to the secondary signal portion of eech of the first and 
second measured signals. The secondary reference is used to remove the secondary portion of each of the 

30 first and second measured signals via a correlation canceter, such as an adaptive noise cancefer. The 
correlation canceter is a device which takes a first and second input and removes from the first input aH 
signal components which are correlated to the second input. Any unh which performs or nearly performs 
this function is herein considered to be a correlation canceter. 
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An adaptive correlation canceler can be described by analogy to a dynamic multiple notch filter 
which dynamically changea its transfer function in response to a reference signel and the measured signals 
to remove frequencies from the measured signals that are also present in the reference signal. Thus, a 
typical adaptive correlation canceler receives the signal from which it is desired to remove a component 
5 and receives a reference signal of the undesired portion. The output of the correlation canceler is a good 
approximation to the desired signal with the undesired component removed. 

Alternatively, the first and second measured signals may be processed to generate a primary 
reference which does not contain the secondary signal portions from either of the first or second measured 
signals. The primary reference may then be used to remove the primary portion of each of the first and 
10 second measured signals via a correlation canceler. The output of the correction canceler is a good 
approximation to the secondary signal with the primary signal removed and may be used for subsequent 
processing in the same instrument or an auxiliary instrument. In this capacity, the approximation to the 
secondary signal may be used as a reference signal for input to a second correlation canceler together with 
either the first or second measured signals for computation of, respectively, either the first or second 

15 primary signal portions. 

Physiological monitors can benefit from signal processors of the present invention. Often in 
physiological meesurements a first signal comprising a first primary portion and a first secondary portion 
and a second signal comprising a second primary portion and a second secondary portion are acquired. The 
signals may be acquired by propagating energy through a patient's body (or a material which is derived from 

20 the body, such as breath, blood, or tissue, for example) or inside a vessel and measuring an attenuated 
signal after trensmission or reflection. Alternatively, the signal may be acquired by measuring energy 
genereted by a petienfa body, such es in electrocardiography. The signals are processed via the signal 
processor of the present invention to acquire either a secondary reference or a primary reference which 
is input to a correlation canceler. such as an adaptive noise canceler. 

25 One physiological monitoring apparatus which benefits from the present invention is a monitoring 

system which determines a signal which is representative of the arterial puke, called a plethysmography 
wave This signal can be used in blood pressure calculations, blood constituent measurements, etc A 
specific example of such a use is in pulse oxmnmy. Pulse oximetry involves determining the saturation of 
oxygen m the blood. In this configuration, the primary portion of the signal is the arterial blood contribufon 
30 to attenuation of energy as it passes through a portion of the body where blood flows dose to the skin. 
The pumping of the heart causes Wood flow to increase and decrease in the arteries in a period* f ashon. 
casing periodic attenuation wherein the periodic waveform is the pkrthysmographic waveform 
representative of the arterial pulse. The secondary portion is noise. In accordance with the present 
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invention, the measured signals are modeled such that this secondary portion of the signal is related to the 
venous blood contribution to attenuation of energy as it passes through the body. The secondary portion 
also includes artifacts due to patient movement which causes the venous blood to flow m an unpredictable 
manner, causing unpredictable attenuation and corrupting the otherwise periodic plethysmographic waveform. 
S Respiration also causes the secondary or noise portion to vary, although typically at a lower frequency than 
the patients pulse rate. Accordingly, the measured signal which forms a plethysmogrephic waveform is 
modeled in accordance with the present invention such that the primary portion of the signal is 
representative of arterial blood contribution to attenuation and the secondary portion is due to several other 
parameters. 

10 A physiological monitor particularly adapted to pulse oximetry oxygen saturation measurement 

comprises two fight emitting diodes (LEO's) which emit light at different wavelengths to produce first and 
second signals. A detector registers the attenuation of the two different energy signals after each passes 
through an absorptive media, for example a digit such as a finger, or an earlobe. The attenuated signals 
generally comprise both primary (arterial attenuator) and secondary (noise) signal portions. A static filtering 

15 system, such as a bandpass filter, removes a portion of the secondary signal which is outside of a known 
bandwidth of interest, leaving an erratic or random secondary signal portion, often caused by motion and 
often difficult to remove, along with the primary signal portion. 

A processor in accordance with one embodiment of the present invention removes the primary 
signal portions from the measured signals yielding a secondary reference which is a combination of the 

20 remaining secondary signal portions. The secondary reference is correlated to both of the secondary signal 
portions: The secondary reference and at least one of the measured signals are input to a correlation 
canceler, such as an adaptive noise canceler, which removes the random or erratic portion of the secondary 
signal. This yields a good approximation to a primary plethysmographic signal as measured at one of the 
measured signal wavelengths. As is known in the art, quantitative measurements of the emount of 

25 oxygenated arterial blood in the body can be determined from the ptethysmonraphic signal in a variety of 
ways. 

The processor of the present invention may also remove the secondary signal portions from the 
measured signals yielding a primary reference which is a combination of the remaining primary signal 
portions. The primary reference is correlated to both of the primary signal portions. The primary reference 
30 and at least one of the measured signals are input to a correlation canceler which removes the primary 
portions of the measured signals. This yields a good approximation to the secondary signal at one of the 
measured signal wavelengths. This signal may be useful for removing secondary signals from an auxfliary 
instrument as weB as determining venous blood oxygen saturation. 
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In accordance with the signal model of the present invention, the two measured signals each hawing 
primary and secondary signal portions can be related by coefficients. By relating the two equations with 
respect to coefficients defined in accordance with the present invention, the coefficients provide information 
about the arterial oxygen seturation end about the noise (the venous oxygen saturation and other 
5 parameters). In accordance with this aspect of the present invention, the coefficients can be determined 
by rninimizing the correlation between the primary and secondary signal portions as defined in the model 
Accordingly, the signal model of the present invention can be utilized in many ways in order to obtain 
information about the measured signals as wil be further apparent in the detailed description of the 
preferred embodiments. 

10 One aspect of the present invention is e method for usb in a signal processor in e signal processor 

for processing at least two measured signals S 1 and Sj each containing a primary signal portion s and a 
secondary signal portion n. the signals S, and S 2 being in accordance with the following relationship: 

51 = $1 + n i 

5 2 = s 2 + n z 

where s, and s 2 , and nj and n 2 are related by. 

$1 = r a s 2 and =r/7 2 

end where r, and r w are coefficients. 
15 The method comprises a number of steps. A value of coefficient r a is determined which minimize 

correction between s, and n,. Then, at least one of the first and second signals is processed using the 
determined value for r, to significantly reduce n from at toast one of the first or second measured signal 
to form a dean signal. 

In one embodiment the cleen signal is displayed on a display. In another embodiment, wherein the 
20 first and second signals are physiological signals, the method further comprises the step of processing the 
clean signal to determine a physiological parameter from the first or second measured signals. In one 
embodiment, the parameter is arterial oxygen saturation. In another embodiment the parameter is an ECG 
signal, in yet another embodiment, wherein the first portion of the measured signals is indicative of a heart 
plettrysmograph, the method further comprises the step of calculating the pulse rate. 
25 Another aspect of the present invention involves a physiological monitor. The monitor has a first 

input configured to receive a first measured signal S, having a primary portion. s v and a secondary portion 
n v The monitor also has a second input configured to received a second measured signal S 2 having a 
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primary portion s 2 and a secondary portion n 2 . Advantageously, the first and the second measured signals 
Sj and S 2 are in accordance with the following relationship: 

£, = s t + 
S 2 = Sz + n 2 

where Sj and s 2 , and nj and n 2 are related by: 

s t = r a s 2 and 

and where r 8 and r y are coefficients. 
5 The monitor further has a scan reference processor, the scan reference processor responds to a 

plurality of possible values for r a to multiply the second measured signal by each of the possible values for 
r and for each of the resulting values, to subtract the resulting values from the first measured signal to 
provide a plurality of output signals. A correlation canceler having a first input configured to receive the 
first measured signal, and having a second input configured to receive the plurality of output signals from 

10 the saturation scan reference processor, provides a plurality of output vectors corresponding to the 
correlation cancellation between the plurality of output signals and the first measured signal. An integrator 
having an input configured to receive the plurality of output vectors from the correlation canceler is 
responsive to the plurality of output vectors to determine a corresponding power for each output vector. 
An extremum detector is coupled at its input to the output of the integrator. The extremum detector is 

15 responsive to the corresponding power for each output vector to detect a selected power. 

In one embodiment, the plurality of possible values correspond to a plurality of possible values for 
a selected blood constituent In one embodiment the. the selected blood constituent is arterial blood oxygen 
saturation. In another embodiment the selected blood constituent is venous blood oxygen saturation. In 
yet another embodiment the selected blood constituent is carbon monoxide. 

20 Another aspect of the present invention involves a physiological monitor. The monitor has a first 

input configured to receive a first measured signal S, having a primary portion, s v and a secondary portion, 
n,. The monitor also has a second input configured to received a second measured signal S 2 having a 
primary portion s 2 and a secondary portion The first and the second measured signals S t and S 2 are 
in accordance with the following relationship: 

25 where Sj and s 2 , and n t and n 2 are related by: 
and where r a and r v are coefficients. 
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n 2 


and 


"1 



A transform module is responsive to the first and the second measured signals and responsive to 
a plurality of possible values for r, to provide at least one power curve as an output. An extremum 
calculation module is responsive to the at least one power curve to select a value for r a which minimizes 
the correlation between s and n. and to calculate from the value for r, a corresponding saturation value 
5 as an output A display module is responsive to the output of saturation calculation to display the 
saturation value. 

Rriof n»«rri ptinn nf the Drawirms 
Figure 1 illustrates an ideal plethysmography waveform. 
Figure 2 schematically illustrates a typical finger. 
10 Figure 3 illustrates a plethysmography waveform which includes a motion-induced erratic signal 

portion. 

figure 4a illustrates a schematic diagram of a physiological monitor to compute primary 
physiological signals. 

Figure 4b illustrates a schematic diagram of a physiological monitor to compute secondary signals. 
15 Figure Sa illustrates an example of an adaptive noise canceler which could be employed in a 

physiological monitor, to compute primary physiological signals. 

Figure 5b illustrates an example of an adaptive noise canceler which could be employed in a 
physiological monitor, to compute secondary motion artifact signals. 

Figure 5c illustrates the transf er function of a multiple notch filter. 
20 Rgure 6a illustrates a schematic of absorbing, material comprising N constituents within the 

absorbing material. 

figure 6b llustrates another schematic of absorbing material comprising N constituents, including 

one mixed layer, within the absorbing material 

figure 6c illustrates another schematic of absorbing material comprising N constituents, including 

25 two mixed layers, within the absorbing material 

figure 7a illustrates a schematic diagram of a monitor, to compute primary and secondary signs* 
in accordance with one aspect of the present invention. 



PCT/US95/13469 

WO 96/12435 



30 



-9- 



Figure 7b illustrates the ideal correlation canceler energy or power output as a function of the 
signal coefficients r v r 2 . ~r„. In this particular example, r 3 - r, and r 7 - t r 

Figure 7c illustrates the non ideal correlation canceler energy or power output as a function of the 
signal coefficients r,. r 2 . -r„. In this particular example. r 3 - r„ and r 7 - r w . 
S Figure 8 is a schematic model of a joint process estimator comprising a least-squares lattice 

predictor and a regression filter. 

Figure 8a is a schematic model of a joint process estimator comprising a QRO least-squares lattice 

{LSU predictor and a regression filter. 

Figure 9 is a flowchart representing a subroutine for implementing in software a joint process 

10 estimator as modeled in Figure 8. 

Figure 9a is a flowchart representing a subroutine for implementing in software a joint process 

estimator as modeled in Figure 8a. 

Figure 10 is a schematic model of a joint process estimator with a least-squares lattice predictor 

and two regression filters. 

15 Figure 10a is a schematic model of a joint process estimator with a QRO least-squares latt.ce 

predictor and two regression filters. 

Figure 11 is an example of a physiological monitor in accordance with the teachings of one aspect 

of the present invention. 

Rgure 1 1a illustrates an example of a low noise emitter current driver with accompanying digital 

20 to analog converter. 

Rgure 12 illustrates the front end analog signal conditioning circuitry and the analog to digital 

conversion circuitry of the physiological monitor of Rgure 11. 

Figure 13 illustrates further detail of the digital signal processing circuitry of Figure 11. 
Figure 14 illustrates additional detail of the operations performed by the digital signal processing 

25 circuitry of Rgure 11. 

Rgure 15 fllustrates additional detail regarding.the demodulation module of Rgure 14. 

Figure I B illustrates additional detail regarding the decimation module of Rgure 14. 

Rgure 17 represents a more detailed block diagram of the operations of the statistics module of 

Fgure 18 illustrates a block diagram of the operations of one embodiment of the saturation 

transform module of Figure 14. 

figure 19 illustrates a block diagram of the operation of the saturation calculation module of Rgure 



Rgure 14. 



14. 
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Figure 20 illustrates a block diagram of the operations of the pulse rate calculation module of 
Figure 14. 

Figure 21 illustrates e block diagram of the operations of the motion artifact suppression module 
of Figure 20. 

5 Figure 21a illustrates an alternative block diagram for the operations of the motion artifact 

suppression module of Figure 20. 

Figure 22 illustrates a saturation transform curve in accordance with the principles of the present 

invention. 

Figure 23 illustrates a block diagram of an alternative embodiment to the saturation transform in 
10 order to obtain a saturation value. 

Figure 24 illustrates e histogram saturation transform in accordance with the alternative 
embodiment of Figure 23. 

Figures 25A-25C illustrate yet enother alternative embodiment in order to obtain the saturation. 
Figure 26 Hlustrates a signal measured at a red wavelength yla->Ued-660nm for use in a 
15 processor of the present invention for determining the secondary reference n'(t) or the primary reference 
s'(t) and for use in a correlation canceler. The measured signal comprises a primary portion s^ a {t) and a 
secondary portion n^t). 

Figure 27 illustrates a signal measured at an infrared wavelength >ib->i| R -910nm for use in a 
processor of the present invention for determining the secondary reference n'tt) or the primary reference 
20 s'(t) and for use in a correlation canceler. The measured signal comprises a primary portion s^ft) and a 

secondary portion n^ b (t). 

figure 28 illustrates the secondary reference n'M determined by a processor of the present 

invention. 

Figure 29 Hlustrates a good epproximation s^ t (t) to the primary portion s^ a (t) of the signal S^t) 
25 measured at ^a->lred-660nm estimated by correlation cancellation with a secondary reference n'(t). 

figure 30 illustrates a good approximation s u Ab [X) to the primary portion z Ab i\) of the signal S Ab [t) 
measured at >lb-yUR-910nm estimated by correlation cancellation with a secondary reference n'(t). 

Figure 31 depicts a set of 3 concentric electrodes, Le^ a tripolar electrode sensor, to derive 
electrocardiography (ECGJ signals, denoted as S v S 2 and S 3 , for use with the present invention. Each of 
30 the ECS signals contains a primary portion and a secondary portion. 
Figure 32 shows a tapped delay fine FIR filter. 

Figure 33 shows phasor summation of FIR filter output for two frequencies, 
figure 34 shows an M-th order ad-pass subfOter. 
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Figure 35 shows complementary lowpass and highpass filters from a 2-path all-pass network. 
Figure 36 shows a Hilbert transform aU-pass filter reaction as a transformed half band filter. 
Figure 37 shows a resampling 2-path all-pass filter. 

Figure 38 shows the spectrum of an iterated 2-path filter with polynomials in Z^. 

Figure 39 shows the spectra of four channel quadrant centered and cardinal centered filter banks. 

Figure 40 shows a lowpass transformation of an aO pass filter. 

Rgure 41 A shows an aH-pass filter structure. 

Rgure 41 B shows a polyphase aH-pass structure. 

Figure 42 shows the pole-zero rfotribution for M-th order aH-pass subfilters. 

Rgures 43A-430 show typical phase responses of all-pass filter paths, and corresponding magnitude 

responses. ■ , _ . 

Figures 44A and 44B show the pole-zero plot and magnitude response for 2-path and 5-path 

polyphase all-pass networks. 

Rgure 45 shows a bandpass transformation of an aH-pass filter. 

notailpri Descr int'"" n f thB i™***™ 

The present invention involves a system which utilizes first and second measured signals that each 
contain a primary signal portion and a secondary sign,, portion, .n other words, given a first and second 
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composite signals S^t) - s,(t) ♦ n^t) and S 2 (t) - s 2 (t) ♦ n 2 tt>, the system of the presEnt invention can 
be used to isolate either the primary signal portion s(t) or the secondary signal portion nft). Followmfl 
processing, the output of the system provides a good approximation n"(t) to the secondary signal portion 
n(t) or a good approximation s'(t) to the primary signal portion aft). 

5 The system of the present invention is particularly useful where the primary andjor secondary signal 

portion n|t) may contain one or more of a constant portion, a predictable portion, an erratic portion, a 
random portion, etc The primary signal approximation s"(t) or secondary signal approximation n"(t) is 
derived by removing as many of the secondary signal portions n(t) or primary signal portions s(t) from the 
composite signal S(t) as possible. The remaining signal forms either the primary signal approximation s"(t) 

10 or secondary signal approximation n'(t), respectively. The constant portion and predictable portion of the 
secondary signal n(t) are easily removed with traditional filtering techniques, such as simple subtraction, 
low pass, band pass, and high pass filtering. The erratic portion is more difficult to remove due to its 
unpredictable nature. If something is known about the erratic signal even stetisticaDy, it could be removed, 
at least partially, from the measured signal via traditional filtering techniques. However, often no 

15 information is known about the erretic portion of the secondary signal nit). In this case, traditional filtering 

techniques are usually insufficient 

In order to remove the secondary signal n(t). a signal model in accordance with the present 
invention is defined as follows for the first and second measured signals S t and S 2 : 

5 1 = Si+n, 

5 2 = s 2 +n 2 

with s A = r a s 2 and = r v n 2 

or r a = — and = — 

a &, n 2 



20 



where s t and n, are at least somewhat (preferably substantially) uncorrected and s 2 and n 2 are at least 
somewhat (preferably substantially) uncorrelated. The first and second measured signals S, and S 2 are 



WO 96/12435 PCT/US95/13469 

•12- 

related by correlation coefficients r 8 and r v as defined above. The use and selection of these coefficients 
is described in further detail below. 

In accordance with one aspect of the present invention, this signal model is used in combination 
with a correlation canceler, such as an adaptive noise canceler, to remove or derive the erratic portion of 

5 the measured signals. 

Generally, a correlation canceler has two signal inputs and one output. One of the inputs is either 
the secondary reference n'lt) or the primary reference s*H which are correlated, respectively, to the 
secondary signal portions n(t> and the primary signal portions s(t| present in the composite signal S(t). The 
other input is for the composite signal S(t). Ideally, the output of the correlation canceler s"(t) or n"(t) 

10 corresponds, respectively, to the primary signal s(t) or the secondary signal n(t> portions only. Often, the 
most difficult task in the application of correlation cancelers is determining the reference signals n'(t) and 
s'(t) which are correlated to the secondary n(t) and primary s(t| portions, respectively, of the measured 
signal S(t) since, as discussed above, these portions are quite difficult to isolate from the measured signal 
S(t). In the signal processor of the present invention, either a secondary reference n'(t) or a primary 

15 reference s'{t) is determined from two composite signals measured simultaneously, or nearly simultaneously, 
at two different wavelengths. Aa end Ab. 

A block diagram of a generic monitor incorporating a signal processor according to the present 
invention, and a correlation canceler is shown in Figures 4a and 4b. Two measured signals, S^lt) and 
Sjufil are acquired by a detector 20. One skilled in the art will realize that for some physiological 

20 measurements, more than one detector may be advantageous. Each signal is conditioned by a signal 
conditioner 22a and 22b. Conditioning includes, but is not limited to, such procedures as filtering the 
signals to remove constant portions and amplifying the signals for ease of manipulation. The signals are 
then converted to digital data by an analog to-digital converter 24a and 24b. The first measured signal 
S^lt) comprises a first primary signal portion, labeled herein s,, a ltt. and a first secondary signal portion, 

25 lateled herein n^ltt. The second measured signal S Ab W is at least partially correlated to the first 
measured signal S^t) and comprises a second primary signal portion, labeled herein s Ab W. and a second 
secondary signal portion, labeled herein a Ab {t). Typically the first end second secondary signal portions, 
n>la (t, and n Ab H). are uncorrected and/or erratic with respect to the primery signal portions s Aa it) and 
s^(t). The secondary signal portions n^t) and n^tt) are often caused by motion of a patient in 

30 physiological measurements. 

The signals S^W and S Ab W are input to a reference processor 26. The reference processor 26 
multiplies the second measured signal S Ab W by either a factor r, - s^W/s^tt, or a factor r v - 
n^W/n^lt) and then subtracts the second measured signal S Ab M from the first measured signal S^lt). 
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The signal coefficient factors r a and r, are determined to cause either the primary signal portions s^t) 
and s„ b W or the secondary signal portions n At H) and n^ b (t) to cancel, respectively, when the two signals 
SylB ,„ and s^(t) are subtracted. Thus, the output of the reference processor 26 b either a secondary 
reference signal n'W - n^t) - in figure 4a. which is correleted to both of the secondery signal 

5 portions n^t) and n^ b (t) or a primary reference signal s'lt) - * Aa it) - r^ b W. in Figure 4b. which is 
correlated to both of the primary signal portions t At [\) and s^t). A reference signal n'W or s'(t) is input, 
along with one of the measured signals' S^t) or S^t). to a correlation cenceler 27 which uses the 
reference signal n'W or s'W to remove either the secondary signal portions n^t) or n^t) or the primary 
signal portions s^W or s^t) from the measured signal S^t) or S„ b W. The output of the correlation 
10 canceler 27 is a good primary signal approximation s'W or secondary signal approximation n'W- In one 
embodiment, the approximation s"W or n"W is displayed on a display 28. 

In one embodiment, an adaptive noise canceler 30. en example of which is shown in block diagram 
form in figure 5a. is employed es the correlation canceler 27. to remove either one of the erratic, 
secondary signal portions n,, 8 W and n^t, from the first and second signals S,, a W and S^t). The 
15 adaptive noise canceler 30 in figure Sa has as one input a sample of the secondary reference n'W which 
is correlated to the secondary signal portions n^W and n^t). The secondary reference n'W is determined 
from the two measured signals S^W and S^t) by the processor 26 of the present invention as described 
herein. A second input to the adaptive noise canceler. is a sample of either the first or second composite 
measured signals S^.W - ^.W ♦ n^W or S M M - s^W ♦ n^t). 
20 The adaptive noise canceler 30. in figure 5b. may also be employed to remove either one of 

primary signal portions s^W and s^t) from the first and second measured signals S^W and S^W. The 
adaptive noise cenceler 30 has as one input a sample of the primary reference s'W which is correlated to 
the primary signal portions s„ 8 W and s^t). The primary reference s'W is determined from the two 
measured signals S^t) and S^t) by the processor 26 of the present invention as described herem. A 
25 second input to the adaptive noise canceler 30 is a sample of either the first or second measured signals 

The adaptive noise canceler 30 functions to remove frequencies common to both the reference n'W 
or s'W and the measured s^nal S^.W or S^W. Since the reference signals are correlated to either the 
secondary signal portions n At [t) and or the primary signal portions s At [X) and $>lb W. the reference 
30 signals wiO be correspondingly erratic or weD behaved. The adaptive noise canceler 30 ects in a manner 
which may be analogized to a dynamic multiple notch filter based on the spectral distribution of the 
reference signal n'W or s'W- 
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Flgure 5c illustrates an exemplary transfer function of a multiple notch filter. The notches, or dips 
in the amplitude of the transfer function, indicate frequencies which are attenuated or removed when a 
signal passes through the notch filter. The output of the notch filter, is the composite signal having 
frequencies at which a notch is present removed. In thB analogy to an adaptive noise canceler 30, the 
5 frequencies at which notches are present change continuously based upon the inputs to the adaptive noise 
canceler 30. 

The adaptive noise canceler 30 "(figures 5a and 5b) produces an output signal, labeled herein as 
* m A* M ' *Ab®' n "jla w or ■%i W which '* fed D8ck t0 an internal P rocessor 32 within the adaptive noise 
canceler 30. The internal processor 32 automatically adjusts its own transfer function according to a 

10 predetermined algorithm such that the output of the internal processor 32 labeled b^tl in figure 5a and 
c A it) in figure 5b. closely resembles either the secondary signal portion n Aa ii) or n^t) or the primary signal 
portion Sy , a (t) or s,, D (t). The output b^t) of the internal processor 32 in figure 5a is subtracted from the 
measured signal, S M M or S^, b (t|. yielding a signal output s' Ab H) - s^ 8 W ♦ n^tt) - b^t) or a signal 
output s-^t) - s^ltKn^ltM^lt). The internal processor optimizes s^.lt) or s" Ab M such that s^W 

15 or s-^W is approximately equal to the primary signal s A Jt) or s Ah (t). respectively. The output c^t) of 
the internal processor 32 in figure 5b is subtracted from the measured signal, S^W or S^t), yielding a 
signal output given by *' At M - s^ 8 W ♦ n„ a (t) - c Aa [l) or a signal output given by n^W - s^tt) ♦ 
n^ B (t) - c^ b (t). The internal processor optimizes n^.W or n",, b (t) such that n",, a W or n' Ab M is 
approximately equal to the secondary signal portion n^ a (t) or n^lt), respectively. 

20 One algorithm which may be used for the adjustment of the transfer function of the internal 

processor 32 is a least-squares algorithm, as described in Chapter 6 and Chapter 12 of the book MSEliXS 
SinnM Processing bv Bernard Widrow and Samuel Stearns, published by Prentice Hall, copyright 1985. This 
entire book, including Chapters 6 and 12, is hereby incorporated herein by reference. 

Adaptive processors 30 in figures 5a and 5b have been successfully applied to a number of 

25 problems including antenna sidelobe canceling, pattern recognition, the elimination of periodic interference 
in general, end the elimination of echoes on long distance telephone transmission fines. However, 
considerable ingenuity is often required to find a suitable reference signal n'tt) or s'(t) since the portions 
■UW' I4.W and s^lt) cannot easily be separated from the measured composite signals S^M and 
S^lt). If eliher the actual secondary portion n^W or n^W or the primary signal portion s Aa [\) or s Ah <t) 

30 were a priori available, techniques such as correlation cancellation would not be necessary. 

RPMFBAII7ED DETFRMIWATIOM OF P RIMARY AND SECONDARY RFFERENCE SIGNA1 
An explanation which describes how the reference signals n'(t) and s'(t) may be determined follows. 
A first signal is measured at, for example, a wavelength As. by a detector yielding a signal S^lt): 
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(1) 



where s^ a (t) is the primary signal portion and n Aa M is the secondary signal portion. 

A similar measurement is taken simultaneously, or nearly simultaneously, at a different wavelength, 

5 Ab, yielding: 

*Ab® " *Ab® * n Ab®. l2) 

Note that as long as the measurements, S^W and S^ftL are taken substantially simultaneously, the 
10 secondary signal components, n^ a (t) and n^ftl. are correlated because any random or erratic functions 

affect each measurement m nearly the same fashion. The substantially predictable primary signal 

components. s At W and s^t), arB also correlated to one another. 

To obtain the reference signals n'lt) and s'ft). the measured signals S^ 8 (t) and S Ab i\) are 

transformed to eliminate, respectively, the primary or secondary signal components. In accordance with 
15 the present invention one way of doing this is to find proportionality constants, r, and r r between the 

primary signal portions s Aa i\) and s Ab W and the secondary signal portions n Aa M and n Ah [t) such that the 

signals can be modeled as follows: 

W« " 'a *Ab® 

(3) 

20 n A* {t) " r v n Ab ltl 

In accordance with the inventive signal model of the present invention, these proportionality relationships 
can be satisfied in many measurements, including but not limited to absorption measurements and 
physiological measurements. Additionally, in accordance with the signal model of the present invention, in 
most measurements, the proportionality constants r a and r, can be determined such that 

25 n A*® * r a ">lb W 

(4) 

Multiplying equation (2) by r, and then subtracting equation (2) from equation (1) results in a single 
equation wherein the primary signal terms s^ft) and s^tt) cancet 



30 



n'W - S^ffl • r.S^lt) - n^W - VV, D tt); <5a) 
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a non-zero signal which is correlated to each secondary signal portion n At it) and n Ab i\) and can be used 
as the secondary reference n'(t) in a correlation canceler such as an adaptive noise canceler. 

Multiplying equation (2) by r„ and then subtracting equation 12) from equation (1) results in a single 
equation wherein the secondary signal terms n^lt) and n^lt) caned, leaving: 



s'(t) - S^W • r„ S Ab H) - s Az it) • r ¥ s^lt); 



(5b) 



a non zero signal which is correlated to each of the primary signal portions s^t) and s^t) and can be 
used as the signal reference s'(t) in a correlation canceler such as an adaptive noise canceler. 

lfJ FXAMPLE OF DETERMINATION Q F PRIMARY AND SECONDARY 

REFERENCE SIGN&IS IN AN AB flDRPTIVE SYSTEM 
Correlation canceling is particularly useful in a large number of measurements generally described 
as absorption measurements. An example of an absorption type monitor which can advantageously employ 
correlation canceling, such as adaptive noise canceling, based upon a reference n'(t) or s'<t) determined by 

15 a processor of the present invention is one which determines the concentration of an energy absorbing 
constituent within an absorbing material when the material is subject to change. Such changes can be 
caused by forces about which information is desired or primary, or alternatively, by random or erratic 
secondary forces such as a mechanical force on the material. Random or erratic interference, such as 
motion, generates secondary components in the measured signal These secondary components can be 

20 removed or derived by the correlation canceler if a suitable secondary reference n'(t) or primary reference 
s'(t) is known. 

A schematic N constituent absorbing material comprising a container 42 having N different 
absorbing constituents, labeled A,, A 2 . A 3 ,.. A N , is shown in Hgur. 6a. The constituents A, through A N 
in Figure 6a are arranged in 8 generally orderly, layered fashion within the container 42. An example of 
25 a particular type of absorptive system is one m which fight energy passes through the container 42 and 
is absorbed according to tit generalized Beer-Lambert Law of light absorption. For fight of wavelength A*. 
this attenuation may be a[ proximated by: 



/ = / 0 exp(-f: € Ua cpc) m 

>1 
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InitiaDy transforming the signal by taking the natural logarithm of both sides and manipulating terms, the 
signal is transformed such that the signal components are combined by addition rather than multiplication. 
Le.: 

/=i 



5 where l n is the incident light energy intensity; I is the transmitted Gght energy intensity; e Wa is the 
absorption coefficient of the i ,h constituent at the wavelength Aa; Xjlt) is the optical path length of I 
layer, i.e., the thickness of material of the i th layer through which optical energy passes; and Cj(t) is the 
concentration of the i ,h constituent in the volume associated with the thickness Xj(t). The absorption 
coefficients e, through e H ere known values which are constant at each wavelength. Most concentrations 

10 CjW through c N (t) ere typically unknown, as are most of the optical path lengths Xj(t) of each layer. The 
total optical path length is the sum of each of the individual optical path lengths Xg(t) of each layer. 

When the material is not subject to any forces which cause change in the thicknesses of the 
layers, the optical path length of each layer. Xj{t), is generally constant This results in generally constant 
attenuation of the optical energy and thus, a generally constant offset in the measured signel. TypiceUy. 

15 this offset portion of the signal is of little interest since knowledge about a force which perturbs the 
material is usually desired. Any signal portion outside of a known bandwidth of interest including the 
constant undesired signal portion resulting from the generally constant absorption of the constituents when 
not subject to chenge, is removed. This is easily accomplished by traditional band pass filtering techniques. 
However, when the material is subject to forces, each layer of constituents may be affected by the 

20 perturbation differently than other layers. Some perturbations of the optical path lengths of each layer Xjjt) 
may result in excursions in the measured signal which represent desired or primary information. Other 
perturbations of the optical path length of each layer xj(t) cause undesired or secondary excursions which 
mask primary information in the measured signal. Secondary signal components associated with secondary 
excursions must also be removed to obtain primary information from the measured signal Similarly, the 

25 ability to compute secondary signal components caused by secondary excursions directly allows one to 
obtain primary signal components from the measured signal via simple subtraction, or correlation cancellation 
techniques. 

The correlation can eler may selectively remove from the composite signal meesured after being 
transmitted through or refuted from the absorbing material, either the secondary or the primeiy signal 
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components caused by forces which perturb or change the material differently from the forces which 
perturbed or charged the material to cause respectively, either the primary or secondary signal component. 
For the purposes of illustration, it will be assumed that the portion of the measured signal which is deemed 
to be the primary signal s*W is the attenuation term * 5 c 5 x 6 W associated with a constituent of interest 

5 namely A 5 . and that the layer of constituent A 5 is affected by perturbations different than each of the 
layers of other constituents A, through A 4 and A 6 through A„. An example of such a situation is when 
layer A 5 is subject to forces about which information is deemed to be primary and. additionally, the entire 
material is subject to forces which affect each of the layers. In this case, since the total force affecting 
the layer of constituent A 5 is different than the total forces affecting each of the other layers and 

10 information is deemed to be primary about the forces and resultant perturbation of the layer of constituent 
A 5 . attenuation terms due to constituents A t through A 4 and A 8 through A N make up the secondary signal 
portion n, a (t). Even if the additional forces which affect the entire material cause the same perturbation 
in each layer, including the layer of A 5 , the total forces on the layer of constituent A 5 cause it to have 
different total perturbation than each of the other layers of constituents A, through A 4 and A 6 through 

15 A N . 

It is often the case that the total perturbation affecting the layers associated with the secondary 
signal components is caused by random or erratic forces. This causes the thickness of layers to change 
erratically and the optical path length of each layer, xjW, to change erratically, thereby producing a random 
or erratic secondary signal component However, regardless of whether or not the secondary signal 

20 portion n^t) h erratic, the secondary signal component n^.lt) can be either removed or derived v.a a 
correlation canceler, such as an adaptive noise canceler, having as one input, respectively, a secondary 
reference n'W or a primary reference s'W determined by a processor of the present invention as long as 
the perturbation on layers other than the layer of constituent A 5 is different than the perturbation on the 
layer of constituent A 5 . The correlation canceler yields a good approximation to either the primary s.gnal 

25 s, W or the secondary signal n„ 8 <t). In the event that an approximation to the primary signal is obtained, 
the concentration of the constituent of interest, c 5 (tj. can often be determined since in some physical 
measurements, the thickness of the primery signal component, x 5 U) in this example, is known or can be 
determined. 

The correlation canceler utilizes either the secondary reference n'W or the primary reference s'W 
30 determined from two substantially simultaneously measured signals S„,W and S^lt). S>,W is determined 
as above in equation (7). S^Ct) is determined similarly at a different wavelength ,1b. To find either the 
secondary reference n'W or the primary reference s'W. attenuated transmitted energy is measured at the 
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two different wavelengths A* and Ab and transformed via logarithmic conversion. The signals S At {l) and 
S^M can then be written (logarithm converted) as: 

4 N 

S^«- f 5^a c 5» 5 W + ^ e Ua c i*i + £ *Wa¥i (8) 
M /=6 



(9) 



4 N 

WMlW* ^ f Wb c i«i + £ f Wb c i*i 1101 



/=1 /=6 



10 Further transformations of the signals are the proportionality relationships in accordance with the 

signal model of the present invention defining r, end r ¥> similar to equation (3). which allows determination 
of a noise reference n'(t) and a primary reference s'(t). These are: 

- 'A* n2 " 

15 where 

(13a) 

H. * (13W 
h is often the cai 3 that both equations (12) and (13) can be simultaneously satisfied. Multiplying equation 

(11) by r, and subtracting the result from equation (9) yields a non-zero secondary reference which is a 
20 linear sum of secondary signal components: 

nW - S^lt) - r a S^t) - - V^t) n4a) 



4 N 

- 53 ^eft-It)- II e M8 c^t). 
/=1 /=6 
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/=1 /=6 



(15a) 



/=1 



10 



A/ 
/=6 



Multiplying equation (11) by r„ and subtracting the result from equation (9) yields a primary 
reference which is a linear sum of primary signal components: 

s'lt) - S,, a (t) : r,S Ab i\) - t Aa W • r Ab W < 14l » 
-«= 5 x 5 (t)6 5i>la .r ¥ c 5 x 5 (t)e 5(y , b < 15 « 

15 A sample of either the secondary reference n'(t) or the primary reference s'lt), and a sample of 

either measured signal S^(t) or S^t). are input to a correlation canceler 27, such as an adaptive noise 
canceler 30. an example of which is shown m Figures 5a and 5b and a preferred example of which is 
discussed herein under the heading PREFERRED CORRELATION CANCELER USING A JOINT PROCESS 
ESTIMATOR IMPLEMENTATION. The correlation canceler 27 removes either the secondary portion n M [t) 

20 or n^t). or the primary portions, s Am M or s^t). of the measured signal yielding a good approximation to 
either the primary signals s^W = * Wa c 5 x 5 (t) or s'^t) * ^c 5 x 5 (t) or the secondary s,gnals 
n-^t, „ n ^ a (t) or □",„» = n,, D (t>. In the event that the primary signals are obtained, the concentration 
c 5 (t) may then be determined from the approximation to the primary signal s' Aa it) or s^tt) according to: 
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c 5 W « fjafltovaftfi im) 
As discussed previously, the absorption coefficients are constant at each wavelength Az and Ab and the 
thickness of the primary signal component. x s (t) in this example, is often known or can be determined as 
a function of time, thereby allowing calculation of the concentration c 5 (t) of constituent A s . 

fiFTFRMINATlON OF CPHnFMTRATlOU OR SATURATION 
IM A VOLUME CnWTAIMING MOP f ™AM ONE CONSTITUENT 

Referring to Figure 6b, another material having N different constituents arranged in layers is shown. 
In this material two constituents A 5 and A 6 are found within one layer having thickness x 5 6 (t) - x 5 (t) 
♦ x 6 (t). located generally randomly within the layer. This is analogous to combining the layers of 
constituents A 5 and A 8 in Figure 6a. A combination of layers, such as the combination of layers of 
constituents A s and A 8 . is feasible when the two layers are under the same total forces which result in 
the same change of the optical path lengths Xgft) and x 6 (t) of the layers. 

Often it is desirable to find the concentration or the saturation, i.e.. a percent concentration, of 
one constituent within a given thickness which contains more than one constituent and is subject to unique 
forces. A determination of the. concentration or the saturation of a constituent within a given volume may 
be made with any number of constituents in the volume subject to the same total forces and therefore 
under the same perturbation or change. To determine the saturation of one constituent in a volume 
comprising many constituents, as many measured signals as there are constituents which absorb incident 
fight energy are necessary. It will be understood that constituents which do not absorb fight energy are 
not consequential in the determination of saturation. To determine the concentration, as many signals as 
there are constituents which absorb incident light energy are necessary as well as information about the 
sum of concentrations. 

It is often the case that a thickness under unique motion contains only two constituents. For 
example, it may be desirable to know the concentretion or saturation of A 5 within a given volume which 
contains A 5 and A* In this case, the primary signals s^M and Sylb (t) comprise terms related to both A s 
and A 8 so that a determination of the concentration or saturation of A 5 or A 6 in the volume may be made. 
A determination of saturation is discussed herein. It will be understood that the concentration of A 5 in 
a volume containing both A 5 and A 8 could also be determined if it is known that A 5 <-A 8 -1. Le- that there 
are no constituents in the volume which do not absorb incident light energy at the particular measurement 
wavelengths chosen. The measured signals S^t) and S Ab W can be written (logarithm converted) as: 
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■ 4 «m.wmP + (18al 

- W<> + H.P (18b) 

{19W 

It is also often the case that there may be two or more thicknesses within a medium each 
10 containing the same two constituents but each experiencing a separate motion as in Figure 6c. For 
example, it may be desirable to know the concentration or saturation of A 5 within a given volume which 
contains' A 5 and A 6 as well as the concentration or saturation of A 3 within a given volume which contains 
A 3 and A 4 . A 3 and A 4 having the same constituency as A s and A„ respectively. In this cese, the primary 
signals S/ , 8 (t) and i^W again comprise terms related to both A 5 and A 6 and portions of the secondary 
15 signals n^lt) and n^W comprise terms related to both A 3 and A 4 . The layers. A 3 and A 4 , do not enter 
into the primary equation because they are assumed to be perturbed by a different frequency, or random 
or erratic secondary forces which are uncorrected with the primary force. Since constituents 3 and 5 as 
well as constituents 4 and 6 are taken to be the same, they have the same absorption coefficients <i.e.. 

HM " W - «U. * 'i* and ' « 6ener8,,V SPe3kin9 ' h0W6Ver ' " 3 
20 A 4 wffl have different concentrations than A 5 and A 6 and will therefore have a different saturate 

Consequently a single constituent within a medium may have one or more saturations associated with it. 
The primary and secondary signals according to this model may be written as: 



25 



n Aa W " I f 5^a c 3 - e 6^a c 4l x 3>4 W 
2 N 

/=1 h7. 



(20a) 



(20b) 



n Aa® " l e 5M C 2 * € *M C * ] X 3.4«» + n A^ 



(20c) 
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(21a) 



2 N 

* E.€ Wb C#t)* E CybCiXjtt). 

/=1 /=7 



(21b) 



15 



20 



tl.W " ^5^»b c 3 + e 6^b c 4l *3A® * ^b^ 



(21c) 



where signals n At [l) and ^ b (t) are similar to the secondary signals n^lt) and n^t) except for the 

10 omission of the 3, 4 layer. 

Any signal portions whether primary or secondary, outside of a known bandwidth of interest, 
including the constant undesired secondary signal portion resulting from the generally constant absorption 
of the constituents when not under perturbation, should be removed to determine an approximation to either 
the primary signal or the secondary signal within the bandwidth of interest. This is easily accomplished 
by traditional band pass filtering techniques. As in the previous example, it is often the case that the total 
perturbation or change affecting the layers associated with the secondary signal components is caused by 
random or erratic forces, causing the thickness of each layer, or the optical path length of each layer. Xj(t). 
to change erratically, producing a random or erratic secondary signal component n^.ft). Regardless of 
whether or not the secondary signal portion n ila (t) is erratic, the secondary signal component n^.lt) can 
be removed or derived via a correlation canceler. such as an adaptive noise canceler. having as one input 
a secondary reference n'(t) or a primary reference s'(t) determined by a processor of the present invention 
as long as the perturbation in layers other than the layer of constituents Ag and A 6 is different than the 
perturbation in the layer of constituents Ag and A 6 . Ether the erratic secondary signal components n^W 
and n>lb (t) or the primary components s^ a (t) and s Ab (t) may advantageously be removed from equations (18) 
25 and (19). or alternatively equations (20) and (21). by a correlation canceler. The correlation canceler, again, 
requires a sample of either the primary reference s'(t) or the secondary reference n'(t) and a sample of 
either of the composite signals S^tt) or S Ab it) of equations (18) and (19). 



10 
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nPT f RMIMATIQM tlF PRIMARY AN p Sf MNOARY RFFFRENCE S1GNA|.S 
FOR SATURATin M MEASUREMENTS 
One method for determining reference signals s'(t) or n'(t) from the measured signals Sj, § W and 
S A «t) in accordance with one aspect of the invention is what wfll be referred to as the constant saturation 
approach. In this approach, it is assumed that the saturation of Ag in the volume containing A 5 end A 6 
and the saturation of A 3 in the volume containing A 3 end A 4 remains relatively constent over some period 
of time. Le.: 

Saturation(A 5 |t» - c 5 (t) / (c 5 (t) ♦ c 6 (t)] (22a) 

Saturation(A 3 (t)) - c 3 {t) / [c 3 (t) ♦ c 4 (t)] »22b) 

Saturation(A 5 (t» - {1 ♦ (c 6 (t)/c 5 (t)l } 1 1™ 

15 Saturation(A 3 (t» - {1 + lc 4 (t)/c 3 (tH}- 1 < 23 °J 

are substantially constant over many samples of the measured signals S^ a and S^. This assumption is 
accurate over many samples since saturation generally changes relatively slowly in physiological systems. 
The constant saturation assumption is equivalent to assuming that 

20 

c 5 (t)/cg(t) - constant , <24a) 

Cgttl/c^t) - constant 2 (24b) 

25 since the only other term in equations (23a) and (23b) is a constant, namely the numeral 1. 

Using this assumption, the proportionality constants r a and r„ which allow determination of the 
secondary reference signal n'(t) and the primary reference signal s'(t) in the constant saturation method are: 

HM C 5 »S,6ft > * € 6M c 6 x 5.6tt> (25a) 

30 r, - — 

Hjlb c 5 X 5,6 W * *6^b ^ X 5.6 W 

(26a) 
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(C 5 ,C 6 ) * € 6^b 
« s^W/s^tt) - constant 3 ; where 

and 



20 



30 



15 e 5^b c 3 x 3.4 W) + e 6Mb c 4*3.4< tn 



25 e itifofeJ + Cs^b 

* n" /4 ,W/n" ylb (t) - constant*- where 



(27a) 



f 5^la « c 5' c 6^ * e *M (28a) 



(29a) 
(30a) 



(25b) 



(26b) 



(27b) 



(28b) 

(29b) 
(30b) 



In accordance with the present invention, it is often the case that both equations (26) and (30) 
be simultaneously satisfied to deterrnine the proportionality constants r a and r r Additionaly. the 
absorption coefficients at each wavelength e BM e 6M c^r and € BJk are constant and the centra! 
assumption of the constant saturation method is that c 5 (t)fc 6 (t) and c 3 (t)/c 4 (t) are constant over many 



can 



PCTYUS95/13469 

WO 96/12435 

-26- 

sample periods. Thus, new proportionality constants r. and r v may be determined every few samples from 
new approximations to either the primary or secondary signal as output from the correlation canceler. 
Thus, the approximations to either the primary signals z Aa {t) and s^tt) or the secondary signals n^W and 
n (t) found by the correlation canceler for a substantially immediately preceding set of samples of the 
measured signals S,.(t) and S^t) are used in a processor of the present invention for calculating the 
proportionality constants, r„ and r„ for the next set of samples of the measured signals S, a (t) and S, b (t). 

Multiplying equation (19) by r a and subtracting the resulting equation from equation (18) yields a 
non-zero secondary reference signal: 



10 



n'(t)-S^ a (t)-r a S > , b lt)-n i|a (t)-r a n >lb (t). 



(31a) 



Multiplying equation (19) by r„ and subtracting the resulting equation from equation (18) yields a 
non-zero primary reference signal: 



(31b) 



15 s'(t)-S y , a (t)-r w S >lb (t)-s y|a (t)-r v s >lD (t). 

When using the constant saturation method in patient monitoring, initial proportionality coefficients 
can be determined as further explained below. It is not necessary for the patient to remain motionless even 
for an initialization period. With values for the proportionality coefficients r 8 and r„ determmed. a 
20 correlation canceler may be utilized with a secondary reference n'(t) or a primary reference s'(t). 

D ETERMINATION nf RIRMAl CQEF FIH c "T g *W pRIMARY AND 
CErnainARV REFEREE SIRNALS U S IN fi THf CONS TANT SATURATION MR 
In accordance with one aspect of the present invention, the reference processor 26 of Figures 4a 
and figure 4b of the present invention may be configured to multiply the second measured assumed s.gnal 
25 S .KM-s^O^lt) by each of a plurality of gigna. coefficients r v r 2 . -r„ and then subtract each result 
from the first measured signal S^W-s^n^t) to.obtain a plurality of reference signals 

R'(r. t) - s,, a (t) - r s Ab i\) ♦ "n^W • r n^t) < 32 > 

30 for r - r,. r 2 . -r B as shown in Figure 7a. In other words, a plurality of signal coefficients are chosen 
to represent a cross section of possible signal coefficients. 

,„ order to determine either the primary reference s'(t) or the secondary reference n'(t) from the 
above plurality of reference signa.s of equation (32). signal coefficients r, and r„ are determined from the 
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; !urafity of assumed signal coefficients r v r 2 , - T n . The coefficients r, and r„ are selected such that they 
cause either the primary signal portions s^lt) and s^t) or the secondary signal portions n^ t M and a M \t) 
to cancel or nearly cancel when they are substituted into the reference function R'(r. t). e. g. 

M.W-'vKbM (33W 
n'(t) - R'(r r t) - n^ a (t) - ^n^^t) (33c) 

s'it) - R'(r r t) - s^lt) - r^tt). (33d) 



In other words, coefficients r, and r„ are selected at values which reflect the rninhnum of 
correlation between the primary signal portions and the secondary signal portions. In practice, one does 

15 not usually have significant prior information about either tha primary signal portions s Aa [t) and s Ab H) or 
the secondary signal portions n jta (t) and n^t) of the measured signals S^t) and S^ b (t). The lack of this 
information makes it difficult to determine which of the plurality of coefficients r,, r 2 . ~r„ correspond to 
the signal coefficients r a - s^ a (t)/s^ D (t) and r„ - n^tMn^ft). 

One approach to determine the signal coefficients r a and r„ from the plurality of coefficients r v 

20 r 2 , ~r„ employs the use of a correlation canceler 27, such as an adaptive noise canceler, which takes a 
first input which corresponds to one of the measured signals S^lt) or S^ b (t) and takes a second input 
which corresponds to successively each one of the pluraBty of reference signals R'(r v t), R'(r 2 , t), — -R*(r„» 
t) as shown in Figure 7a. For each of the reference signals R'lr,, t). R'(r 2 , t), - t) the corresponding 
output of the correlation canceler 27 is input to a "squares" operation 28 which squares the output of the 

25 correlation canceler 27. The output of the squares operation 28 is provided to an integrator 29 for forming 
a cumulative output signal (a summation of the squares).. The cumulative output signal is subsequently 
input to an extremum detector 31. The purpose of the extremum detector 31 is to chose signal 
coefficients r a and r v from the set r v r 2 . • r tt by observing which provide a maximum in the cumulative 
output signal as in Figures 7b and 7c In other words, coefficients which provide a maximum integrated 

30 output such as energy or power, from the correlation canceler 27 correspond to the signal coefficients r a 
end r ¥ which relate to a minimum correlation between the primary signal portions and the secondary signal 
portions in accordance with the signal model of the present invention. One could also configure a system 
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geometry which would require one to locate the coefficients from the set r v r 2 . -t n which provide e 
minimum or inflection in the cumulative output signal to identify the signal coefficients r, end r r 

Use of a plurality of coefficients in the processor of the present invention in conjunction with a 
correlation canceler 27 to determine the signal coefficients r, and r„ may be demonstreted by using the 
properties of correlation cancellation. If x, y and z are taken to be eny collection of three time varymg 
signals, then the properties of some correlation cancelers C(x. y) may be defined as follows: 

Property (1) Ox, y) - 0 for x. y correlated (34a) 
Property (2) C(x. y) - x for x. y uncorrected (34b) 
Property (3) C(x + y. z) - C(x. z) + C(y. z) (34c) 

With properties (1). (2) and (3) it is easy to demonstrate that the energy or power output of . correlation 
canceler with a first input which corresponds to one of the measured signals S A Jl) or S^lt) and a second 
input which corresponds to successively each one of a plurality of reference signals R'(r v t), R'(r 2 . t), •-. 
R',r t) can determine the signal coefficients r 8 and r„ needed to produce the primary reference s'(t) end 
secondary reference n'(t). If we take es a first input to the correlation canceler the measured signal S Aa W 
and as a second input the plurality of reference signals R*(r v t). R'(r 2 . t). ~.R'(r n . t) then the outputs of 
the correlation canceler CIS^U). R'^t)) for j - 1. 2, n may be written as 



25 



20 Cfyjlt) ♦ n^s^t)- r^tt) ♦ n^t)- r^lt)) 

where j - 1, 2. -,n and we have used the expressions 

R'(r. t) - S^lt) - xS Ah W 

- s As it) * n An M 

- s Ah [t) + n^t). 

30 The use of property (3) allows one to expand equetion (35) into two terms 

C|S y , a (t).R , (r.t))-C(s /la (t),s > , 8 (t).rs i i b (t)+n / , a (t).m >tb (tH 



(35) 



(36) 
(37a) 
(37b) 
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+c ^aW• s >la^ t )• rs >^bW +n >^8^ t > ^n >^bW , 
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(38) 



SO 



> that upon use of properties (1) and (2) the correlation canceler output is given by 
C[S A9 [\), Wityt))- s^tj^r.) ♦ ^ a (t)tf(r r g 



(39) 



5 



where 6[x) is the unit impulse function 



6M - 0 if x * 0 



(40) 
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<5(x) - 1 if X - 0. 



The time variable, t, of the correlation canceler output C[S At {t), R'frj, t)) may be eliminated by 
computing its energy or power. The energy of the correlation canceler output is given by 



It should be understood that one could, equally well, have chosen the measured signal S Ab it) as 
the first input to the correlation canceler and the plurality of reference signals R'lr,, t). R'(r 2 . t). -.Kir v 
20 t) as the second input. In this event, the correlation canceler energy output is 

E^frjl- J C 2 (Sj, b {t).R'{r,t)dt < 41W 
-<5(r r r,l J s 2 Ah m*Sit^ \ rr^lOdt. 

25 It should also be understood that in practical situations the use of discrete time measurement 

signals may be employed as weO as continuous time measurement signals. A system which performs a 
discrete transform (e.g„ a saturation transform in the present example) in accordance with the present 
invention is described with reference to Figures 11-22. In the event that discrete time measurement signals 
are used/integration approximation methods such as the trapezoid rule, midpoint rule, Tick's rule. Simpson's 

30 approximation or other techniques miry be used to compute the correlation canceler energy or power output 
In the dUcrete time measurement signal case, the energy output of the correlation canceler may be written, 
using the trapezoid rule, as 



15 



E^frj)- J C 2 (S / , a (t),R'(r j .t)dt 



(t)dt. (41a) 
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10 



20 



(42a) 



• 1 B /-0 



/» 

1 /=0 



(42b) 



where t, is the i* discrete time. t 0 is the initial time. t n is the final time and At is the time between 

discrete tone measurement samples. 

The energy functions given above, and shown in figure 7b. indicate that the correlation canceler 
output is usually zero due to correlation between tha measured signal S^ a (t) or SyUj (t) and many of the 
plurality of reference signals Mr,, t). R'(r 2 . t|, -jr|r r ft However, the energy functions ere non zero at 
values of r- which correspond to cancellation of either the primary signal portions s Ab W and s,, b (t) or the 
secondary 'signal portions n„(t) and n, b (t) in the reference signal *»> ft These veh.es correspond to the 

15 signal coefficients r a and r r 

It should be understood that there may be instances in time when either the primary signal portums 
», (t) end s. b (t) or the secondary signal portions n^t) and n, b (t) are identicaPy zero or nearly zero, to 
the'se cases, only one signal coefficient value wul provide maximum energy or power output of the 
correlation canceler. 

Since there may be more than one signal coefficient value which provides maximum correlate 
canceler energy or power output, an ambiguity may arise. It may not be immediately obvious which signal 
coefficient together with the reference function R'(r. t| provides either the primary or secondary reference. 
In such cases, it is necessary to consider the constraints of the physical system at hand. For example. 
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in puke oximetry, it is known that arterial blood, whose signature * the primery plethysmographic wave, 
has greater oxygen saturation than venous blood, whose signature is the secondary erratic or random signal 
Consequently, in pulse oximetry, the retio of the primary signals due to arterial pulsation r, - *U a (t)M iM p) 
is the smaller of the two signal coefficient values while the ratio of the secondary signals due to maHy 
5 venous blood dynamics r„ - n Aa l\)in Ab M is the larger of the two signal coefficient values, assuming >ia 

- 660 nm and Ab - 910 run. 

It should also be understood that'in practical uplementations of the plurality of reference signals 
and cross correlator technique, the ideal features listed as t operties (1). (2) and (3) above wul not be 
precisely satisfied but will be approximations thereof. Thq efore, in practical implementations of this 
10 embodiment of the present invention, the correlation canceler energy curves depicted in figure 7b wM not 
consist of infinitely narrow delta functions but will have finite width essociated with them as depicted in 
Figure 7c 

It should also be understood that it is possible to have more than two signal coefficient vehres 
which produce maximum energy or power output from a correlation cenceler. This situation arises when 
15 the measured signals each contain more than two components each of which are related by a ratio as 
follows: 



n 



(43) 
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n 

E 

/=o 



s^w- E f^o 



where 
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Thus, ref erence signal techniques together with a correlation cancellation, such as an adaptive noise 
canceler. can be employed to decompose a signal into two or more signal components each of which is 
related by a ratio. 

PRFFFRRED COF^'ATIOU CAI^ FIFR USING A 
■IfllNT PROCESS ESTIMATOR IMPLEMENTATION 
Once either the secondary reference n'(t) or the primary reference s'(t) is determined by the 
processor of the present invention, the correlation canceler can be implemented in either hardware or 
software: The preferred implementation of a correlation canceler is that of an adaptive noise canceler using 

a joint process estimator. 

The least mean squares (LMS) implementation of the internal processor 32 described above in 
conjunction with the adaptive noise canceler of Figures 5a and Figure 5b is relatively easy to implement, 
but lacks the speed of adaptation desirable for most physiological monitoring applications of the present 
invention Thus, a faster approach for adaptive noise canceling, called a least-souares lattice joint process 
estimator model, is used in one embodiment. A joint process estimator 60 is shown diagrammatically » 
figure 8 and is described in detail in Chapter 8 of Active Filter Theory by Simon Haykin. published by 
Prentice-Hall. copyright 1986. This entire book, including Chapter 9. is hereby incorporated herem by 

reference. ... 

The function of the joint process estimator is to remove either the secondary s.gnal porfons 
n. It) or n^lt) or the primary signal portions s^tl or s Ab [X) from the measured signals S^W or S Ab it). 
yielding either a primary signal approximation O a (t) or s^t) or a secondary signal approximate 
or n" , b lt) Thus, the joint process estimator estimates either the value of the primary signals s At it) or 
s, D (t) or the secondary signals n A jX) or n M iX). The inputs to the joint process estimator 60 are either 
the secondary reference n'(t) or the primary reference s'(t) and the composite measured signal S„ a M or 
S A 0L The output is a good approximation to the signal S, a (t) or S, D (t) with either the secondary s.gnal 
or the primary signal removed. Le. a good approximation to either s, a (t), $>lD W. n^M or n Ab iX). 

The joint process estimator 60 of Figure 8 utilizes, in conjunction, a least square lattice pred.ctor 
70 and a regression filter 80. Either the secondary reference n'W or the primary reference s'(t) is input 
,o the least square lattice predictor 70 while the measured signd S M W or S^t) is input to the 
regression filter 80. For simplicity in the following description, S^t) wfll be the measured signal from 
which either the primary portion Syla (t) or the secondary portion n, a (t) will be estimeted by the ,ourt 
process estimator 60. However. it wfll be noted that S^t) could also be input to the regression filter 80 
and the primary portion z Ab M or the secondary portion n^ltj of this signal could be estmiated. 
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The joint process estimator 60 removes ell frequencies that are present in both the reference n'tt) 
or s'tt), and the measured signal S^.tt). The secondary signal portion n M H) usually comprises frequencies 
unrelated to those of the primary signal portion s^.tt). It is improbable that the secondary signal portion 
•V| S W would be of exectly the same spectral content as the primary signal portion s,, a tt). However, in the 

S unBkely event that the spectral content of s At it) and n^.ft) are similar, this approach will not yield accurate 
results. Functionally, the joint process esthnetor 60 compares the reference input signal n'tt) or s'tt). which 
is correlated to either the secondary signal portion n^tt) or the primary signal portion s^tt), and input 
signal S^tt) and removes all frequencies which are identical. Thus, the joint process estimator 60 acts 
as a dynamic multiple notch filter to remove those frequencies in the secondary signal component n^ft) 

10 as they change erratically with the motion of the patient or those frequencies in the primery signal 
component s^tt) as they change with the arterial pulsation of the patient. This yields a signal having 
substantially the same spectral content and ampDtuda as either the primary signal s At i\) or the secondary 
signal n>|B (t). Thus, the output s' M l\) or n^.tt) of the joint process estimator 60 is a very good 
approximation to either the primery signal s^tt) or the secondary signal n^,(t). 

15 The joint process estimator 60 can be divided into stages, beginning with a zero-stage and 

terminating in an nAstage. as shown in Figure 8. Each stage, except for the zero-siege, is identical to 
every other stage. The zero stage is an input stage for the joint process estimator 60. The first stage 
through the nAstage work on the signal produced in the immediately previous stege, Le.. the (m-1) tn -stage, 
such that a good primary signal approximation s^.tt) or secondary signal approximation n"^ a (t) is produced 

20 as output from the m^-stage. 

The leestsquares lattice predictor 70 comprises registers 90 and 92, summing elements 100 and 
102. and delay elements 110. The registers 90 and 92 contain multiplicative values of a forward reflection 
coefficient r f<m (t| and a backward reflection coefficient r b Jt) which multiply the reference signal n (tl 
or s'(t) and signals derived from the reference signal n'tt) or s'ft). Each stage of the least-squares lattice 

25 predictor outputs a forward prediction error f m (t) and a backward prediction error b m (t). The subscript m 

is indicative of the stege. 

For each set of samples, ue. one sample of the reference signal n'tt) or s'tt) derived substentiafly 
simultaneously with one sample of the measured signal S^tt), the sample of the reference signal n'tt) 
or s'tt) is input to the least-squares lattice predictor 70. The zero-stage forward prediction error f 0 (t) and 
30 the zero-sj ge backward prediction error b 0 tt) are set equal to the reference signal n'tt) or s'(t). The 
backward , rediction error b„{t) is delayed by one sample period by the delay element 110 in the first stage 
of the fc * -squares lattice predictor 7a Thus, the immediately previous value of the reference n'tt) or s'tt) 
is used h calculations iwolving the first-stage delay element 110. The zero-stage forward prediction error 
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is added to the negative of the delayed zero-stage backward prediction error bntt-1, multiplied by the 
forward reflection coefficient value r f>1 (t) register 90 value, to produce a first-stage forward prediction 
error f , Ufc Additionally, the zero-stage forward prediction error f 0 (t) is multiplied by the backward reflection 
coefficient r„ } (\) register 92 value end added to the delayed zero-stage backward prediction error b 0 (M) 

5 to produce a first-stage backward prediction error b,(t). In each subsequent stage, m, of the least square 
lattice predictor 70, the previous forward and backward prediction error values, f^tt) and b^OU, the 
backward prediction error being delayed by one sample period, are used to produce values of the forward 
and backward prediction errors for the present stege, f m (t) and b m (t). 

The backward prediction error b m (t) is fed to the concurrent stege. m. of the regression filter 80. 

10 There rt is input to a register 96, which contains a multiplicative regression coefficient value ^.(t). For 
example, in the zero-stage of the regression filter 80, the zero-stage backwerd prediction error b 0 W & 
multiplied by the zero stage regression coefficient k 0M IX) register 96 value and subtracted from the 
measured value of the signal S^lt) at a summing element 106 to produce a first stage estimation error 
signal e Ua <t). The firststege estimation error signal e Ua (t) is a first approximation to either the primary 

IS signal or the secondary signal. This first-stage estimation error signal e Ua <t) is input to the first-stage 
of the regression filter 80. The first-stage backward prediction error b,(t). multiplied by the first-stage 
regression coefficient * Ua <t) register 96 value is subtracted from the first-stage estimation error signal 
e Ua lt) to produce the second-stage estimation error e 2 ^ a lt). The second-stage estimation error signal 
e^t) is a second, somewhat better approximation to either the primary signal s AB it) or the secondary 

20 signal n^U). 

The same processes are repeated in the least-squares lattice predictor 70 and the regression filter 
80 for each stage until a good approximation e mM [l). to either the primary signal s^lt) or the secondary 
signal n^W is determined. Each of the signals discussed above, including the forward prediction error 
yt>, the backward prediction error b m (t), the estimation error signal e^.W, is necessary to calculate the 

25 forward reflection coefficient T^t). the backward reflection coefficient r b Jl). and the regression 
coefficient k^M register 90. 92. and 96 values in each stage, m. In addition to the forward prediction 
error f ra (t). the backward prediction error b ra (t). and the estimation error e^W signals, a number of 
mtermetiate variables, not shown in Figure 8 but based on the values labeled in Figure 8. are required to 
calculate the forward reflection coefficient r f m W. the backward reflection coefficient T^t). and the 

30 regression coefficient ft^N register 90,92, end 96 vahies. 

Intermediate variables include a weighted sum of the forward prediction error squares 3 m M. a 
weighted sum of the backward prediction error squares BJ t), a scalar parameter AJt). a conversion factor 
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Km (t). and another scalar parameter p^.ft). The weighted sum of the forward prediction errors a m <t) 
is defined as: 



t 

9JM - E ^ \\J>\h m 



where A without a wavelength identifier, a or b, is a constant multiplicative value unrelated to wavelength 
and is typically less than or equal to one. U, A * 1. The weighted sum of the backward prediction 
errors B m (t) is defined as: 



10 



f 



a" Ml' 



(45) 



where, again. A without a wavelength 'identifier, a or b. is a constant multiplicative value unrelated to 
wavelength and is typically less than or equal to one. le.. A Z 1. These weighted sum intermediate error 
signals can be manipulated such'that they are more easily solved for. as described in Chapter 9. § 9.3 of 

15 the Haykin book referenced above and defined hereinafter m equations (59) and (60). 

nrscmPTiON riF the joimt p ur cess estimator 
The operation of the joint process estimator 60 is ls follows. When the joint pocess estimator 
60 is turned on. the initial values of intermediate variants «, d signals including the parameter A^lt). the 
weigkd sum of the forward prediction error signals ©^IH the weighted sum of the backward prediction 

20 error signals S^tt). the parameter p^.W. and the zero-stage estimation error ^.(t) are outiafiied. 
some to zero and some to a small positive number 6: 



(46) 



25 



(47) 
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WO)-* 

After initiaHzatior, a simultaneous sample of the measured signal S^W or S, b (t) and either the 
secondary reference n'(t> o, the primary reference fld are input to the joint process estimator 60. « 

variables including the weighted sums of the forward and backward error s,gnals S 0 (t) and 6,(4, 
conversion factor y Q W are calculated for the zero-stage according to: 

f 0 (t) - b 0 W - n'W t51a) 
,. (tV ,2 (52a) 



15 



25 



30 



(53a) 



if a secondary reference n'(t| is used or according to: 

f 0 (t, - b 0 ,t, - s',t, (51W 



V*-i»-V 



(52b) 
(53b) 



if a prunary reference s<(t) is used where, again. . without a wavelength identifier, a or n. is a constant 
multinlicative value unrelated to wavelength. 

rmmi c M .f«n« rya, — - WB " 

rt , te pjT^ Th. taw*- ».«*»« r H * *a*w«d ^ 

x ^ -—■■«,> » « - - — * - - *" 
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intermediate values and register values including the parameter A^ft); the forward reflection coefficient 
r f m (t) register 90 value; the backward reflection coefficient r b(in (t) register 92 value; the forward and 
backward error signals f m (t) and b m (t); the weighted sum of squared forward prediction errors $ 1jm [t), as 
manipulated in § 9.3 of the Haykin book; the weighted sum of squared backward prediction errors \JQ. 
as manioc ated in S 9.3 of the Haykin book; the conversion factor yjl): the parameter P^tt); the 
regression coefficient /r^.M register 96 value; and the estimation error e,,,, Wi W value are set according 
to: 



W 1 ' " A ^nvl 1 
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fh .IMIf* «ItlJ|/ i(t-1)l 


(54) 


(t) - -(A iltl/B^ilM)} 


(55) 


1\ (t| - •{A* m . 1 (t)/3F m .|(t)} 


(56) 


f (tl - f ilt) + r*« JtJb-ilM) 
nv 1 m-1 1 ' t # nr * nM 


(57) 




(58) 


OJtt-O^lW-d.^lWl 2 ^!^) 


(59) 




* 

(60) 


K m ft-D - Kn-ift-U • {l"»i(M)| 2 »»ift-1)) 


(61) 


PnvlaW " ^Maft- 1 ' * (^,l^"aUa«'nJ«) 


(62) 




(63) 




(64) 



where a (*) denotes a complex conjugate. 

These equations cause the error signals f m <t). bjt). e^.W to be squared or to.be multiplied by 
another, in effect squaring the errors, and creating new intermediate error values, such as A^W. 
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The error signals and the intermediate error values are recursively tied together, as shown in the above 
equations (54) through (64). They interact to minimfce the error signals in the next stage. 

After a good approximation to either the primary signal s,, a (t) or the secondery signal n A J\) has 
been determined by the joint process estimator 60. a next set of samples, including a sample of the 

5 measured signal S^lt) and a sample of either the secondary reference n'(t) or the primary reference s'W, 
are input to the joint process estimator 60. The reinitialization process does not re-occur, such that the 
forward and backward reflection coefficient r f>m (t) and T M (t) register 90. 92 values and the regression 
coefficient jc^t) register 96 value reflect the multiplicative values required to estimete either the primary 
signal portion sj.W or the secondary signal portion n^M of the sample of S At fX) input previously. Thus. 

10 information from previous samples is used to estimate either the primary or secondary signal portion of a 

present set of samples in each stage. 

In a more numerically stable and preferred embodiment of the above described joint process 
estimator, a normalized joint process estimator is used. This version of the joint process estimator 
normalizes several variables of the above-described joint process estimator such that the normalized 
15 variables fan between -1 and 1. The derivation of the normalized joint process estimator is motivated in 
the Haykin text as problem 12 on page 640 by redefining the variables defined according to the following 
conditions: 
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TOs tionshrmMioi. allow, too comeraw of Equation. |54) • (64! » <"» I*"** 



equations: 
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Initialbation of ftfnrmafaed Joint Process Estimator 

Ut N(t) be defined as the reference noise input at time index n and U(t) be defined as combined 
signal plus noise eiput at time index t the following equations apply |m» Haykin, p. 619): 

1. To initiaBze the algorithm, at time t-0 set 

A».,<P)-0 
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Yo(0)=1 

2. At each instant tfcl. generate the various zeroth-order variables as follow* 

Yo('-D= 1 



shown 



3. Fot rBgiossioii Stoma, inWaBa tin algorithm l>y selling « time Max t-0 

P m (0)=0 

4. ' A^eaer; instant tt1, generate the zeroth-order variable 

e o (/) = ^(0. 



Accordingly, a norniafeed joint process estimator can be used for a more stable system. 
In yet another einbooTment, the correlation cancellation is performed with a QRD algorithm as 
^grammatically in Figure 8a and as described in detail in Chapter 18 of Adaptive Filter Thepry by 



Simon Haykin. published by Prentice-HaH, copyright 1988. 
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The following equations adapted from the Haykin book correspond to the QRD-LSL diagram of 
Figure 8a (also adapted from the Haykin book). 
Computations 

a. Predictions: For time t - 1, 2, and prediction order m - 1, 2, .... M. where M is the 
final prediction order, compute: 
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* 1/2 P 1 *l('- 2 > 



S fc./77-l(' „i/2 



p£i(M) 
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b. Filtering : For order m - 0, 1, _ , M-1; and time t-1. 2, compute 

|2 



1/2/ 



e„.itt-*3SiC?««.itt 
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5, Initialization 

a. Auxiliary parameter initialization: for order m - 1, 2, ~.» M# set 

t/.m-l(0)=" fcm -l(0)=0 
Pm(0)= 0 

b # Sof| constraint initialization : For order m - 0, 1, ...» M, set 

P m (-D=6 
•^m(0) = 8 



where £ is a small positive constant. 

c nata initialization : For t - 1, 2, compute 



e o (0=<M 



YoW-1 



10 where j/(t) is the input and d(t) is the desired response at time L 

PinWCHART OF -iniWT PROCES S ESTIMATOR 
In a signal processor, such as a physiological monitor incorporating a reference processor of the 
present invention to determine a reference n'(t) or s'(t) for input to a correlation canceler. a joint process 
estimator 60 type adaptive noise canceler is generally implemented via a software program having an 
15 iterative bop. One iteration of the loop is analogous to a single stage of the joint process estimator as 
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shown in Figure 8. Thus, if a loop is iterated m times, it is equivalent to an m stage joint process 
estimator 60. 

A flow chart of a subroutine to estimate the primary signal portion s At i\) or the secondary signal 
portion n^ a (t) of a measured composite signal. S Aa it) is shown in Figure 9. The flow chart illustrates the 
5 function of a reference processor for determining either the secondary reference n'(t) or the primary 
reference s'lt). The flowchart for the joint process estimator is implemented in software. 

A one time initialization is performed when the physiological monitor is powered-on, as indicated 
by an •INITIALIZE NOISE CANCELER" action block 12a The initialization sets all registers 90. 92, and 96 
and delay element variables 110 to the values described above in equations (46) through (50). 

10 Next, a set of simultaneous samples of the composite measured signals S^ a (t) and S^ b (t) is input 

to the subroutine represented by the flowchert in Figure 9. Then a time update of each of the deley 
element program variables occurs, as indicated in a "TIME UPDATE OF rZ 1 ] ELEMENTS" action block 130. 
The value stored in each of the delay element variables 110 is set to the value at the input of the delay 
dement variable 1 10. Thus, the zero-stage beckward prediction error DqW is stored as the first-stage delay 

15 element variable, and the first-stage backward prediction error b,{t} is stored as the second-siege delay 
element variable, and so on. 

Then, using the set of measured signal samples S^t) end S^t), the reference signal is obtained 
using the ratiometric or the constant saturation methods described above. This is indicated by a 
-CALCULATE REFERENCE [n'(t) or s'(t)] FOR TWO MEASURED SIGNAL SAMPLES" action block 140. 

20 A zero stege order update is performed next as indicated in a "ZERO-STAGE UPDATE" action block 

150. The zero-stage backward prediction error b 0 (t), and the zero-stage forward prediction error f„(t) are 
set equal to the value of the reference signal n'(t) or At). Adtfrtionelly, the weighted sum of the forward 
prediction errors 9 m (t) and the weighted sum of backv ard prediction errors BJU are set equal to the value 
defined in equations (47) and (481 

25 Next, a loop counter, m. is initialized as imficated in a "m-0" action block 160. A maximum value 

of m, defining the total number of stages to be used by the subroutine corresponding to the flowchart in 
Figure 9, is also defined. Typically, the loop is constructed such that it stops iterating once a criterion for 
convergence upon a best approximation to either the primary signal or the secondary signal has been met 
by the jooit process estimator 60. Additionally, a maximum number of loop iterations may be chosen at 

30 vhich the loop stops iteration, bi a preferred embodiment of a physiological monitor of the present 
: vention. a maximum number of iterations, m-6 to m-10. is advantageously chosen. 

Within the loop, the forward and backward reflection coefficient T^lt) and T^t) register 90 
and 92 values in the least-squares lattice filter ere calculated first es indiceted by the "ORDER I PDATE 
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MTH CELL OF LSL-LATTICE" action block 170 in Figure 9. This requires calculation of mterrnefiate 
variable and signal values used in determining register 90. 92. and 96 values in the present stage, the next 

stage, and In the regression filter 80. 

The calculation of regression filter register 96 value ic^.W is performed next indicated by the 

5 "ORDER UPDATE MTH STAGE OF REGRESSION FILTER(S)" action block 180. The two order update action 
blocks 170 and 180 are performed in sequence m times, until m has reached its predetermined nuuimum 
fin the preferred embodiment, m-6 to .-10) or a solution has been converged upon, as indicated by . 
YES path from a "DONE" decision block 190. In a computer subroutine, convergence is deterrruned by 
checking if the weighted sums of the forward and backward prediction errors QJ* and B m W are less than 

10 a small positive number. An output is calculated next, es indicated by a "CALCULATE OUTPUT" action 
block 200. The output is a good approximation to either the primary signal or secondary s.gnal, as 
determined by the reference processor 26 and jobt process estimator 60 subroutine corresponding to the 
flow chart of Figure 9. This is displayed (or used in a calculation in another subroutine), as md.c.ted by 

a "TO DISPLAY" action block 210. 

15 A new set of samples of the two measured signals S/la (t) and S, D (t) is input to the processor 

and joint process estimator 60 adaptive noise canceler subroutine corresponding t. the flowchart of 
Figure 9 and the process reiterates for these samples. Note, however, that the initialization process does 
not re-occur. New sets of measured signal samp.es S>la (t) and S Ab M are continuously input to the 
reference processor 26 and joint process estimator adaptive noise cenceler subroutine. The output forms 

20 a chain of samples which is representative of a continuous wave. This waveform is a good approximator, 
to either the primary signal waveform s, a (U or the secondary waveform n, a (t> at wavelength Az. The 
waveform may also be e good approximation to either the primary signal wavefonn s, B (t) or the secondary 

waveform fC Ah [\) at wavelength A*. ^ 

A corresponding flowchart for the ORD algorithm of Fioure 8. is depicted m Figure 9a. wrth 

25 reference numeral corresponding in number with an V extension 

raj nil ATION OF saturation FP "M rnRRELATlON CANCflFR OUTPUT 
Physiological monitors may use the approximation of the primary signals s^U) or s'^W or the 
secondary signals n-, 8 (t, or n-, h <t) to calculate another quantity, such as the saturation of one constituent 
in a volume containing that constituent plus one or more other constituents. GeneraOy. such calcu.at.ons 
30 require ^formation about either a primary or secondary signal a, two wavelengths. For example, the 
constant saturation method requires a good approximation of the primary signal portions s^t) and s^t) 
of both measured signals S Aa iX) and S, D (t|. The arterial saturetion is determined from the epproxemmon. 
to both signals. Le. s^t) and s-, D (t). The constant saturation method also requires a good approxunat^n 
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of the secondary signal portions n^t) or n^ttf. An estimate of the venous saturation may be determined 
from the approximations to these signals L e. n"^ a (t) and rf^ft). 

A joint process estimator 60 having two regression filters 80a and 80b is shown in Figure 10. 
A first regression filter 80a accepts a measured signal S^t). A second regression filter 80b accepts a 
5 measured signal S^ft) for a use of the constant saturation method to determine the reference signal n'W 
or s'(t). The first end second regression filters 80a and 80b are independent. The backward prediction 
error bjt) is input to each regression filter 80a and 80b. the input for the second regression filter 80b 
bypassing the first regression fflter 80a. 

The second regression filter 80b comprises registers 98, and summing elements 108 arranged 
10 similarly to those in the first regression filter 80a. The second regression filter 80b operates via an 
additional intermediate variable in conjunction with those defined by equations (54) through (64), Le.: 

/WbM " + {" m W»*nUbt«»'Kjt)}; and «*) 

15 P<Mb«»- 0 - m 

The second regression filter 80b has an error signal value defined similar to the first regression filter error 
signal values, e^^t), Le.: 

20 ^ M0 (t) - e^it) - /c^fcWb^t); and «*7) 

w« - ,or * * °- 1681 

The second recession filter has a regression coefficient tc^it) register 98 value defined similarly to the 
25 first regression fiter error signal values. Leu 

^nUbW " {*rUbW»mW>: or 169) 

These values are used in coraunction with those intermediate variable values, signal values, register and 
30 register values defined in equations (46) through (64). These signals are calculated in an order defined by 
placing the additional signals immediately adjacent a similar signal for the wavelength An. 

For the constant saturation method, S^t) is input to the second regression tUter 80b. The output 
is then a good approximation to the primary signal sT Ab W or secondary signal s^W. 
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The addition of the second regression filter 80b doe, not substantially change the computer 
program subroutine represented by the flowchart of Figure 9. Instead of . order update of the m stege 
of only one regression filter, an order update of the m» stage of both regression filters 80s and 80b « 
performed. This is characterized by the plural designation in the "ORDER UPDATE OF m* STAGE OF 
REGRESSION FILTERISJ" activity block 180 in figure 9. Since the regression filters 80a and 80b operate 
independently, independent calculations can be performed in the reference processor and ,o,nt process 
estimator 60 adaptive noise canceler subroutine modeled by the flowchart of Figure 9. 

An alternative diagram for the joint process estimator of figure 10. using the QRD algorithm and 
having two regression filters is shown in Figure 10a. This type of joint process estimator would be used 
for correlation cancelation using the QRD elgorithm described in the Haykin book. 

rat mi ATIOM OF SATURATION 

Once good approximations to the primary signal portions s^tt) end s^ b (t) or the secondary signal 
portion. „V«I and n^t), have been determined by the joint process estimator 60. the saturation of A 5 
in a volume containing A 5 and A B . for example, may be celculated according to various known methods. 
Mathematical^ the approximations to the primary signals can be written in terms of and At. as 

= *5^b c 5 x 5.6W ♦ e 6iAb e 6 x 5 B it). I71) 

Equations (70) and (71) are equivalent to two equations having three unknowns, namely c 5 (t). c 6 (t) and 
x RR (t) The saturation can be determined by acquiring approximations to the primary or secondary s.gn.1 
portions at two different, yet proximate times t, and t 2 over which the saturation of A 5 in the volume 
containing A 5 and A 6 and the saturation of A 3 in the volume containing A 3 ami A 4 does not change 
i substantially. For example, for the primary signals estimated at times l, and t 2 : 

hi (72) 



D 



(73) 

s W t l> " e 5>UrVMM * e 6^b c 6 x S.Bt t 1 ) 

(74) 
(75) 



WO 96/12435 PCT/US95/13469 

-49- 

Then, difference signals may be determined which relate the signals of equations (72) through (75). Le.: 



10 



(76) 
(77) 



where Ax - x 5 6 (t,) • x 56 (t 2 ). The average saturation at time t-lVt 2 )J2 is: 

Saturation^) - c 5 (t)flc 5 (t) ♦ c 6 (t)l (78) 



(79) 



15 It win be understood that the Ax term drops out from the saturation calculation because of the division. 
Thus, knowledge of the thickness of the primary constituents is not required to calculate saturation. 

PHLSE OXIMETR Y MEASUREMENTS 
A specific example of a physiological monitor utilizing a processor of the present invention to 
determine a secondary reference n*(t) for input to a correlation canceler that removes erratic motion-induced 
20 secondary signal portions is a pulse oximeter. Pulse oximetry may also be performed utifaing a processor 
of the present invention to determine a primary signal reference s'(t) which may be used for display 
purposes or for input to a correlation canceler to derive information about patient movement and venous 

blood oxygen saturation. 

A pube oximeter typically causes energy to propagate through a medium where blood flows close 
25 to the surface for example, an ear lobe, or a digit such as a finger, a forehead or a fetusT scalp. An 
attenuated signal is measured after propagation through or reflected from the medium. The pulse oximeter 
estimates the saturation of oxygenated blood. 

Freshly oxygenated blood is pumped at high pressure from the heart into the arteries for use by 
the body. The volume of blood in the arteries varies with the heartbeat giving rise to a variation in 
30 absorption of energy at the rate of the heartbeat, or the pulse. 

Oxygen depleted, or deoxygenated. blood is returned to the heart by the veins along with unused 
oxygenated blood. The volume of blood in the veins varies With the rate of breathing, which is typicaly 
much slower than the heartbeat Thus, when there is no motion induced variation in the thickness of the 
veins, venous blood causes a low frequency variation in absorption of energy. When there is motion 
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induced variation in the thickness of the veins, the low frequency variation in absorption is coupled with 
the erratic variation in absorption due to motion artifact 

In absorption measurements using the transmission of energy through a medium, two light emitting 
diodes (LED's) are positioned on one side of a portion of the body where blood flows close to the surface. 

5 such as a finger, and a photodetector is positioned on the opposite side of the finger. Typically, in pulse 
oximetry measurements, one LED emits a visible wavelength, preferably red, and the other LED emits an 
infrared wavelength. However, one skilled in the art will realize that other wavelength combinations could 
be used. The finger comprises skin, tissue, muscle, both artepal blood and venous blooiL fat, etc 
each of which absorbs light energy differently due to different absorption coefficients, different 

10 concentrations, different thicknesses, and changing optical pathlengths. When the patient is not moving, 
absorption is substantially constant except for the flow of blood. The constant attenuetion can be 
determined and subtracted from the signal via traditional filtering techniques. When the patient moves, this 
causes perturbation such as changing optical pathlength due to movement of background fluids (a.*, venous 
blood having a different saturation than the arterial blood). Therefore, the measured signal becomes erratic. 

15 Erratic motion induced noise typically cannot be predetermined andfor subtracted from the measured signal 
via traditional filtering techniques. Thus, determining the oxygen saturation of arterial blood and venous 

blood becomes more difficult . 

A schematic of a physiological monitor for pulse oximetry is shown in Figures 11-13. figure 11 
depicts » general hardware block diagram of a pulse oximeter 299. A sensor 300 has two light emitters 
20 301 and 302 such as LED's. One LED 301 emitting fight of red wavelengths end another LED 302 em.tt.ng 
Bght of infrared wavelengths are placed adjacent a finger 310. A photodetector 320. which produces an 
electrical signal corresponding to the attenuated visible and infrared fight energy sign* « located opposite 
the LED's 301 and 302. The photodetector 320 is connected to front end analog signal conditioning 
circuity 330. 

25 The front end analog signal conditioning circuitry 330 has outputs coupled to analog to digital 

conversion circuit 332. The analog to digital conversion, circuitry 332 has outputs coupled to a digital 
signal processing system 334. The digital signal processing system 334 provides the desired parameters 
as outputs for a display 338. Outputs for display are. for example, blood oxygen saturation, heart rate, 
and a clean pfethysmographic waveform. 

30 The signal processing system also provides an emitter current control output 337 to a digital-to- 

analog converter circuit 338 which provides control mforrnation for light emitter driver* 340. The fight 
emitter drivers 340 couple to the fight emitters 301. 302. The digital signal processing system 334 also 
provides a gain control output 342 for the front end analog signal conditioning circuitry 330. 
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Figure 11a illustrates a preferred embodiment for the combination of the emitter drivers 340 and 
the digital to analog conversion circuit 338. As depicted in figure 11a, the driver comprises first and 
second input latches 321, 322. a synchronbing latch 323. a voltage reference 324. a digital to analog 
conversion circuit 325. first and second switch banks 326, 327, first and second voltage to current 
5 converters 328. 329 and the LED emitters 301, 302 corresponding to the LEO emitters 301. 302 of Figure 
11. 

The preferred driver depicted in Figure 1 la is advantageous in that the present inventors recognized 
that much of the noise in the oximeter 299 of Figure 11 is caused by the LED emitters 301. 302. 
Therefore, the emitter driver circuit of Figure 11a is designed to rranimiie the noise from the emitters 301, 

10 302. The first and second input latches 321. 324 are connected directly to the DSP bus. Therefore, these 
latches significantly mmirnizes the bandwidth (resulting in noise) present on the DSP bus which pssses 
through to the driver circuitry of Figure 11a. The output of the first and second input latches only changes 
when these latched detect their address on the DSP bus. The first input latch receives the setting for the 
digital to analog converter circuit 325. The second input latch receives switching control data for the 

15 switch banks 326. 327. The synchronizing latch accepts the synchronizing pulses which maintain 
synchronization between the activation of emitters 301. 302 end the analog to digital conversion circuit 
332. 

The voltege reference is also chosen as a low noise DC voltage reference for the digital to analog 
conversion circuit 325. In addition, in the present embodiment the voltage reference has an lowpass 

20 output filter with a very low comer frequency (e.g,, 1Hz in the present embodiment). The digital to analog 
converter 325 also has a lowpass filter at its output with a very low corner frequency |e.g., 1 Hz). The 
digital to analog converter provides signals for each of the emitters 301. 302. 

In the present ernbodinent the output of the voltage to current converters 328, 329 ere switched 
such that with the emitters 301, 302 connected in back-to-back configuration, only one emitter is active 

25 an any given time. In addition, the voltage to current converter for the inactive emitter is switched off 
at its input as wen, such that it is completely deactivated. This reduces noise from the swHching and 
voltage to current conversion circuitry. In the present embodiment, low noise voltage to current converters 
are selected (e.g. Op 27 Op Amps), end the edback loop is configured to have e low pass filter to reduce 
noise. In the present embodiment, the low Pj ss fE.ering function of the voltage to current converter 328. 

30 329 has a corner frequency of just above 625 to. . ;hicb is the switching speed for the emitters, as further 
discussed below. Accordingly, the preferred flhrer circuit of Figure 11a, mmirmzes the noise of the emitters 
301, 302. 
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In general, the red and infrared light emitters 301. 302 each emits energy which is absorbed by 
the finger 310 and received by the photodetector 320. The photodetector 320 produces an electrical signal 
which corresponds to the intensity of the light energy striking the photSa tector 320. The front end analog 
signal conditioning circuitry 330 receives the intensity signals and filters and conditions these signals as 

5 further described below for further processing. The resultant signals are provided to the analog-todigital 
conversion circuitry 332 which converts the analog signals to digital signals for further processing by the 
digital signal processing system 334. The digital signal processing system 334 utilizes the two signals fa, 
order to provide a what will be called herein a "saturation transform," It should be understood, that for 
parameters other than blood saturation monitoring, the saturation transform could be better termed as a 

10 concentration transform, in-vivo transform, or the Eke, depending on the desired parameter. The term 
saturation transform is used to describe an operation which converts the sample data from time domain 
to saturation domain values as will be apparent from the discussion below. In the present embodiment, 
the output of the digital signal processing system 334 provides clean plethysmography waveforms of the 
detected signals and provides values for oxygen saturation and pulse rate to tha display 336. 

15 It should be understood that in different embodiments of the present invention, one or more of the 

outputs may be provided. The digital signal processing system 334 also provides control for driving the 
light emitters 301. 302 with an emitter current control signal on the emitter current control output 337. 
This value is a digital valua which is converted by the digital to-analog conversion circuit 338 which 
provides a control signal to the emitter current drivers 340. The emitter current drivers 340 provide the 

20 appropriate current drive for the red emitter 301 and the infrared emitter 302. Further detail of the 
operation of the physiological monitor for pulse oximetry is explained below. In the present 

embodiment, the light emitters are driven via the emitter current driver 340 to provide light transmission 
with digital modulation at 625 Hi. In the present embodiment the fight emitters 301. 302 are driven at 
a power level which provides an acceptable intensity for detection by the detector and for convening by 

25 the front end analog signal conditioning circuitry 330. Once this energy level is determined for a given 
patient by the digital signal processing system 334. the current level for the red and infrared emitters is 
maintained constant. It should be understood, however, that tha current could be adjusted for changes in 
the ambient room Fight and other changes which would effect tha voltage input to the front end analog 
signal conditioning circuitry 330. In the present invention, the red and infrared light emitters are modulated 

30 as follows: for one complete 625 Hz red cyde. tha red emitter 301 is activated for the first quarter cycle, 
and off for the remaining three-quarters cycle; for one complete 625 Hi infrared cycle, the infrared light 
emitter 302 is activated for one quarter cycle and is off for the remaining three-quarters cyde. In order 
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to only receive one signal at a time, the emitters ere cycled on and off alternatively, in sequence, with each 
only active for a quarter cycle per B25 Hi cycle and a quarter cycle separating the active times. 

The light signal is attenuated (amplitude modulated) by the pumping of blood through the finger 310 
(or other sample medium). The attenuated (amplitude modulated) sigru: is detected by the photodetector 
5 320 at the 625 Hz cerrier frequency for the red and infrared light. B. cause only a single photodetector 
is used, the photodetector 320 receives both tha red and infrared signals to form a omposite time division 
signal 

The composite time division signal is provided to the front analog signal conditioning circuitry 330. 
Additional detail regarding the front end analog signal conditioning circuitry 330 and the analog to digital 

10 converter circuit 332 is ftlustreted in Figure 12. As depicted in Figure 12. the front end circuity 302 has 
a preamplifier 342. a high pass filter 344, an amplifier 346. a progremmable gain amplifier 348, and a tow 
pass filter 350. The preamplifier 342 is a transimpedance amplifier that converts the composite current 
signal from the photodetector 320 to a corresponding voltage signal, and amplifies the signal. In the 
present embodiment, the preamplifier has a predetermined gain to boost the signal amplitude for ease of 

15 processing. In the present embodiment, the source voltages for the preamplifier 342 are -15 VOC and ♦ IS 
VDC. As will be understood, the ettenuated signal contains a component representing ambient Eght as wed 
as the component representing the infrared or the red light as the case may be in time. If there is fight 
in the vicinity of the sensor 300 other than the red end infrared light, this ambient light is detected by the 
photodetector 320. Accordingly, the gain of the preamplifier is selected in order to prevent the ambient 

20 Ignt in the signal from saturating the preampfifier under normal and reasonable operating conditions. 

In the present embodiment, the preamplifier 342 comprises en Analog Devices AD743JR OpAmp. 
This transimpedance amplifier is particularly advantageous in that it exhibits several desired features for 
the system described, such as: low equivalent input voltage noise, low equivalent input current noise, low 
input bias current, high gain bandwidth product, low total harmonic distortion, high common mode rejection, 

25 high open loop gain, and a high power supply rejection ratio. 

The output of the preamplifier 342 couples as an input to the high pass filter 344. The output 
of the preempOfier also provides e first input 346 to the analog to digital conversion circuit 332. In the 
present embodiment, the high pass filter is a single-pole filer with a comer frequency of about 1/2 • 1 Hz. 
However, the corner frequency is readily raised to about 90 Hz in one embodiment. As wi be understood, 

30 the 625 Hz carrier frequency of the red and infrared signals a weB above a 90 Hz corner frequency. The 
high-pass filter 344 has an output coupled as an input to en ampifier 346. In the present enibodknent. 
the arnpfifier 346 comprises a unity gain amplifier. However, the gain of the amplifier 346 is adjustable 
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by the variation of e single resistor. The gam of the amplifier 34B would be increased if the gain of the 
preamplifier 342 is decreased to compensate for the effects of ambient fight 

The output of the amplifier 346 provides an input to a programmable gain amplifier 348. The 
programmable gain amplifier 348 also accepts a programming input from the digital signal processing system 

5 334 on e gain control signal line 343. The gain of the programmable gain amplifier 348 is digitally 
programmable. The grin is odjusted dynamically at Wtiafaation or sensor placement for changes in the 
medium under test from patient to patient. For example, the signal from different fingers differs somewhat. 
Therefore, a dynamically adjusteble amplifier is provided by the programmable gain amplifier 348 in order 
to obtain a signal suitable for processing. 

10 The programmable gain ampDfier is also advantageous in en ehernative embodiment in which the 

emitter drive current is held constant. In the present embodiment the emitter drive current is edjusted for 
each patient in order to obtain the proper dynamic renge at the input of the analog to digital conversion 
circuit 332. However, changing the emitter drive current can alter the emitter wavelength, which in turn 
affects the end result in oximetry calculetions. Accordingly, it would be advantageous to fix the emitter 

15 drive current for en patients. In en alternative embodiment of the present invention, the programmable gain 
empnfier can be edjusted by the OSP in order to obtein e signal at the input to the analog to digital 
conversion circuit which is properly within the dynamic range {♦ 3v to -3v in the present embodiment) of 
the analog to digital conversion circuit 332. In this manner, the emitter drive current could be fixed for 
an patients, eOmineting the wevelength shift due to emitter current drive chenges. 

20 The output of the programmable gain amplifier 348 couples as an input to a low-pass filter 350. 

Adventegeously, the low pass filter 350 is a single-pole filter with e corner frequency of eporoximately 10 
Khz in the present embodiment. This low pass filter provides anti-efiasing in the present embodiment 

The output of the low-pass filter 350 provides a second input 352 to the analog-to-digitel 
conversion circuit 332. Figure 12 also depicts additional defect of the enalogto-digrtal conversion circuit 

25 In the present embodiment the anelog-to-digital conversion circuit 332 comprises a first analog-to-digital 
converter 354 end a second enelog-to-digital converter .35B. Advantageously, the first analog-to-digital 
converter 354 eccepts input from the first input 346 to the andog-to-digHel conversion circuit 332, and 
the second analog to digital converter 356 accepts input on the second input 352 to the analog-to-digital 
conversion circuitry 332. 

30 In one advantageous embodiment the first endofrto-digital converter 354 is a diegnostic analog-to- 

digital converter. The diagnostic task Iperformed by the digital signal processing system) is to read the 
output of the detector es ampfified by the preemplifier 342 in order to determine if the signal is seturetimj 
the input to the high-pass filter 344. In the present enfcodiment if the input to the high pass filter 344 
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becomes saturated, the front end analog signal conditionmB circuits 330 provides a '0* output. 
Alternatively, the first analog-to-digital converter 354 remains unused. 

The second analog-to-digital converter 352 accepts the conditioned composite analog signal from 
the front end signal conditioning circuitry 330 and converts the signal to digital form. In the present 

5 embodiment the second analog to digital converter 356 comprises a single-channel delte-sigma converter. 
In the present embodiment, a Crystal Semiconductor CS5317KS delta-sigma analog to digital converter is 
used. Such a converter is advantageous in that it is low cost, and exhibits lot noise characteristics. More 
specifically, a delta-sigma converter consists of two major portions, a noise modulator and a decimation 
rater. The selected converter uses a second order analog delta-sigma modulator to provide noise shaping. 

10 Noise shaping refers to changing the noise spectrum from a flat response to e response where noise at the 
lower frequencies has been reduced by increasing noise at higher frequencies. The decimation filter then 
cuts out the reshaped, higher frequency noise to provide IB-hit performance at a lower frequency. The 
present converter samples the data 128 times for every 16 bit data word that it produces. In this manner, 
the converter provides excellent noise rejection, dynamic range and low harmonic distortion, that help in 

15 critical measurement situations Gke low perfusion and electrocautery. 

In addition, by using a single-channel converter, there is no need to tune two or more channels to 
each other. The delta-sigmB converter is also advantageous in that it exhibits noise shaping, for improved 
noise control. An exemplary analog to digital converter is a Crystal Semiconductor CS5317. In the pk sent 
embodiment, the second analog to digital converter 356 samples the signal at a 20 Khz sample rate. The 

20 output of the second analog to digital converter 356 provides data samples at 20 Khz to the digital signal 
processing system 334 (Figure 11). 

The digital signal processing system 334 is illustrated in additional detal in Figure 13. In the 
present embodiment, the digital signal processing system comprises s microcontroller 360, e digital signal 
processor 362, a program memory 364. a sample buffer 366, a data memory 368, a read only memory 

25 370 and communication registers 372. hi the present embodiment the digital signal processor 362 is an 
Analog Devices AD 21020. In the present embodiment the microcontroller 360 comprises a Motorola 
68HC05, with built in program memory. In the present embodiment the sample buffer 366 is a buffer 
which accepts the 20 Khz sample data from the analog to digital conversion circuit 332 for storage in the 
data memory 368. In the present embodiment the deta memory 368 comprises 32 KWords (words being 

30 40 bits m the present embodiment) of static random access memory. 

The microcontroller 360 is connected to the DSP 362 via a conventional JTAG Tap foe. The 
microcontroller 360 ttarisnuts the boot loader for the DSP 362 to the program memory 364 via the Tap 
line, and then allows the DSP 362 to boot from the program memory 364. The boot loader in program 
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memory 364 then causes the transfer of the operating instructions for the DSP 362 from the read only 
memory 370 to the program memory 364. Advantageously, the program memory 364 is a very high speed 

memory for the DSP 362. 

The microcontroller 360 provides the emitter current control and gain control signals via the 

communications register 372. 

Figures 14-20 depict functional block diagrams of the operations of the pulse oximeter 299 carried 
out by the digital signal processing system 334. The signal processing functions described below are 
carried out by the DSP 362 in the present embodiment with the microcontroller 360 pmMnp system 
management. In the present embodiment, the operation is software/firmware controlled. Figure 14 depicts 
a generalized functional block diagram for the operations performed on the 20 Khz sample data enter** 
the digital signal processing system 334. As illustrated in Figure 14, a demodulation, as represented m a 
demodulation module 400, is first performed. Decimation, as represented in a decimation module 402 « 
then performed on the resulting data. Certain statistics are calculated, as represented in a stat.st.cs 
module 404 and a saturation transform is performed, as represented in . saturation transform module 406, 
on the data resulting from the decimation operation. The data subjected to the statistics operations and 
the data subjected to the saturation transform operations ere forwarded to saturation operefons. as 
represented by a saturation calculation module 408 and pulse rate operations, as represented in a pulse 

rate calculation module 410. 

In general, the demodulation operation separates the red and infrared signals from the composite 
signal and removes the 625 Hi carrier frequency, leaving raw data points. The raw data point, are 
provided at 625 Hz intervals to the decimation operation which reduces the samples by en order of 10 to 
samples at 62.5 Hz. The decimation operation also provides some filtering on the samples. The result** 
data is subjected to statistic and to the saturation transform operations in order to calculate a saturat.cn 
value which » very to.er.nt to motion artifacts and other noise in the signal. The saturation value * 
ascertained in the saturation calculation modulo 408. and a pulse rate and a dean plethysmograptac 
waveform is obtained through the pulse rate module 410. Additional detail regarding the variou. operate 
is provided in connection with Figures 15-21. 

Figure 15 aiustrates the operation of the demodulation module 400. The modulated signal format 
i, depicted In Figure 15. One full 625 Hz cycle of the composite signal is depicted in Rgure 15 with the 
first quarter cycle being the active red fight plus ambient light signal the second quarter cycle bemg an 
ambient light signal the third quarter cycle bemg the active infrared plus ambient Gght signal .nd the 
fourth quarter cycl. being an ambient light signal. As depicted in Rgure 15. with . 20 KHz semphng 
frequency/the single full cycle at 625 Hz described above comprises 32 samples of 20 KHz data. e.ght 
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samples relating to red phis ambient light, eight samples relating to ambient light, eight samples relating 
to infrared plus ambient light, and finally eight samples related to ambient fight 

Because the signal processing system 334 controls the activation of the light emitters 300, 302, 
the entire system is synchronous. The data is synchronously divided (and thereby demodulated) into four 
8-sample packets, with a time division demultiplexing operation as represented in a rJemultiplexirig module 
421. One eight-sample packet 422 represents the red phis ambient Dght signal; a second eight-sample 
packet 424 represents an ambient fight signal; a third eight-sample packet 426 represents the attenuated 
infrared light plus ambient light signal; and a fourth eight-sample packet 428 represents the ambient light 
signal A select signal synchronously controls the demurring operation so as to dhrida the time-division 
multiplexed composite signal at the input of the demultiplexer 421 into its four subparts. 

A sum of the last four samples from each packet is then calculated, as represented in the summing 
operations 430. 432. 434. 436 of Figure 15. In the present embodiment, the last four samples are used 
because a low pass filter in the analog to digital converter 356 of the present embodiment has a settDng 
time. Thus, collecting the last four samples from each 8-sampla packet allows the previous signal to clear. 
Tlus summing operation provides an integration operation which enhances noisB nummity. The sum of the 
respective ambient fight samples is then subtracted from the sum of the red and infrared samples, as 
represented m the subtraction modules 438. 440. The subtraction operation provides some attenuation of 
the ambient fight signal present in the data. In the present embodiment, it has been found that 
approximately 20 dB attenuation of the ambient Hght is provided by the operations of the subtraction 
modules 438. 440. The resultant red and infrared sum values are divided by four, as represented in the 
divide by four modules 442, 444. Each resultant value provides one sample each of the red and infrared 
signals at 625 Hz. 

It should be understood that the 625 Hi carrier frequency has been removed by the demodulation 
operation 400. The 625 Hz sample data at the output of the demodulation operation 400 is sample data 
without the carrier frequency. In order to satisfy Nyquist sampling requirements, tess than 20 Hz is needed 
(understanding that the human puisa is about 25 to 250 beats per minute, or about .4 Hz • 4 Hz). 
Accordingly, the 625 Hz resolution is reduced to 62.5 Hz in the decimation operation. 

Figure 16 illustrates the operations of the decimation module 402. The red and infrared sample 
data is provided at 625 Hz to respective red and infrared bufferltHters 450. 452. In the present 
embedment, the red and infrared buffer/filters are 519 samples deep. Advantageously, the buffer filers 
450, 452 function as continuous first-in, first-out buffers. The 519 samples are subjected to low-pass 
ffterwig. Preferably, the low-pass filtering has a cutoff frequency of approxanateh/ 7.5 Hz with attenuation 
of approximately -110 dB. The buffer/filters 450, 452 form a Finite Impulse Response IFIR) filter with 
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coefficients for 519 taps. In order to reduce the sample frequency by ten. the tow-pass filter calculation 
is performed every ten samples, as represented in respective red and infrared decimation by 10 modules 
454, 456. In other words, with the transfer of each new ten samples into the buffer/fflters 450 f 452. 
a new low pass filter calculation is performed by multiplying the impulse response (coefficients) by the 519 
5 filter taps. Each filter calculation provides one output sample for respective red and infrared output buffers 
458. 460. In the present embodiment, the red end infrared output buffers 458, 460 are also continuous 
FIFO buffers that hold 570 samples of data. The 570 samples provide respective infrared and red samples 
or packets (also denoted -snapshot" herein) of samples. As depicted in Figure 14, the output buffers 
provide sample data for the statistics operation module 404. saturation transform module 406, and the 

10 pulse rate module 410. 

Figure 17 illustrates additional functional operation details of the statistics module 404. In 
summary, the statistics module 404 provides first order oximetry calculations and RMS signal values for 
the red and infrared channels. The statistics module also provides a cross-corretetion output which 
indicates a cross-correlation between the red and infrared signals. 
15 As represented in Figure 17. the statistics operation accepts two packets of samples (e.g,, 570 

samples at 62.5 Hz in the present embodiment) representing the attenuated infrared and red signals, with 
the carrier frequency removed. The respective packets for infrared and red signals are normalized with • 
log function, as represented in the Log modules 480, 482. The normalizetion is followed by removal of 
the OC portion of the signals, es represented m the DC Removal modules 484, 486. In the present 
20 embodiment, the DC removal involves ascertaining the DC value of the first one of the samples (or the 
mean of the first several or the mean of an entire snapshot) from each of the respective red and infrared 
snapshots, and removing this DC value from all samples in the respective packets. 

Once the DC signal is removed, the signals are subjected to bendpass filtering, as represented m 
red and infrared Bandpass Filter modules 488, 490. In the present embodiment, with 570 samples in each 
25 packet, the bandpass filters are configured with 301 taps to provide a FIR filter with a linear phase 
response and little or no distortion. In the present embodiment, the bandpass filter has a pass band from 
34 beats/minute to 250 beatsjminute. The 301 taps sT.de over the 570 samples in order to obtein 270 
filtered samples representing the filtered red signal and 270 filtered samples representing the filtered 
infrered signal In an ideal case, the bandpass filters 488, 490 remove the DC in the signal. However, 
30 the DC removal operations 484, 486 assist in DC removal m the present einboaTmeni. 

After filtering, the last 120 samples from each packet (of now 270 samples in the present 
embodimenU «* selected for further processing as represented in Select Last 120 Sample, modules 492, 
494. The last 120 samples are selected because, in the present embodiment, the first 150 samples fad 
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within the settling time for the Satuution Lensfer rnotkle 406, which processes the seme data packets, 

as further discussed below. 

Conventional saturation equation calculations are performed on the red and infrared 120-sampte 
packets. In the present embodiment, the conventional saturation calculations are performed in two different 

5 ways. For one celculation, the 120-sample packets are processed to obtein their overall RMS value, as 
represented in the first red end infrared RMS modules 496. 498. The resultant RMS values for red and 
infrared signals provide input values to a first REO_RMS/IR_RMS ratio operation 500. which provides the 
RMS red value to RMS infrared value ratio as an input to a saturation equation module 502. As well 
understood in the art, the ratio of the intensity of red to infrared ettenuated fight as detected for known 

10 red and infrared wavelengths (typically - 650 nm and A m - 910 nm) relates to the oxygen saturation 
of the patent Accordingly, the saturation equation module 502 represents e conventional look-up table 
or the like which, for predetermined ratios, provides known saturation values at hs output 504. The red 
and infrared RMS values are also provided as outputs of the statistics operations module 404. 

In addition to the conventional saturation operation 502, the 120-sample packets are subjected to 

15 a cross-correlation operation as represented in a first cross-correlation modute 506. The first cross- 
correlation module 506 determines if good correlation exists between the infrered and red signals. This 
cross correlation is advantageous for detecting defective or otherwise malfunctioning detectors. The cross 
correlation is also advantageous in detecting when the signal model (i.e* the model of Equations (1) - (3)) 
is satisfied. If correlation becomes too low between the two channels, the signal model is not met. In 

20 order to determine this, the normalized cross correlation can be computed by the cross-correlation module 
506 for each snapshot of data. One such correlation function is as follows: 



If the cross correction is too low. the oximeter 299 provides s warning {e.g., audMe, visual, etc.) 
to the operator. In the present, embodiment, if a selected snapshot yields a normalized correlation of less 
25 than 0.75. the snapshot does not qualify. Signals which satisfy the signal model wiB have a correlation 

greater than the threshold. 

The red and infrared 120-sample packets are also subjected to a second saturation operation and 
cross correlation in the same manner as described above, except the 120 samples are divided into 5 equal 
bins of samples (Le., 5 bins of 24 samples each). The RMS. ratio, saturation, and cross correlation 
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operations are performed on a bin-by-bin basis. These operations ere represented in the Divide Into five 
Equal Bins modules 510. 512, the second red end infrared RMS modules 514, 516. the second RED- 
RMS/IRRMS ratio module 518. the second saturation equation module 520 and the second cross correlation 
module 522 in Figure 17. 

5 Figure 18 illustrates additional detail regarding the saturation transform module 406 depicted in 

Figure 14. As illustrated in Figure 18, the saturation transform module 406 comprises a reference 
processor 530. e correlation canceler 531. e master power curve module 554. and a bin power curve 
module 533. The saturation transform module 406 can be correlated to figure 7a which has a reference 
processor 26 and a correlation canceler 27 and an integrator 29 to provide a power curve for separate 

10 signal coefficients as depicted in figure 7c The saturation transform module 406 obtains a saturation 
spectrum from the snapshots of data. In other words, the saturation transform 406 provides information 
of the saturation values present in the snapshots. 

As depicted in figure 18. the reference processor 530 for the saturation transform module 406 
has a saturation equation module 532. a reference generator module 534. a DC removal module 536 and 

15 a bandpass filter module 538. The red end infrared 570.sample packets from the decimation operation are 
provided to the reference processor 530. In addition, a plurality of possible saturation values (the 
"saturation axis scan") are provided as input to the saturation reference processor 530. In the present 
embodiment. 117 saturation values era F ovided as the saturation axis scan. In a preferred embodiment, 
the 117 saturation values range uniformly from a blood oxygen saturation of 34.8 to 105.0. Accordingly. 

20 in the present embodiment, the 1 17 saturation values provide an axis scan for the reference processor 530 
which generates a reference signal for use by the correlation canceler 531. In other words, the reference 
processor is provided with each of the saturation values, end a resultant reference signal is generated 
corresponding to the seturation value. The correlation canceler is formed by a joint process estimetor 550 
and a tow pass filter 552 in the present embodiment 

25 It should be understood that the scan values could be chosen to provide higher or lower resolution 

than 117 scan values. The seen values could also be non-uniformiy spaced. 

As fflustreted in figure 18, the saturation equation module 532 accepts the seturation axis scan 

values as an input and provides e ratio V " - ,B compaii$0n t0 ° f 

figure 7a-7c, this ratio "r n " corresponds to the plurality of scan value discussed above in general. The 
30 saturation equation simply provides a known ratio V Mftnfrared) corresponding t. the seturation value 
received as an input 

The ratio "r n " " P f0vided 88 an input t0 *■ re,er8nCB 534 ' " "* ^ 

sample packets. Tha" reference generator 534 multiplies either the red or infrared samples by the ratio V 
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«nd subtracts the value from the infrared or red samples, respectively. For instance, in the present 
embodiment the reference generator 534 multiplies the red samples by the ratio V «* trtCts *» 
value from the infrared samples. The resulting values become the output of the reference generator 534. 
This operation is completed for each of the saturation scan values H* 117 possible values in the present 
5 embodiment). Accordingly, the resultant data can be described as 1 17 reference signal vectors of 570 data 
points each, hereinafter referred to as the reference signal vectors. This data can be stored in an array 
or the Eke. 

In other words, assuming that the red and infrared sample packets represent the red S^X) and 
infrared S |R lt) measured signeb which have primary sit) and secondary nit) signal portions, the output of 
10 the reference generator becomes the seconder* reference signal n'(t), which complies with the signal model 
defined above, as follows: 

In the present enrolment, the reference signal vectors and the infrared signal are provided as 
input to the DC removal module 536 of the reference processor 530. The DC removal module 536. like 
the DC removal modules 484. 486 in the statistics module 404. ascertains the DC value of the first of the 
15 samples for the respective inputs (or mean of the first several or all samples in a packet) and subtracts 
the respective DC baseline from the sample values. The resulting sample values are subjected to a 
bandpass filter 538. 

The bandpass filter 538 of the reference processor 530 performs the same type of filtering as the 
bandpass filters 488. 490 of the statistics module 404. Accordingly, each set of 570 samples subjected 
20 to bendpass filtering results in 270 remaining samples. The resulting deta at a first output 542 of the 
bandpass filter 538 is one vector of 270 samples (representing the filtered infrared signal in the present 
enibodknent). The resulting data at a second output 540 of the bandpass filter 538. therefore, is 117 
reference signal vectors of 270 data points each, corresponding to each of the saturation axis scan values 
provided to the saturation reference processor 530. 

ft should be understood that the red and infrared sample packets may be switched in their use in 
the reference processor 530.. In addition, it should be understood that the DC removal module 536 and 
the bandpass flier module 538 can be executed prior to input of the data to the reference processor 530 
huaus. the calculations performed in the reference processor are linear. This results in a significant 
processing economy. 

30 The outputs of the reference processor 530 provide first and second inputs to a joint process 

estimator 550 of the type described above with reference to Figure 8. The first input to the joint process 
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estimator 550 is the 270-sampie packet representing the infrared signal in the present embodiment This 
signal contains primary and secondary signal portions. The second input to the joint process estimator » 
the 117 reference signal vectors of 270 samples each. 

The joint process estimator also receives a lambda input 543. a minimum error input 544 and a 
5 number of cells configuration input 545. These parameters are well understood in the art. The lambda 
parameter is often called the -forgetting parameter" for a joint process estimator. The lambda .nput 543 
provides control for the rate of cancellation for the joint process estimator, m the present embodiment, 
lanbda is set to a low value such as 0.8. Because statistics of the signal are non-stationary, a low value 
improves tracking. The minimum error input 544 provides an initialization parameter (conventionally known 
10 as the 'inftiefoation value") for the joint process estimator 550. In the present embodiment the m.mmum 
error value is 10' 8 . This initialization parameter prevents the joint process estimator 500 from d.v.dmg by 
zero upon initial calculations. The number of cells input 545 to the joint process estimator 550 configures 
the number of cells for the joint process estimator. In the present embodiment the number of cells for 
the saturation transform operation 406 is six. As weO understood in the ert for each sine wave, the ,omt 
15 process estimator requires two cells. If there are two sine waves m the 35 - 250 beats/minute range. s« 
. cells ellows for the two heart beat sine waves and ona noise sine wave. 

The joint process estimator 550 subjects the first input vector on the first input 542 to a 
correlation cancellation based upon each of the plurality of reference signal vector, provided in the second 
input 540 to the correlation canceler 531 (an 117 reference vectors in sequence in the present 
20 embodiment). The correlation cancellation result, in a single output vector for each of the 117 reference 
vectors. Each output vector represents the information that the first input vector and the corresponding 
reference signal vector do not have in common. The resulting output vectors ere provided as an output 
to the joint process estimator, and subjected to the low pass filter module . 552. In the present 
embodiment, the tow pass filter 552 comprise, a FIR filter with 25 taps and with . corner frequency of 
25 10 Hz with the sampling frequency of 62.5 Hz (Le. at the decimation frequency). 

The joint process estimator 550 of the present embodiment ha, a settling time of 150 data pomts. 
Therefore, the last 120 data point, from each 270 point output vector are used for further processmg. 
In the present erriodiment. the output vector, are further processed together as a whole, end are dnaded 
into . plurality of bins of equal number of data points. As depicted in figure 18. the output vectors are 
30 provided to . master power curve module 554 and to a Divide into five Equal Bins module 556. The Dhnde 
into five Equal Bins module 556 divides each of the output vector, into five bins of equal number of data 
points fr* with 120 data points per vector, each bin has 24 data points). Each bin is then provided to 
the Bin Power Curves module 558. 
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The Master Power Curve module 554 performs e saturation transform as follows: for each output 
vector, the sum of the squares of the data points is ascertained. This provides a sum of squares value 
corresponding to each output vector (each output vector corresponding to one of the saturation scan 
values). These values provide the besis for a master power curve 555, as further represented in Figure 
5 22. The horizontal exis of the power curve represents the saturation axis scan values end the vertical axis 
represents the sum of squares velue (or output energy) for each output vector. In other words, as depicted 
in Figure 22, each of the sum of squares could be plotted with the magnitude of the sum of squares value 
plotted on the vertical "energy output" exis at the point on the horizontal axis of the corresponding 
saturation scan value which generated that output vector. This results in a master power curve 558, an 
10 example of which is depicted in Figure 22. This provides e saturation transform m which the spectral 
content of the attenuated energy is examined by looking at every possible saturation value and examining 
the output value for the assumed saturation value. As will be understood, where the first and second 
inputs to the correlation canceler 531 are mostly correlated, the sum of squares for the correspond^ 
output vector of the correlation canceler 531 wBI be very low. Conversely, where the correlation between 
15 the first and second inputs to the correlation canceler 531 are not significantly corrected, the sum of 
squares of the output vector wl be high. Accordingly, where the spectral content of the reference signal 
end the first input to the correlation canceler are made up mostly of physiological (e*. movement of 
venous blood due to respiration) end non-physiological (e.g., motion induced) noise, the output energy wffl 
be low. Where the spectral content of the reference signal and the first input to the correlation canceler 
20 ere not correlated, the output energy wiB be much higher. 

A corresponding transform is completed by the Bin Power Curves module 558. except a saturation 
transform power curve is generated for each bin. The resulting power curves are provided as the outputs 
of the saturation transform module 406. 

In general in accordance with the signal model of the present invention, there wH be two peaks 
25 in the power curves, as depicted in figure 22. One peak corresponds to the arterial oxygen saturation of 
the blood, and one peak corresponds to the venous oxygen concentration of the blood. With reference to 
the signal model of the present invention, the peek corresponding to the highest saturation value (not 
necessarily the peak with the greatest magnitude) corresponds to the proportionefity coefficient r r In other 
words, the proportionality coefficient r, corresponds to the red/infrared ratio which wiD be measured for 
30 the arterial saturation. SWlarly. peak that corresponds to the lowest saturation value (not necessarily the 
peak with the lowest magnitude) wl generally correspond to the venous oxygen saturation, which 
correS poniis t8 the proportionality coefficient r t in the signel mode) of the present invention. Therefore, 
the proportionality coefficient r v wul be a red/infrared ratio corresponding to the venous oxygen saturation. 
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| n order to obtain arterial oxygen saturation, the peak in the power curves corresponding to the 
highest saturation value could be selected. However, to improve confidence in the value, further processing 
is completed. Figure 19 illustrates the operation of the saturation calculation module 408 based upon the 
output of the saturation transform module 406 and the output of the statistics module 404. As depicted 
in figure 19. the bin power curves and the bin statistics are provided to the saturation calculation module 
408 In the present embodiment, the master power curve, are not provided to the saturation module 408 
but can be displayed for a visual check onsystem operation. The bin statistics contain the red and infrared 
RMS values, the seed saturation value, and a value representing the cross-correlation between the red and 
infrared signals from the statistics module 404. 

The saturation calculation module 408 first determines a plurality of bin ettributes as represented 
by the Compute Bin Attributes module 560. The Compute Bin Attribute, module 560 collects a data bin 
from the information from the bin power curve, and the information from the bin static In the present 
embodiment, this operation involves placing the saturation value of the peak from each power curve 
corresponding to the highest saturation value in the data bin. In the present embodiment, the select™ of 
15 the highest peak i, performed by first computing the first derivative of the power curve in questum by 
convolving the power curve with a smoothing differentiator filter function. In the present embodiment, the 
smoothing differentiator filter function (using a FIR filter) has the following coefficients: 
0.014964670230367 
0.098294046682706 
20 0.204468276324813 
2.717182664241813 
5.704485606695227 
0.000000000000000 
•5.704482606695227 
25 -2.717182664241813 
-0.204468276324813 
-0.098294046682706 
0.014964670230367 

This filter performs the differentiation and smoothing. Next, each point in the original power curve m 
30 question is evaluated and determined to be a possible peak if the following conditions are met: (1) the pomt 
i, at toast 2% of the maximum value in the power curve; (2) the value of the first derivative changes from 
greater than zero to less than or equal to zero. For each point that is found to be a possible peak, the 
neighboring points are examined and the largest of the three points i, considered to be the true peak. 
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The peak width for these selected peaks is also calculated. The peak width of a power curve in 
question is computed by summing, all the points in the power curve and subtracting the product of the 
minimum velue in the power curve and the number of points in the power curve. In the present 
embodiment the peak width calculation is applied to each of the bin power curves. The maximum value 

5 is selected as the peak width. 

In addition, the infrared RMS value from the entire snapshot the red RMS vak». the seed 
saturation value for each bin. and the cross correlation between the red and infrared signals from the 
statistics module 404 ere also placed in the date bin. The attributes are then used to determine whether 
the data bin consists of acceptable date, as represented in a Bin Qualifying Logic module 562. 

10 If the correlation between the red and infrared signals is too low. tha bin is discarded. If the 

saturation value of the selected peek for a given bin is lower than the seed seturation for the seme bin. 
the peek is repleced with the seed saturation value. If either red or infrared RMS value is below a very 
smalt threshold, the bins are all discorded, and no saturation value is provided, because the measured 
signals ere considered to be too small to obtain meaningful data. If no bins contain acceptable data, the 

15 exception handling module 563 provides e menage to the display 336 that the data is erroneous. 

If some bins qualify, those bins that qualify as having accepteble data are selected, and those that 
do not qualify are replaced with the overage of the bins that are eccepted. Each bin is given a time stamp 
in order to meintein the time sequence. A voter operation 565 examines each of the bins and selects the 
three highest seturation values. These values are forwarded to e dip end smooth operation 566. 

20 The cBp and smooth operation 566 basically performs everaging with a low pass filter. The low 

pass filter provides sdjustable smoothing as selected by a Select Smoothing Filter module 568. The Select 
Smoothing Filter module 568 performs its operation based upon a confidence determinetion performed by 
a High Confidence Test module 570. The high confidence test is an examination of the peak width for the 
bin power curves. The width of the peeks provides some indication of motion by the patient - wider peeks 

25 indicating motion. Therefore, if the peaks are wide, the smoothing filter « slowed down. If peaks are 
narrow, the smoothing filter speed is increased. AccordTngly, the smoothing filter 566 is adjusted besed 
on the confidence leveL The output of the clip and smooth module 566 provides the oxygen saturation 
values in accordance with the present invention. 

In the presently preferred embodiment the dip and smooth filter 566 takes each new saturation 

30 value and compares it to the current saturating value. If the magnitude of the difference is less than IB 
(percent oxygen saturation) then the value is pass. Otherwise, if the t, w saturation value is less than the 
filtered saturation velue, the new saturation value is changed to 16 lessen the filtered saturation value. 
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„ tte new saturation value Is greater than the fltered saturation value, then the new saturation value is 
changed to IB more than the filtered saturation van*. 

During high confidence (no motion), the smoothing filter is a simple one-pole or exponent 
smoothing filter which is computed as follows: 

5 y(n) - 0 6 # x|n) * 0 4 * Y(n ' 1, 

where x(n) is the clipped new saturation value, and y(n) is the filtered saturation value. 

During motion condition, a three-pole IIR (infinite impulse response) filter is used. Its characterurtic. 
are controLe- by three ft. constants t,. V an- « c P» ^ - «* °- 900 ' 
The coefficients for a direct form ., IIR filter are computed from these time constant, usmg the foHowmg 

10 relationships: 

a 0 - 0 

e, - t b + (t c Ht a + t b ) 

a 2 - (t^gft, + t b «■ (tjltj) 

a 3 - igViy 



15 



20 



25 



30 



h 0 -1.t t -(t e l(t, + tgftfc»' 
b, - 2(t h )(t c )(t a -1) 

b 2 - (t b )(g(t b * w-itbHgtg-t,) 

Rnure 20 and 21 illustrate the pulse rate module 410 (Figure 14, in greater detail. As Undated 
i„ Figure 20. the heart rate module 410 has a transient remove, and bandpass futer module 578. a n»«on 
artifact suppress module 580. a saturation equation module 582. a motion status module 584. first ami 
second spectral estimation modules 586, 588. a spectrum analysis module 590, . slew rate fining moduta 
592. an output filter 594, end en output filter coefficient module 598. 

As further depicted in Figure 20. the heart rate module 410 accept, the infrared ^an - 570- 
samp* snapshots from tl. output of the delation mod* 402. The heart rate 410 furt^ 

accepts the saturation value which « output from the saturation calculation module 408. .n .ddim the 
max,™ pea, width as calculated by tte confidence test module 570 (same as ^ c**- 
desert above) i, elso provided es en mput to the heart rate module 410. The *M * red s^ 
packets, the saturation value and the output of the motion statu, module 584 are provmed to the mofon 
artifact suppression module 580, 
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The average peak width value provides an input to a motion status module 584. In the present 
embodiment H the peaks are wide, this is taken as an indicaticn of motion. H motion a not detected, 
spectral estimation on the signals is carried out directly without motion artifact suppression. 

In the case of motion, motion artifacts are suppressed using tha motion artifact suppression module 
5 580. The motion artifact suppression module 580 is nearly identical to the saturation transform module 
406. The motion artifact suppression module 580 provides an output which connects as an input to the 
second spectral estimation module 588. tha first and second spectral estimation modules 586. 588 have 
outputs which provide inputs to the spectrum analysis module 580. The spectrum analysis module 580 
also receives an input which is the output of the motion status module 584. The output of the spectrum 
10 analysis module 580 is the initial heart rate determination of the heart rate module 410 and is provided 
as input to the slew rate Smiting modulo 592. The slew rate limiting module 592 connects to the output 
filter 594. The output filter 594 also receives en input from the output filter coefficient module 596. The 
output filter 594 provides the filtered heart rate for the display 336 {Figure 11). 

In the case of no motion, one of the signals (the infrared signal in the present embodiment) is 
15 subjected to DC removal and bandpass filtering as represented in the DC removal and bandpass filter 
module 578. The DC removal and bandpass filter module 578 provide the same filtering es the DC removal 
end bandpass filter modules 536. 538. During no motion conditions, the filtered infrared signal is provided 
to the first spectral estimation module 586. 

In the present embodiment the spectral estimation comprises a Chirp Z transform that provides 
20 a frequency spectrum of heart rate information. The Chirp Z transform is used rather than a conventional 
Fourier Transform because e frequency range for the desired output can be designated in a Chirp Z 
transform. Accordingly, in the present embodiment a frequency spectrum of the heart rate is provided 
between 30 and 250 beats/minute. In the present embodiment the frequency spectrum is provided to • 
spectrum analysis modute 590 which selects the first harmonic from the spectrum as the pulse rate. 
25 Usually, the first harmonic is the peak in the frequency spectrum that has the greatest magnitude and 
represents the pulse rate. However, in certain conditions, the second or third harmonic can exhibit the 
greater magnitude. With this understanding, in order to select the first harmonic the first peak that has 
an amplitude of at least l/20th of the largest peak in the spectrum is selected). This minimizes the 
possibility of selecting es the heart rate a peak in the Chirp Z transform caused by noise. 
30 In the case of motion, a motion artifact suppression is completed on the snapshot with the motion 

artifact suppression module 580. The motion artifact suppression module 580 is depicted in greater detai 
in Figure 21. As can be seen in Figure 21, the motion artifact suppression module 580 is nearly identical 
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to the saturation transform module 406 (figure 1* Accordingly, the motion artifact suppression module 
has a motion artifact reference processor 570 and a motion artifact correlation canceler 571. 

The motion artifact reference processor 570 is the same as the reference processor 530 of the 
aaturetion transform module 406. However, the reference processor 570 utifces the saturation value from 

5 the saturation module 408. rather than completing an entire saturation transform with the 117 saturate 
scan values. The reference processor 570, therefore, has a saturation equation module 581, a reference 
generator 582, a DC removal module 583, and a bandpass titer module 585. These modules are the same 
as corresponding modules in the saturation transform reference processor 530. In the present embodiment, 
the saturation equation module 581 receives the arterial saturation value from the saturation ce«culat,on 

10 module 408 rather than doing a saturation axis scan as in the saturation transform module 406. Tms « 
because the arterial saturation has been selected, end there is no need to perform an a«s scan. 
Accordingly, the output of the saturation equation module 581 corresponds to the proportionality constant 
r He the expected red to infrared retio for the arterial saturation value). Otherwise, the reference 
processor 570 performs the same function as the reference P^or 530 of the saturation transform 

15 module 406. 

The motion artifact correlation canceler 571 is also similar to the saturation transform correlate 
canceler 531 (Figure 18). However, the motion artifact suppression correlation canceler 571 uses a slightly 
different motion artifact joint process estimator 572. Accordingly, the motion artifact suppress^ 
correlation canceler 571 has a joint process estunator 572 and a low-pass filter 573. The motion artifact 
20 joint process estimator 572 differs from the saturation transform joint process estimator 550 in that there 
are a different number of ceto (between 6 and 1 10 in the present embodiment), as selected by the Number 
of Cens input 574. m that the forgetting parameter differs W in the present embodiment), and a, that 
the time delay due to adaptation differs. The low-pass filter 573 is the same as the low pass futer 552 
of the ssturation transform correlation canceler 531. 
25 Beceuse only one saturation value is provided to the reference processor, only one output vector 

of 270 samples results at the output of the motion artifact suppression correlation canceler 571 for each 
input packet of 570 samples, hi the present embodiment, where the infrared wmrelength is provuled as 
a first input to the correlation canceler. the output of the correlation canceler 571 provides a dean infrared 
wmref orm. It should be understood that, a, described above, the infrared end red wavelength signals could 
30 be switched such that a clean red waveform is provided at the output of the motion artifact suppression 
correlation canceler 571. The output of the correlation cancel 571 is . clean waveform because the 
actua. saturation value of the patient is known which allows the reference processor 570 f generate - 
noise reference for inputting to the corre.af.on canceler 571 as the reference signal. The clean waveform 
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tt the output of the motioii artifact suppression modulo 580 is o clean plethysmograph waveform which 
Ccii be forwarded to the display 336. 

As described ebove, an alternative joint process estimator uses the QRD least squares lattice 
pproach (Figure 8a, 9a and 10a). Accordingly, the joint process estimator 573 (as wel as the joint 
process estimator 550) could be replaced with a joint process estimator executing the QRD least squares 
lattice operation. 

Figure 21a depicts an alternative embodiment of the motion artifact suppression module with a joint 
process estimator 572a replacing the joint process estimator 572. The joint process estimator 572a 
comprises a QRD least squares lattice system as in Figure 10a. In accordance with this embodiment, 
different initialization parameters are used as necessary for the QRD algorithm. 

The initialization parameters are referenced in Figure 21a as -Number of CeOs." "Lambda/ 
-MinSumErr." "Gamslnit" and "SumErrmit." Number of Cells and Lambda correspond to Eke parametefs 
in the joint process estimator 572. Gamslnit corresponds to the y initialization variable for aR stages 
except the zero order stage, which as set forth in the QRD equations above is initialized to 'V. 
SummErrlnit provides the 6 initialization parameter referenced above in the ORD equations. In order to 
avoid overflow, the larger of the actual calculated denominator in each division in the QRD equations and 
MinSumErr is used. In the present embodiment, the preferred initialization parameters are as follows: 

Number of Cels - 6 

Lambda - .8 

MinSumErr - 10* 20 

Gamslnit - 10* Z 

SumtrrlnH - Iff 8 . 

The dean waveform output from the motion artifact suppression module 580 also provides an input 
to the second spectral estimation module 588. The second spectral estimation module 588 performs the 
same Chirp Z transform as the first spectral estimation module 586. In the case of no motion, the output 
from the first spectral estimation module 586 is provided to the spectrum analysis module 586; in the case 
of motion, the output from the second spectral estimation module 588 is provided to e spectrum analysis 
module 590. The spectrum analysis module 5S0 examines the frequency spectrum from the appropriate 
spectral estimation module to determine the pulse rate. In the case of motion, the spectrum arialysis 
module 590 selects the peak in the spectrum with the highest amplitude, because the motion artifact 
suppression module 580 attenuates al other frequencies to s value below the actual heart rate peak. In 
the case of no motion, the spectrum analysis module selects the first harmonic m the spectrum as the heart 
rate as descr&ed above. 
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The output of the spectrum analysis module 590 provides the raw heart rate as an input to the 
slew rate limiting module 592. which provides an input to an output filter 594. In the present embodiment 
the slew rete firniting module 592 prevents changes greater that 20 beats/minute per 2 second interval 

5 The output filter 594 comprises en exponential smoothing filter similar to the exponential smoothing 

filter described ebove with respect to the clip and smooth filter 566. The output filter is controlled via 
an output filter coefficient module 596. If motion is large, this filter is slowed down, if there is little or 
no motion, this filter can sample much faster and still maintain a clean value. The output from the output 
filter 594 is the pulse of the patient, which is edvantageously provided to Uie display 336. 

10 Alte rnative To Saturation Tr ansform Mndufr-Benk Of FlttBl 

An alternative to the saturation transform of the saturation transform module 406 can be 
implemented with a bank of filters as depicted in figure 23. As seen in figure 23. two benks of filters, 
a first filter bank 600 and a second filter bank 602 are provided. The first filter bank 600 receives a f,rst 
measured signal S^t) (the infrared signal samples in the present embodiment) on a corresponding first filter 

15 bank input 604, and the second filter bank 602 receives a second measured signel S^t) (the red samples 
ft, the present embodiment) on a corresponding second filter bank input 606. In a preferred embodiment, 
the first and second filter banks utilira static recursive polyphase bandpass filters with fixed center 
frequencies and comer frequencies. Recursive polyphase filters are described in an article Harris, et el 
-Digital Signal Processing With Efficient Polyphase Recursive All-Pass filters" attached hereto es Append,* 

20 A However, adaptive implementations are also possible. In the F «ent implementation, the recursive 
polyphase bandpass filter elements are each designed to indude a specific center frequency and bandwidth. 

There are N filter elements in each filter bank. Each of the fitter elements in the first filter bank 
600 have e matching (i.e. same center frequency end bandwidth) filter element in the second filter bank 

25 602. The center frequencies and the comer frequencies of N elements are each designed to occupy N 
frequency ranges, 0 to F,. F, - F 2 , F 2 - F 3 . F 3 - F 4 . ... F^ • F„ as shown in figure 23. 

ft should be understood that the lumber of filter elements can range from 1 to infinity. However, 
in the present embodiment, there ere approximately 120 separate filter elements with center frequencies 
spread evenly across a frequency range of 25 beats/inmute - 250 Ijeatsmiimrte. 

30 The outputs of the filters contain information about the primary and secondary signals for the first 

and secon d measured signals (red end infrared in the present example) at the specified frequencies. The 
outputs for each oak of matching filters (one in the first filter bank 600 and one in the second filter bank 
602) are provided to saturation determination modules 610. figure 23 depicts only one saturation 
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determination module 610 for esse of illustration. However, a saturation determination module can be 
provided for each matched pair of filter dements for parallel processing. Each saturation determination 
module has a ratio module 616 and a saturation equation module 818. 

The ratio module 616 forms a ratio of the second output to the first output For instance, in the 

5 present example, a ratio of each red RMS value to each corresponding infrared RMS value (Red/IR) is 
completed in the ratio module 616. The output of the ratio module 616 provides an input to the saturation 
equation module 618 which references a corresponding seturation value for the input ratio. 

The output of the seturation equation modules 618 are collected (as represented in the histogram 
module 620) for each of the matched filter pairs. However, the data collected is initially a function of 

10 frequency and saturation. In order to form a saturation transform curve similar to the curve depicted h 
Figure 22, a histogram or the like is generated as in Figure 24. The horizontal axis represents the 
saturation value, and the verticel exis represents a summation of the number of points (outputs from the 
saturation equation modules 618) collected at each saturation value. In other words, if the output of the 
saturation equation module 618 for ten different matched filter pairs indicates a saturation value of 98%, 

IS then a point in the histogram of Figure 24 would reflect a value of 10 at 98% saturation. This results 
in s curve similar to the saturation transform curve of Figure 22. This operation is completed in the 
histogram module 620. 

The results of the histogram provide a power curve similar to the power curve of Figure 22. 
Accordingly, the arterial saturation can be calculated from the histogram by selecting the peak (greatest 
20 number of occurrences in the area of interest) cwrespondkig to the highest saturation value (e.g, the peek 
V in Figure peaks corresponding to the highest saturation value peak. Similarly, the venous or background 
saturation can be determined from the histogram by selecting the peak corresponding to the lowest 
saturation value (e.g, the peak 'd' in Figure 24). in e manner similar to the processing in the saturation 
calculation module 408. 

25 It should be understood that as an alternative to the. ' istogram, the output saturation (not 

necessarily a peak in the histogram) corresponding to the highest;*.' turation value could be selected as the 
arterial saturation with the corresponding ratio representing r r Similarly, the output saturation 
corresponding to the lowest seturation value could be selected as the venous or background saturation with 
the corresponding ratio representing r r For example, in this embodiment the entry 'a' in the histogram 

30 of Frgure 24 would be chosen as the arterial saturation and the entry in the histogram V ^.h the lowest 
saturation value would be chosen as the venous or background saturation. 
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ALTERNATIVE nETERMINATIQ M HF MEmCIPHTS r n AMD r T 
As explained above, in accordance with the present invention, primary and secondary signal 
portions, particularly for pulse oximetry, can be modeled as follows: 

5 ^ " »1 * "I ™ (89} 

S |R - s 2 + n 2 (infrared) W 

■*1 " r a s 2 9t ^ n l " r v°2 ^ 
Substituting Equation (91) into Equation (89) provides the following: 

Sred " r a»2 + r v"2 (red > (92> 
10 Note that S^ and S,„ are used in the model of equations (89) • (92). This is because the 

discussion below is particularly directed to blood oximetry. S^, and S |R correspond to S, and S 2 in the 
preceding text, and the discussion that follows could be generalized for any measure signal S, and S 2 . 

As explained above, determining r a and r„ (which correspond to arterial and venous blood oxygen 
saturation via a saturation equation) can be accomplished using the saturetion transform described above 
15 doing a scan of many possible coefficients. Another method to obtain r a and r„ based on red and infrared 
data is to look for r a and r„ which minimize the correlation between h and n k , assuming s* is at least 
somewhat (and preferably substantially) uncorrected with n k (where k - 1 or 2). These values can be 
found by minimizing the following statistical calculation function for k - 2: 



20 



(93) 



where i represents time. 

It should be understood that other correlation functions such as a normalized correlation could also 

be used. 

Minimizing this quantity often provides a unique pair of r, and r„ if the noise component is 
25 ^correlated to the desired signal component. Minimizing this quantity can be accomplished by solving 
Equations (90) and (92) for s 2 and n 2 . and finding the minimum of the correlation for possible values of 
r a and r r Solving for s 2 and n 2 provides the fonowing: 



inverting the two-by-two matrix provides: 

30 
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Preferably, the correlation of equation (93) is enhanced with a user specified window function as follows: 

N 
E 



(93a) 



10 



The Btockman Window is the presently preferred embodiment. It should be understood that there 
are many additional functions which rninimize the correlation between signal and noise. The function above 
is simply one. Thus, 
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Correlation n£ = j£[ ^ - r v$m)(- s redi 

(93b) 

In order to implement the minimization on a plurality of discrete data points, the sum of the 
squares of the red sample points, the sum of the squares of the infrared sample points, and the sum of 
5 the product of the red tones the infrared sample points are first calculated deluding the window function, 



RR - hSredfWi 

M 

IRR - ^(S lf i(S red )w i 

M 

These values are used in the correlation equation |93b). Thus, the correlation equation becomes 
10 en equation m terms of two variables, r, and r r To obtain r a and r r an exhaustive scan is executed for 
a good cross-section of possible values for r, and r v <e. B .. 20-50 values each corresponding to saturate 
values ranging from 30 - 105% The minimum of the correlation function is then selected and the values 
of r a and r v which resulted in the minimum are chosen as r, and r r 
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Once r B and r v have been obtained, arterial oxygen saturation and venous oxygen saturation can 
be determined by provided r a and r f to a saturation equation, such as the saturation equation 502 of the 
statistics module 404 which provides an oxygen saturation value corresponding to the ratios r a and r r 

In a further implementBtion to obtain r, end r,, the same signal model set forth above is again 
used. In order to determine r t and r ¥ in accordence with this implementation, the energy in the signal Sj 
is inaxirnized under the constraint that Sj is uncorrelated with n 2 . Again, this implementation is based upon 
minimizing the correlation between s and n and on the signal model of the present invention where the 
signal s relates to the arterial pulse end the signal n is the noise (containing information on the venous 
blood, as wefl as motion artifacts and other noise); r, is the ratio (RED/1R) related to arterial saturation and 
r, is the ratio (REDIIR) related to venous saturation. Accordingly, 'm this implementBtion of the present 
invention, r a end r ¥ are determined such that the energy of the signal Sj is meximized where s 2 end n 2 
ere uncorrelated. The energy of the signal Sj is given by the following equation: 



1 N 

ENERGY = — (S^-r^ 2 



— 1— \R^rJR^r 2 R^ 



(94) 



(95) 



(95) 



where R t is the energy of the red signal, R 2 is the energy of the infrared signal end R u is the correlation 
between the red and infrared signals. 

The correlation between. s 2 and n 2 is given by 

- — ^— [-tfi+('V^)*i.2- / V r <A] < 9 * 
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1 N O 

Correlation{s 2 ,n^= ^(S^-r^^-S^+r^; 

As explained above, the constraint is that \ and m ft - 2 for the present example) are 
uncorrelated. This "decorrelation constraint" is obtained by setting the correlation of Equation (97) to zero 
as follows: 



5 In other words, the goal is to maximize equation (94) under the constraint of equation (98). 

In order to obtain the goal, a cost function is defined (e.*. a Lagrangian optimization in the present 
embodiment) as follows: 

(r a -r v ) 2 

(89) 



where fi is the Lagrange multiplier, finding the value of r,, r ¥ and // that solve the cost function can be 
10 accomplished using a constrained optimizetion method such as described in luenberger. linear & Nonlinear 
Programming. Addison-Wesley, 2d E4. 1984. 

Along the same lines, if we assume that the red and infrared signals and S IR are non-static 
the functions R,. R 2 and R 12 defined above era time dependent Accordingly, with two equations, two 
unknowns can be obtained by expressing the decorrelation constraint set forth in equation (98) et two 
different times. The decorrelation constraint can be expressed at two different times, f, and f 2 . as 



15 



follows: 
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Because equations (100) and (101) are non-linear in r, and r r a change of variables allows the use of 
linear techniques to solve these two equations. Accordingly, with x - r, ♦ y y - r,r v equations (100) 
and (101) become 



/? 12 (yx - /? 2 (yy - 



(103) 



5 These equation (102) and (103) can be solved for x and y. Then, solving for r a and r ¥ from the changes 
of variables equations provides the following: 

r + 2- = x rl - xr v +y = O oo4> 
W 

Solving equation (104) results in two values for r r In the present embodiment, the r ¥ value that results 

in x 2 -r„y >° « selected - ,f both vaIues of r « result m x2 f v y> °' the r w th8t mMimizM thB merm of *2 
10 (Enerjjjsj)) at f 2 is selected. r ¥ is then substituted into the equations above to obtain r r Alternatively 
r B can be found directly in the same manner r t was determmed. 
A hTn„tiv Ta Satura tinn Transform-Complex FFT 

The Mood oxygen saturation, pulse rate and a clean plethysmography waveform of a patient can 
ebo be obtained using the signal model of the present invention using a complex FFT. as explained further 
15 with reference to Figures 25A-25C. In general by utilizing the signal model of equations (89M921 with 
two measured signals, each with a first portion and a second portion, where the first portion represents 
i desired portion of the signal and the second portion represents the undesired portion of the signal, and 
where the measured signals can be correlated with coefficients r, and y a fast saturation transform on 
a discrete basis can be used on the sample points from the output of the decimation operation 402. 
20 Rgure 25A corresponds generally to Figure 14, with the fast saturation transform replacing the 

previously described saturation transform. In other words, the operations of Rgure 25A can replace the 
operations of Rgure 14. As depicted in Rgure 25A. the fast saturation transform is represented in a fast 
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saturation transforrrjpulse rate calculation module 630. As in Figure 14. the outputs are arterial oxygen 
saturation. . clean plethysmography waveform, and pulse rate. Hour. S3 and 25C illustrate additional 
detail regarding the fast saturation ttansform/pulse rateLalculation n:doH 630. As depicted in Figure 25B, 
the fast saturation transform module 630 has infrared log and red log modules 640, 642 to perform a log 

5 normalization as m the infrared and red log modules 480. 482 of figure 17. Similarly, there are infrared 
DC removal and red DC removal modules 644, 646. In addition, there are infrared and red high-pass filter 
modules 645, 647, window function modules 648. 640. complex FFT modules 652. 654. select modules 
653, 655. magnitude modules 656. 658. threshold modules 660. 662. a point-by-point ratio module 670, 
a saturation equation module 672. end a select saturation module 680. There are also phase modules 690. 

10 692. a phase difference module 694, and a phase threshold module 696. The output of the select 
saturation module 680 provides the arterial saturation on an arterial saturation output line 682. 

In this alternative embodiment, the snapshot for red end infrared signals is 562 samples from the 
decimation module 402. The infrared DC removal module 644 and the red DC removal module 646 are 
slightly different from the infrared and red DC removal modules 484, 486 of figure 17. In the infrared 

15 and red DC removal modules 644. 646 of figure 25B, the mean of all 563 sample points for each 
respecth* channel is calculated. This mean is then removed from each individual sample point m the 
respective snepshot in order to remove the baseline DC from each sample. The outputs of the infrared and 
red DC removal modules 644. 646 provide inputs to respective infrared high-pass filter module 645 and rod 
high-pass filter module 647. 

20 The high-pass filter modules 645. 647 comprise FIR filters with 51 taps for coefficumts. 

Preferably, the high-pass filters comprise Chebychev filters with a side-lobe level parameter of 30 and . 
comer frequency of .5 Hz (Le., 30 beats/minute). It will be understood that this filter could be vaned for 
performance. With 562 sample points entering the high-pass filters, and with 51 taps for coefficients, 
there are 512 samples provided from these respective infrared and red snapshots et the output of the high. 

25 pass filter modules: The output of the high-pass filter modules provides an input to the window function 
modules 648. 650 for each respective charmd. 

The window function modules 648, 650 perform a conventional windowing function. A Kaiser 
windowing function is used in the present embodiment The functions throughout figure 25B maintain a 
point-by-point analysis. In the present embodiment the time bandwidth product for the Kaiser window 

30 function is 7. The output of the window function modules provides an input to the respective complex Fast 
Fourier Transform {FFT) modules 652, 654. 

The complex FFT modules 652. 654 perform complex FFTs on respective infrared and red channels 
on the data snapshots. The data from the complex FFTs » then analyzed in two paths, once which 
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examines the magnitude and one which examines the phase from the complex FFT data points. However, 
prior to further processing, the date is provided to respective infrared and red select modules 653, 655 
because the output of the FFT operation wffl provide repetitive information from 0-1/2 the sampling rate 
end from 1/2 the sampling rate to the sampling rate. The select modules select only samples from 0-1/2 

5 the sampling rate (e.g., 0 31.25 Hz in the present embodiment! and then select from those samples to cover 
a frequency renge of the heart rate and one or more harmonics of the heart rate. In the present 
embodiment samples which fall in the frequency range of 20 beats per minute to 500 beats per minute 
ere selected. This value can be veried in order to obtain harmonics of the heart rate as desired. 
Accordingly, the output of the select modules results in less than 258 samples. In the present embodiment. 

10 the sample points 2-68 of the outputs of the FFTs ere utilized for further processing. 

In the first path of processing, the output from the select modules 653, 655 are provided to 
respective infrered and red magnitude modules 656, 658. The megnitude modules 656. 658 perform • 
magnitude function wherein the magnitude on a point-by-point basis of the complex FFT points is selected 
for each of the respective channels. The outputs of the magnitude modules 656. 658 provide an input to 

15 infrared and red threshold modules 660, 662. 

The threshold modules 660. 662 exemine the sample points, on a point-by-point basis, to select 
those points where the magnitude of en individual point is above e particular threshold which is set at a 
percentage of the meximurn magnitude detected among all the remaining points in the snapshots. In the 
present embodiment the percentage for the threshold operation is selected as 1% of the maximum 

20 magnitude. 

After thresholding, the dete points are forwarded to a point-by point ratio module 671 The point- 
by-point retio module tekes the red over infrared ratio of the values on s poirrt-by-point basis. However, 
a further test is performed to quaDfy the points for which a ratio is taken. As seen in Figure 25B, the 
sample points output from the select modules 653. 655 are also provided to infrared and red phese modules 

25 690. 692. The phase modules 690. 692 select the phase value from the complex FFT points. The output 
of the phase modules 690, 692 is then presented to a phase difference moduto 694. 

The phese difference module 694 calculates the difference in phase between the corresponding 
data points from the J phi: e modules 690, 692. If the magnitude of the phase difference between any two 
corresponding points' is less than a particular threshold (e.g. 0.1 radians) in the present enAodiment). then 

30 the sample points qualify. If the phase of two corresponding sample points is too far apart then the 
sample points ere not used. The output of the phase threshold module 696 provides en enable input to 
the REDJIR rate module 670. Accordingly, in order for the retio of e particular pair of sample points to 
be taken, the three tests are executed: 
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6. the red sample must pass the red threshold 660; 

7. the infrared sample must pass the Infrared threshold 662; and 

8. the phase between the two points must be less than the predefined threshold as 

determined in the phase threshold 696. 

5 For those sample points which qualify, a ratio is taken in the ratio module 670. For those points 

which do not qualify, the saturation is set to zero at the output of the saturation equation 672. 

The resulting ratios are provided to e saturation equation module which is the seme as the 
saturation equation modules 602, 520 in the statistics module 504. In other words, the saturation equation 
module 672 accepts the ratio on a point-by-point bash and provides as an output a corresponding, saturate 

tO va.ua corresponding to the discrete ratio points. The saturation points output from the saturation equate 
module 672 provide a series of saturation point, which could be plotted as saturation with respect to 
frequency. The frequency reference was entered into the points at the complex FFT stage. 

The arterial (and the venous) saturation can then be selected, as represented in the select arterial 
saturation module 680, in one of two methods according to the present invention. According to one 

15 method, the arterial saturation value can be selected simply as the point corresponding to the largest 
seturation vahae for an pomts output from the saturation equation module 672 for a packet Artematrvely. 
a histogram similar to the histogram of figure 22 <*n be generated in which the number of saturation 
values at different frequencies (points) are summed to form a histogram of the number of occurrences for 
each particular saturation value. In either method, the arterial saturation can be obtained end provided as 

20 an output to the select arterial saturation module on the arterial saturation output fine 682. In order to 
obtain the venous saturation, the rmrdmum arterial saturation value, of points that exhibit non-zero value, 
is selected rather than the maximum arterial saturation value. The saturation can be provided to the 

display 336. ' . 

The fast saturation transform information can also be used to provide the pulse rate and the clean 

25 plethysmographic wave form as further illustrated in Figure 25C. In order to obtain the pulse rate and a 
clean ptethysmographic wave form, several additional function, are necessary. As seen in Figure 25C, the 
pulse rate and clean plethysmography wave form are determined using . window function module 700. a 
spectrum analysis module 702 and an inverse window function module 704. 

As depicted in figure 25C. the input to the window function module 700j ,b opined from the 

30 output of the complex FFT module, 652 or 654. In the present embodiment, only one ensured signal ,s 
necessary. Another input to the window function module 700 is the arterial saturation obtained from the 
output of the select erterial saturation module 680. 
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Tha window function module perfonns a windowing function selected to pass those frequencies 
that significantly correlate to the frequencies which exhibited saturation values very close to the arterial 
saturation value. In the present embodiment, the following windowing function is selected: 



1- 



SAT^-SAT, 



115 



(105) 



100 



10 



15 



where SAT„ equals the saturation value corresponding to each particular frequency for the sample points 
and SAT ^represents the arterial saturation as chosen at the output of the select arterial saturation 
module 680. This window function is applied to the window function input representing the complex FFT 
of either the red or the infrared signal. The output of the window function module 700 is a red or infrared 
signal represented with a frequency spectrum as determined by the FFT, with motion artifacts removed by 
the windowing function. It should be understood that many possible window functions can be provided. 
In addition, with the window function described above, it should be understood that using a higher power 
will provide more noise suppression. 

In order to obtain pulse rate, the output points from the window function module 700 are provided 
to a spectrum analysis module 702. The spectrum analysis module 702 is the seme as the spectrum 
analysis module 590 of Figure 20. In other words, the spectrum analysis module 702 deternines the pulse 
rate by deterrnining the first harmonic in the frequency spectrum represented by the output points of the 
windowing function 700. The output of spectrum analysis module 702 is the pulse rate. 

In order to obtain a clean plethysmography waveform, the output of the windowing function 700 
i, applied to an inverse window function module 704. The inverse window function module 704 completes 
an inverse of the Kaiser window function of the window function module 648 or 650 of Figure 25B. In 
other words, the inverse window function 704 does a pomt-bypoint inverse of the Kaiser function for 
points that are stiB defined. The output is a dean ptethysmographic waveform. 

Accordingly, by using a complex FFT and windowing functions, the noise can be suppressed from 
the ptethysmographic waveform ki order to obtain the arterial saturation, the pulse rate, and a deaf. 
25 ptethysmographic waveform. It should be understood that although the above description relates to 
operations prirnarify in the frequency domain, operations that obtain similar results could else be 
accompished in the time domain. 
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(105c) 



-piatinn to P ffrtfr" 1 E°"ations 

The measurements described for pulse oximetry above are now related back to the more 
generalized discussion above. The signals (logarithm converted) transmitted through the finger 310 at each 
wavelength Aa end Ab ere: 

5 

S / , > (t)-S^ redl (t)-eHb02V A HhO2 ,tA ^ t > +e Hb>la cA Hb xA^t, 

W^vW^ . (105a) 

10 S yl8 (t)-f H b02>l. cA Hb02 3tA W +€ HbMa cA Hb xA W * 

SylaW" *A* {%) * n ^« (t,; 
S JUl (t)-S y , red2 W-e HD0W b cA Hh02 JtA W +e Hb^b cA Hb xA W 

•^-••ii^bAi^* l10BW 

(106c) 

20 S lh (t)- s, h (t) ♦ ni b (t) . 

The variables above are best understood as correlated to figure 6c as follows: assume the layer m figure 
6c containing A 3 and A 4 represents venous blood in the test medium with A3 representmg deoxygenet«l 
hemoglobin (Hb) and A 4 representing oxygenated hemoglobin (HB02) in the venous blood. Simflarfy. assume 
that the layer in Hgure 6c contenting Ag end A 6 represents arterial blood in the test medium, wrth A* 
25 representing deoxygenated hemoglobin (Hb) end A 6 representing oxygenated hemoglobin (HB02) m the 
arteriel blood. Accordingly. ^02 represents the concentration of oxygenated hemoglobin in the venous 
blood. cVHb represents the concentration of deoxygenated hemoglobin in the venous blood, x represents 
the thickness of the venous blood (e.g. the thickness the layer contouring A3 and A 4 ) Similarly, c Hb02 
represent, the concentration of oxygenated hemoglobin to the erterial blood. c*Hb represents the 
30 concentration of deoxygeneted hemoglobin in the erterial blood, and x* represents the ttackness of the 
arterial blood (e.g- the thickness of the layer containing A5 and Ag) 

The wavelengths chosen ere typically one in the visible red range, U. An. end one in the nfrared 
range. Le..>lb. Typical wavelength vah.es chosen ere Am - 660nm and Am - 910nm. m accordance with 
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the constant saturation method, ft is assumed that c^Wlc^t) - consent, and c^WfcVt) - 
constant^ The oxygen saturation of arterial and venous blood changes slowly, if « al. with respect to 
the sample rate, making this a valid assumption. The proportionaBty coefficients for equations (105) and 
(106) can then be written as 



r ^ = G HbdZXa C HbO^ +*Hb.\a C HtA f ) (107 ) 



W» ' '.WW** 
10 * 'a^iW 

s M [t) * r w (t) $Jlb (t) 



15 



20 



25 



(108a) 
(108a) 
(108b) 
(109b) 



in pulse oximetry, it is typically the cue that both equations (108) and (109) can be satisfied 

simultaneously. 

Multiplying equation (106) by r,(t> and then subtract** equation (106) from equation (105), a non- 
zero secondary reference signal n'(t) is determined by: 

n'(t) - S^tt) - r.WS^tt) l"0«) 



Multiplying equation (106) by r^O and then subtracting equation (106) from equation (105). a non-iero 

prinarv reference signal s'(t) is determined by: 

sW-S^W-r^W < 110W 
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The constant saturation assumption does not cause the venous contribution to the absorption to 
be canceled along with the primary signal portions s„,tt) and s^t). Thus, frequencies associated with 
both the low frequency modulated absorption due to venous absorption when the patient is still and the 
modulated absorption due to venous absorption when the patient is moving are represented in the secondary 
5 reference signal n'tt). Thus, the correlation canceter or other methods described above remove or derive 
both erratically modulated absorption due to venous blood in the finger under motion and the constant low 
frequency cyclic ebsorption of venous blood. 

To illustrate the operation of the oximeter of Figure 11 to obtain clean waveform, figures 26 and 
27 depict signals measured for input to a reference processor of the present invention which employs the 
10 constant saturation method. Le.. the signals S^W - S At Jt) and S^ft) - S^ IR (t|. A first segment 26a 
end 27a of each of the signals is relatively undisturbed by motion artifact. Le, the patient did not move 
' substantially during the time period in which these segments were measured. These segments 26a and 
27a are thus generally representative of the primary plethysmogrephic waveform at each of the measured 
wavelengths. A second segment 26b and 27b of each of the signals is affected by motion ertifect Le.. 
15 the patient did move during the time period in which these segments were measunML Each of these 
segments 26b end 27b shows large motion induced excursions in the measured signal. A third segment 
26c and 27c of each of the signals is again relatively unaffected by motion artifact and is thus generally 
representative of the primary plethysmogrephic waveform at each of the measured wavelengths. 

figure 28 shows the secondary reference signal n'tt) - n^t) • r^ft). as determined by a 
20 reference processor of the present invention. Again, the secondary reference sigmri n'tt) is correlated to 
the secondary signal portions n,. and n>Ur Thus, a first segment 28. of the secondary reference signal 
n'tt) is generally flat, corresponding to the fact that there is very little motion induced muse m the first 
segments 26a and 27a of each signaL A second segment 28b of the secondary reference signal n'tt) 
exhibits large excursions, corresponding to the large motion induced excursions in each of the measured 
25 signals. A third segment 28c of the noise reference signal n'W is generally flat again corresponding to 
the lack of motion artifact in the third segments 26c end 27c of each measured signaL 

It should also be understood that a reference processor could be utfiHed in order to obtain the 
primary referent signel s'(t) - «^ A bW- ^ primary reference signal s'(t) would be generaBy indiartive 
of the plethysmogrsph waveform. 
30 figures 29 and 30 show the approximations s^tt) and s-^tt) to the primary signals s^tt) and 

s^ltl as estimated by a correlation cancder using a secondary reference signal n'tt). Note that the scale 
of Figures 26 through 30 is not the same for each figure to better illustrate changes, in each signaL 
figures 29 and 30 illustrate the effect of correlation cancellation using the secondary reference signal n'tt) 
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as determined by the reference processor. Segments 29b and 30b are not dominated by motion induced 
noise as were segments 26b and 27b of the measured signals. Additionally, segments 29a, 30a, 29c and 
30c have not been substantially changed from the measured signal segments 26a, 27a, 26c, and 27c where 
there was no motion induced noise. 
S It should be understood that approximation n^,(t) and n% b (t) to the secondary signals n A JA and 

n M M es estimated by a correlation canceler using a primary reference signal s'lt) can also be determined 

in accordance with the present invention! 

MFTHflD FOR ESTtMftTIMB PRIMARY ™T\ SFMUDARY SlfiNAl PORTIONS OF 
MFASIIRFD SIGHft l^ W A pl " Sg OXIMETER 

10 Implemeriting the various embodiments of the correlation canceler described above in software is 

relatively straightforward given the equations set forth above, and the detailed description above. However, 
e copy of a computer program subroutine, written in the C programming language, which calculates a 
primary reference s'lt) using the constant saturation method and. using a joint process estimator 572 which 
implements a joint process estimator using the equations (54) • (64) is set forth in Appendix B. This joint 

15 process estimator estimates a good approximation to the primary signal portions of two measured signals, 
each having a primary portion which is correlated to the primary reference signal s'(t) and a secondary 
portion which is correlated to the secondary reference signal n'(t). This subroutine is another way to 
implement the steps illustrated in the flowchart of Figure 9 for a monitor particularly adapted for pulse 
oximetry. The two signals are measured at two different wavelengths A* and Ab. where At is typically 

20 in the visible region and Ah is typically in the infrared region. For example, in one ernbodknent of the 
present invention, tailored specifically to perform puke oximetry using the constant saturation method. Am 
- 660nm and Ah - 940nm 

The correspondence of the program variables to the variables defined in equations (54) . 164) in 

the discussion of the joint process estimator is as follows: 

25 

A^ft) - ntfrnlDelta 

r tm W " IM * n l- fref 
r b Jt) - nclmJJjref 

f m (t) - nclmlferr 

30 bjt) - ntimlberr 

QJti - itdml.Fswsqr 

B^ft) - nclmLBswsqr 

K ro W - ndmlGamma 
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PnjJA - nclm].Roh_a 
PnUbW " nc N.Roh_b 
e^gW - nc[m].BTT_8 

5 *nvlB W - "W-M 

*nvlb w " n rt n l- K - b 

A first portion of the program performs the initialization of the registers 90. 92. 96. and 98 and 
intermediate variable values as in the "INITIALIZED CORRELATION CANCELED action block 120. A second 

,0 portion of the program performs the time updates of the delay element variables 1 10 with the value at the 
input of each delay element variable 110 is stored in the delay element variable 110 as in the "TIME 
UPDATE OF LEFT IZ' 1 l ELEMENTS" action block 130. The calculation of saturation is performed m a 
separate module. Various methods for calculation of the oxygen saturation are known to those skilled in 
the ert One such calculation is described in the articles by G.A. Mook. et .1 and Michael R. Neuman cited 

,5 above Once the concentration of oxygenated hemoglobin and deoxygenated hemoglobin are determined, 
the value of the saturation is determined similarly to equations (721 through (79) wherein measurements 
at times t, and t 2 are mede at different, yet proximate times over which the saturation is relate 
constant For puke oximetry, the average saturation at time t - (t 1+ t 2 W2 is then determined by: 

<a _ CHboAfi (112a) 

^hrteriafJ) A /A A /A 



e H B*a-eHbOZXa<eHB,Xb^Hb02.Xb)(* S J* S *b) 



12b) 



(A = C HbodA (113a) 

CHbOiM * C Hb\V 
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^Hb.Xa^Hb.Xb^ n xJ^ n Xb) J* 

*HB,Xa-*Ht02,Xa-^HB.Xb- e Ht>O2.Xb)(* n J* n \b) 



13b) 



A third portion of the subroutine calculates the primary reference or secondary reference, as in the 
-CALCULATE PRIMARY OR SECONDARY REFERENCE (s'(t> or n'(t» FOR TWO MEASURED SIGNAL 
SAMPLES' ection block 140 for the signals S^t) and S^t) using the proportionality constants r,(t) end 
r (t) determined by the constant saturation method as in equation (31 The saturation is calculated in a 
S separate subroutine and a value of r,W or r^t) is imported to the present subroutine for estimating either 
the primary portions s^.W and s^ft) or the secondary portions n^W and n Ab [t) of the composite measured 

signets S^lt) and S^ft). 

A fourth portion of the program performs Z-stage update as in the "ZERO STAGE UPDATE" action 
block 150 where the Z stagB forward prediction error F 0 (t) and Z-stage backward prediction error b 0 (t) are 
10 set equal to the value of the reference signal n'Jt) or s'W just calculated. Additionally zero-stege values 
of intermediete variables 9f 0 and 0 o ltMnclro].Fswsqr and nclm].Bswsqr in the program! are calculated for 
use in setting registers 90. 92. 96. and 98 values in the least-squares lattice predictor 70 in the regression 
filters 80a and 80b. 

A fifth portion of the program is an iterative loop wherein the loop counter, M. is reset to zero 
15 with a maximum of m-NC.CELLS. as in the "m-O" action block 160 in Figure 9. NC_ CELLS is a 
predetermined maximum value of iterations for the loop. A typical value for NC.CELLS is between B and 
10. for example. The conditions of the loop ere set such that the loop iterates e minimum of five times 
and continues to iterate until a test for conversion is met or m-NC_ CELLS. The test for conversion is 
whether or not the sum of the weighted sum of four prediction errors plus the weighted sum of backward 
20 prediction errors is less than a smal number, typically 0.00001 (Le.. 9Jt)+#n(t) sJ 0JW0O1). 

A sixth portion of the program calculates the forward and backward reflection coefficient T^tJ 
and r^W register 90 and 92 values (nclm].fref and nclmJJiref in the program) as m the "ORDER UPDATE 
m m -STAGE OF LSL-PREOICTOR - action Mock 170. Then forward and backward prediction errors f ro (t) and 
b^t) mdmlferr and nclmlberr in the program) are calculated. Additionally, mtermedtote variables SJti. 
25 /ftnffl. and *t) (nclmLFswsqr, nclmLBswsqr, ndm]. gamma in the program) are calculated. The first cycle 
of the loop uses the value for nclOLFswsqr and nrfOLBswsqr calculated in the ZERO STAGE UPDATE 
portion of the program. 

A seventh portion of the program, stll within the loop begun in the fifth portion of the program, 
calculates the regression coefficient register 96 and 98 values and k^U) WmLKa end 
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„dm].K.b in the program) In both regression filters, as in the "ORDER UPDATE m* STAGE OF 
REGRESSION FILTERIS)" action block 180. Intermediate error signals and variables e^W. e Mb W, 
PnvuM' and 'oUb<« <nclm]^rr_a and ncM.err_b, nclm].roh_a, and ncMrohJ. in the subroutine) are also 

calculated. . 

The loop iterates untfl the test for convergence is passed. The test for convergence of the pant 
process estimator is performed each time the loop iterate, analogously to the "DONE* action block 190. 
If the sum of the weighted sums of the forward and backward prediction errors Qfjt) ♦ M*) i» less than 
or equal to 0.00001. the loop terrrrinetes. Otherwise, sixth and seventh portions of the program repeat 

The output of the present subroutine is a good approximation to the primary signals s'^t) and 
s-.„(t) or the secondary signals n^lt) and n^t) for the set of sample. S A JX) and S Ab M input to the 
program After approximation, to the primary signal portions or the secondary signals portions of many 
sets of measured signal sample, are estimated by the joint process estimator, a compilation of the outputs 
provides waves which are good approximations to the plethysmograpnic wave or motion artifact at each 

wavelength, M and >lb. • . . 

It should be understood that the subroutine of Appendix B is merely one embodiment when 
implements the equetions (54) - (64). Although implementation of the normalized and QRD-LSL equation, 
i, also straightforward, a subroutine for the normalized equations is attached a. Appendix C. end . 
subroutine for the QRO-LSL algorithm is attached es Appendix D. 

While one embodiment of a physiological monitor incorporating a processor of the present invent™ 
for determining e reference signal for use in a correlation canceler, such as an adaptive noise canceler. to 
remove or derive primary and secondary components from a physiological measurement ha. been described 
in the form of a pulse oximeter, ft wfll be obvious to one skilled in the art that other types of physiology 
monitors may also employ the above described techniques 

Furthermore, the signal processing techniques described in the present invention may be used to 
i confute the arterial and venous blood oxygen saturation, of a physiological system on a continuous or 
neerry continuous time be*. These calculation, may. be performed, regardless of whether or not the 
physiological system undergoes voluntary motuui. 

Furthermore, it will be understood that transformations of measured signals other than logarithm* 
conversion end determination of a proportionality factor which allow, removal or derivation of the pranary 
* or secondary signal portions for deterrrunation of a reference signal are possible. AdmtioneOy, although the 
proportionality factor r has been described herein a, a ratio of a portion of . first signal to a portion of 
a second signal, a similar proportionality constant determined as a ratio of a portion of «r second agnd to 



WO 96/12435 PCT/US95/13469 

-89- 

■ portion of a first signal could equally well be utilized in the processor of the present invention. In the 
latter case, a secondary reference signal would generally resemble n'(t) - n^(t) • m^ft). 

Furthermore, it will be understood that correlation cancellation techniques other than joint process 
estimation may be used together with the reference signals of the present invention. These may include 
5 but are not limited to least mean square algorithms, wavelet transforms, spectral estimation techniques, 
neural networks, Weiner and Kabnan filters among others. 

One skilled in the art will realize that many different types of physiological monitors may employ 
the teachings of the present invention. Other types of physiological monitors include, but are in not imited 
to, electro eerdiogrephs, blood ps. ssure monitors, blood constituent monitors (other than oxygen seturetion) 
10 monitors, capnographs. heart rate monitors, respiration monitors, or depth of anesthesia monitors. 
Additionally, monitors which measure the pressure end quantity of a substance within the body such es 
e breathaBzer. a drug monitor, a cholesterol monitor, a glucose monitor, a carbon dioxide monitor, a glucose 
monitor, or a carbon monoxide monitor may also employ the above described techniques. 

Furthermore, one skilled in the art wOl realize that the above described techniques of primary or 
15 secondary signal removal or derivetion from a composite signal including both primary and secondary 
components con elso be performed on electrocardiography (ECG) signals which are derived from positions 
on the body which are close and highly correlated to each other. It should be understood that a tripolar 
Laplacian electrode sensor such es that depicted in Figure 31 which is a modification of a bipolar Laplaeian 
electrode sensor discussed in the erticle "Body Surface Laplacian ECG Mapping- by Bin He end Richerd J. 
20 Cohen contained n the journal IEEE Transactions on Biomedicel Engineering, VoL 39, No. 11, November 
1992 could be used es an ECG sensor. It must also be understood that there ere a myried of possible ECG 
sensor geometry's that may be used to satisfy the requirements of the present invention. The same type 
of sensor could elso be used for EEG and EMG measurements. 

Furthermore, one sklled in the art will realize that the above described techniques can also be 
25 performed on signals made up of reflected energy, rather than transmitted enerflv. One skilled m the art 
wi abo realize that a primary or secondary portion of a measured signal of any type of energy, including 
but not hmited to sound energy. X-ray energy, gamma ray energy, or light energy can be estimated by the 
techniques described above. Thus, one skied in the art wB realize that the techniques of the present 
invention can be epplied in such monitors as those using ultrasound where a signal is trensmitted through 
30 a portion off the body and reflected beck from within the body back through this portion off the body. 
Additionally, monitors such as echo cardiographs may also utilize the techniques of the present invention 
skies they toe rely on transmission and reflection. 
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WhQe the present invention has been described in terms of a physiological monitor, one skilled in 
the art will realize that the signal processing techniques of the present invention can be applied in many 
areas, including but not limited to the processing of a physiological signal The present invention may be 
applied in any situation where a signal processor comprising a detector receives a first signal which includes 
a first primary signal portion and a first secondary signal portion and a second signal which include* . 
second primary signal portion and a second secondary signal portion. Thus, the signal processor of the 
present invention is readily appDcable to numerous signal processing areas. 
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APPENDIX A 
DIGITAL SIGNAL PROCESSING WITH EFFICIENT 
POLYPHASE RECURSIVE ALL-PASS FILTERS 

Fred Harris*, Maximilien d'Oreye de Lantremange*, and A.G. 
Constantinides** 

. Department of Electrical and Computer Engineering San Diego 
State University, San Diego, California 92182-0190, USA. 

* * Siqnal Processing Section, Department of Electrical Engineering 
fmperial College of Science, Technology and Medicine, Exh.bmon 
Road, London DW7-2BT, England. 

ABSTRACT: Digital filters can be realized with remarkably small 
computational burden as polyphase recursive all-pass networks. 
These networks are formed as a sum of M parallel all-pass subfilters 
with phase shifts selected to add constructively in the passband and 
destructively in the stopband. The design technique we present here 
starts with prototype M-path recursive all-pass filters obtained by a 
new optimizing algorithm (described in detail elsewhere). These filters 
are operated with combinations of all-pass transformations, 
resampling, and cascading options. We demonstrate a number of 
designs using the new algorithm atong with examples of filtering 
25 systems obtained by the techniques presented. 



Presented at International Conference on Signal Processing 
Florence, Italy, Sept. 4-6, 1991 
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1. jMTPnni 1CTION 

A standard finite impulse response (FIR) filter can be modeled 
as indicated in Fig. 32, as the weighted summation of the contents of 
a tapped de.ay line. This summation is shown in (1) and the transfer 
5 function of this filter is shown in (2). 



10 
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M-1 ' (1) 

y(n) = £ M m ~ n) x</77 ' 

m-0 



M-1 (2) 

Y[Z\ = h\n) z n 

n*0 



We note the frequency selective characteristics of the Liter 
output is due to the phase shift between successive samples * .the 
tapped deiay Hne. This Phaser summation is shown ,n F,g. 33 

two different frequencies. 

Note that it is the phase shift, and not the we.gh.ng terms, 
which is responsible for the frequency dependent gain of the f..^ 
The weighing terms are used orimariiy to ^T^t-S 
stopband attenuation of the filter. An idea denved from spat,a 
beamforming , k nown as beam spoiHng, (1. and from — 
.peak-to-rms control, <2. is to use structured phase sh.fts <w.th 
arbitrary amplitudes) to accomplish the same goai. 

The class of filters we describe here replaces the we.ghts of the 
FIR flltB r with aii-pass subfiiters which exhibit unity gain w.«h 
Lquency dependent phase shift, .n anticipation of the. polyphase 



SUBSTITUTE SHEET (RULE 26) 



WO 96/12435 PCT/US95/13469 



-93- 

structure, the all-pass subfilters are, as shown in (3) and Fig. 34, first 
order polynomials in Z n where M is the number of taps in the 
structure. This composite form is indicated in Fig. 41a with a 
modified form reflecting a polyphase structure in 41 b. The transfer 
5 function of this structure is shown in (4). 

Figs. 44a and 44b show examples of equal ripple designs 
obtained from this new algorithm for a 2-path and for a 5-path filter 
with two and three all-pass stages per path respectively. Interacting 
constraints on the optimal design (described in the aforementioned 

1 0 paper) limit the possible number of stages per path for the 5-path filter 
to the sequences (1,1,1,1,0), (2,2,2,1,2), (3,3,2.2,2), etc., so that 
the optimum 5-path filter uses three less stages than allocated to the 
design. Similarly the 2-path filter is limited to the sequences (1,0), 
(1,1), (2,1), (2,2), (2,3). (3,3), etc. 

15 2. TWO PATH FILTERS 

We find a wealth of interesting properties and clever applications 
for this structure even if we limit our discussion to a 2-path filter. As 
designed, the 2-path filter is a half bandwidth filter with the 3-dB 
bank-width at 0.25 f_. If the zeros are restricted to the half sampling 

20 frequency, the filter identical to a half bandwidth Butterworth obtained 
by the standard arped bilinear transform. For this special case, the 
real parts of the root locations go to zero. If the zeros are optimized 
for equal ripple stop band behavior the filter becomes a constrained 
Elliptic filter. The constraints are related to a property of 

25 complementary all-pass filters. We define the all-pass sections for the 
2-path filter at H 0 (0) and H t (0) and, as indicated in Fig. 35. the scaled 
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by one-half sum and difference of these paths by A(0) and B(0), the 
lowpass and highpass filters, respectively. We know that the all-pass 
sections satisfy (5) from which we derive the power relationship 
between the lowpass and. highpass filters shown in (6). 

\H 0 (e)\ 2 = \h,(9)\ 2 = 1 (5) 



\AW\ 2 +\B{0)\ 2 = \0.5[H o (e)+H,ie)}\ 2 (6 a) 

+ \o.5iH 0 {e)-H,{e))\ 2 



-o.5[|// 0 («r+i//i(«ii a i. <6b> 



=1 



(6c) 



Now to interpret how this relationship impacts the 2-path filter. 
For the complementary filters we define the minimum passband gain 
of by A-€ U and the peak stopband gain by e 2 . Substituting these gains 
in (6) we obtain (7). 



2_i (7a) 



(1 -*,) 2 =1 -i€ 2 ) 2 (7b) 
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y-2€,+(€y=i-{e z ) 2 {7c) 

For small e v we can '9 nore ^J* which ,eads to (8) " 

5 Thus, if the stopband attenuation is selected to be 0.001 (60.) 

dB), the passband ripple is 0.000 000 5 (1 26.0 dB) or the passband 
minimum gain is 0.999 999 5 (-0.000 000 22 dB). As can be seen 
these filters have a remarkably flat passband. We conclude that the 
2-path equal ripple stopband filters are elliptic filters with coupled 
10 passband and stopband ripple with -3.0 dB gain at 0.25 f . There are 
two significant differences between the 2-path and the standard 
implementation of the equivalent elliptic filter. The first is a four-to- 
one savings in multiplications; a fifth order half bandwidth 2-path filter 
requires only two coefficients as opposed to the direct implementation 

1 5 which requires eight (scaling included). The second is that all-pass 
structures exhibit unity gain to internal states hence do not require 
extended precision registers to store internal states as do standard 

implementations. 

As with any filter design, the transition bandwidth can be traded 
20 for out of band attenuation for a fixed order or transition bandwidth 
can be steepened for a fixed attenuation by increasing filter order. 
While nomographs exist for the elliptic and Butterworth filters which 
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can be implemented by the 2-path structure, a simple approximate 
relationship for the equal ripple filter is given in (9), where A(dB) is 
attenuation in dB, If Af is transition bandwidth, and N is the total 
number of all-pass segments in the filter. 

(9) 

A{dB) = (72Af+10)/V 

2A. HILBERT TRANSFORM FILTERS 
The Hilbert Transform (HT) can be thought of as a particular type 
of half band filter. It is modelled and implemented by a wide band 90' 
phase shifting network. The HT is often used in DSP applications to 
form the analytic signal a(n) as shown in (10), a signal with spectra 
limited to the positive (or negative) frequency band. 

• , * (10) 

a(n)= x{n)+/xln) 

We can cast the 2-path filter as a HT in a number of ways but 
a simple method, akin to the transformation performed on a ha.f band 
FIR filter (4), is a heterodyne of the ha.f band filter to the quarter 
sample frequency, .f the original fi.ter response and transformer are 
denoted by h(n) and H(Z) respectively, then the heterodyned 
expressions can be easily seen to be those in (1 1). 

hin) = >h{n)e^/2 = hWW (11a) 



20 
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H[Z) = > H[Ze*/2) = H<JZ) ( 1 1 b) 

The transformation of a half band filter to an HT filter is 
accomplished by replacing each Z 1 of the transfer function with -jZ\ 
Since the polynomials of the all-pass subfilters are second order, this 
transformation can be accomplished by changing the sign of the 
5 coefficient used in each all-pass sub filter and by associating a -j with 
the lower path delay. These operations are indicated in Fig. 36. 
2B. INTERPOLATION AND DECIMATION FILTERS 
Any half bandwidth filter can be used to perform 2-to-1 
resampling (up or down, popularly identified as interpolation and 

1 0 decimation). Finite impulse response (FIR) filters can be operated with 
a polyphase partition which does not process inserted zero input 
points (up-sampling) or compute the discarded output points (down 
sampling). While this option is not available for the general recursive 
filter, it is an option for the recursive all-pass M-path filter. The M- 

1 5 path filter can be partitioned into polyphase segments because of the 
interaction of the internal delays of the M-th order subfilters and the 
delay line of each path. This relationship is particularly easy to see for 
the 2-path filter. An impulse applied at index n Q makes a contribution 
to the output at the same time via the upper path while an output is 

20 not available from the lower path till the next index due to its extra 
delay Z\ During that same next index, the upper path offers no 
contribution to the output because there is no Z 1 path associated with 
its Z 2 polynomial. Thus the filter delivers successive samples of its 
impulse response from alternate paths of the filter. 
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Down sampling is accomplished by delivering alternate input 
samples to each path of the filter operating at the reduced input rate. 
Since the delay of the second path is accomplished by the .nput 
commutation the physical delay element in that path is removed. The 

filtered and down sampled out-put is generated by the sum or 
difference of the two path responses. Note that the computation rate 
per input sample point is half that of the non resampled filter, thus ,f 
the 2-path filter of Fig. 44a is used for down sampling we would 
perform 3 multiplications per input point. 

Up sampling is accomplished by reversing the down samphng 
process. This entails delivering the same input to the two paths and 
commutating between their outputs. Both forms of resampHng are 
shown in Fig. 37. As in the previous example, when the 2-pa,h f.lter 
of F i g . 44a is used for up sampling we would perform 3 

1 5 multiplications per output point. 

The cascade operation of down sampling and up sampling of the 
complementary outputs performs two band maximally decimated 
quadrature mirror filtering and reconstruction (5). 

Cascading multiple stages of up sampling or down sampl.ng 
fitters can achieve high order sample rate change with remarkably 
small workload per data point. For instance, a 1-to-16 up sampler 
with 96 dB dynamic range for a signal with cutoff frequency of 0.4 
f requires 2-path filters with successively smaller number of all-pass 
stages of orders (3.3). .2.2). (2.1). and (1.1). This requires a total o 
42 multiplications for 16 outputs, for an average workload of 
approximately 2.7 operations per output point. 



20 



25 
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2C. ITERATED POLYPHASE ALL-PASS FILTERS 

The all-pass networks described earlier are formed by 
polynomials in Z n and specifically for the 2-path filters Z 2 . A useful 
transformation is to replace each Z 2 with Z 2k , to obtain higher order 
5 filters exhibiting K-fold spectral replication of the original filter. For K- 
= 2, the half band filter centered at DC becomes a pair of quarter 
band filters centered at DC and f,/2. An example of spectral 
replication for the iterated filter is shown in Fig. 38. 

Note that this transformation scales the prototype filter's 

1 0 spectrum by 1 1K, reducing both passband width and transition width. 
Thus, very steep transition bandwidths and narrow filters can be 
realized by the use of low order polynomials of degree 2K obtained via 
delay elements. Energy in the replicated spectral regions can be 
removed by a sequence of filters incorporating various degrees of 

1 5 resampling and iterated filters of reduced degree. 
3. FFFIHIENT Al \ -PASS F1LTFR BANKS 

Efficient spectral partitioning can be obtained by cascading and 
resampling the outputs of iterated lowpass and HT filters. For 
instance, the spectra of a four-channel filter bank centered on the four 

20 quadrants can be formed with a complementary half-band filter 
followed by a pair of resampled HT filters, as shown on the left side 
of Fig. 39. The workload for this 70 dB attenuation filter set is three 
multiplications per output channel. 

A four-channel filter bank straddling the previous set (i.e., 

25 centered on the four cardinal directions) can be formed with an 
iterated by two complementary half-band filter followed by both a 
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complementary half band resampled filter and a resampled HT f.lter. 
This spectra of this filter set is shown on the right side of Fig. 39. 
The workload for this 70 dB attenuation filter set is two 
multiplications per output channel. 
1 p-r^E R ^M -PASSTP AMCcnfiMATiQNS 

The prototype complementary all-pass filter set can also be 
transformed to filters of arbitrary bandwidth and of arbitrary center 
frequency by standard all-pass transformations [6]. The lowpass 
transformation for the half -band filter is shown in ( 1 2) . 

6+z i _ 1 -tan(0 o /2) (12 ) 

" 1+tan(0 o /2) 

This transformation warps the frequency variable of the filter so that 
the -3 dB point resides at frequency * 0 . When 9 Q is different from 77/2, 
the second-order polynomials of the all-pass subfilters acquire 
coefficients for the ^ term. The resulting filter requires two 
multiplications per pole zero pair, which is still half the workload of 
the traditional canonic forms. In addition, the pole at the origin (due 
to the delay of the second path) transforms to a first-order all-pass 
filter thus converting an inactive pole to an active one requ.r.ng an 
additional first-order stage. Fig. 40 presents the frequency-warped 
version of the filter presented in Fig. 44a. Both this filter and ,ts 
complement are available from the warped structure. 
The bandpass transformation is presented in (13). 
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z -i => - z -y C = COSi*,) < 13) 

1 -c Z' y 

This transformation warps the frequency variable of the filter so that 
the center frequency resides at frequency 9,. If 6, is equal to 0.0, th.s 
transformation is equivaient to the iterated transformation descnbed 
in 2 C For any other frequency, this transformation converts the 
second-order all-pass filters to fourth-order all-pass structures. This 
requires four multiplications to form four pole-zero pairs of the 
resultant filter, which is still half the worktoad of the traditional 
canonic forms. In addition, the transformation converts the real pole 
(which is an image of the pole originally at the origin! to a pa,r of 
complex poles requiring a second-order all-pass stage. Rg. 45 
presents the frequency warped version of the filter presented .n F,g. 
45. As earlier, both this filter and its complement are available from 

the warped structure. 

Other transformations and geometries 17] can be used with the 
digital all-pass structure, and the reader is directed to the bibliography 
of the listed papers for a profusion of options. 



-M 



HJZ-") 



-Ml = 



(3) 
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M-i (4) 

n«0 



The roots of (3). the M-th roots of -a„ and its reciprocal, are 
ecua.lv spaced about the origin as shown in Fig. 42 for the indicated 
deg ees M. Note that we reaiize M po.es and M zeros per 
multiplication in the aH-pass stage. This means, for instance, that a 
rl-path fiiter offers two po.es and two zeros per mu,tip.icat,on as 
opposed to one pole (or zero, per multiplication for the standard 
(factored or unfactored) canonic forms. _ 

As we go around the unit circle, we visualize that the all-pass 
subf liters exhibit rapid change in phase in the neighborhood of each 
pole-zero pair. By proper choice of the pole positions, the pha* sh,ft 
of each .eg of the M path filter can be made to match or to d,ffer by 
m JU.es of 2*/M over selected spectra, interval This is suggested 
T r«A-43D for M - 2 and M = 4. The coefficients of the a..-pass 
u Lrs can be determined by standard a.gorithms 

. , -a hw Harris and d'Oreye and reported in a paper 
algorithm developed by Harris ana a u»«y 

recently submitted for publication (Signal Processing). 

Wel^evtewed the form and utility of the aU-pass po.yphase 
filter structure. We then presented a number of a..-pass 
Eolations -ich cou.d be easi, a PP .ied to «£J 
sets . The emphasis has been on two-pat, n^ut .mater, 
is easily extended to arbitrary M. We aliuoea 
aVorithms for designing M-path filters, which wi.. be reported ,n an 
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upcoming Signal Processing paper. We have demonstrated 
capabilities of this algorithm by using it to generate the examples 

presented in this paper. 

The primary advantage of these structures is very low workload 
for a given filtering task. Other papers 18,9,10] have discussed the 
low sensitivity to finite arithmetic exhibited by these filters. The 
primary impediment to the use of these filters is their newness and 
the lack of availability design methods for computing filter weights. 
This paper (along with excellent surveys listen in the bibliography) 
address the first problem, and a forthcoming paper addresses the 
second. 
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APPENDIX B 
LEAST SQUARE LATTICE NOISE CANCELING 

Example for ratiometric approach to noise canceling 



/define LAMBDA 0.95 
void OxiLSL_NC( 



int reset, 

int passes, 

int *signal_l, 

int *signal_2, 

int *signal_3, 

int *target_l, 

int *target_2) 



{ 



static 
static 

if ( reset 



int 
int 
float 



i, ii, k, m, n f contraction; 
*s_a, *s_b, *s_c, *out_a, *out_c; 
Delta_sqr, scale, noise_ref; 



TRUE) { 

= signal_l; 
s ~~b = signal_2; 

s c = signal 3 ; 

out_a » target_l; 

out_c = target_2; 

factor * 1.5; 
scale * 1.0 / 4160.0; 

* noise canceller initialization at time t=0 */ 



nc tO] .berr 
nc [0] .Gamma 



« 0.0; 
« 1.0; 



for (m=0; m<NC_ CELLS ; m++) { 



nc [m] 
nc [m] 
nc [ra] 
nc [m] 
nc [m] 
nc [m] 
nc [m] 



. err_a 
. err_b 
. Roh_a 
. Roh_ c 
.Delta 
. Fswsqr 
. Bswsqr 



0.0; 
0.0; 
0.0; 
0.0; 
0.0, 
0.00001; 
0.00001; 
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== ==== ==== = = = = = === END INITIALIZATION 

for (k«0; k<passes; k++) { 

contraction - FALSE; 

for (m=0; m< NC_CELLS; m++) { /* Update delay elements 

nc [in) .berrl = nc [m] .berr ; 

nc[m] .Bswsqrl = nc[m] .Bswsqr; 



noise_ref 

nc [0] . err_a 
nc (0) .err_b 

++s_a; 
++s_b; 
++s_ c; 

nc[0] .ferr 
nc [0] .berr 
nc [0] . Fswsqr 
nc{0] .Bswsqr 



= factor * logd.O - (*s_a) * scale) 

- logd.O - <<s_b) * scale) ; 

= logd.O - (*s_b) * scale); 

= logd.O - <*s_c) * scale); 



= noise_ref ; 

= . noise_ref ; # 

= LAMBDA * nc[0].Fswsqr + noise_ref * noise_ref; 

= nc [0] .Fswsqr; 



30 



/* Order Update */ \ \ / 

for (n=l; ( n < NC_CELLS && (contraction == FALSE) ; n++) \ 

/* Adaptive Lattice Section */ 



m = n-1; 
ii= n-1; 

35 nc[m] -Delta 
nc [m] .Delta 
Delta_sqr 

nc [n] . f ref 
40 nc[n] .bref 

nc [n] . ferr 



* = LAMBDA; 

nclm] .berrl * nc[m].ferr / nc [m] . Gamma ; 
- nc[m] .Delta * nc [ml. Delta; 

= -nclm] ".Delta / nc [m] -Bswsqrl ; 
= -nc[m] .Delta / nc [m] . Fswsqr; 

= nc [ml . ferr + nc [n] . f ref * nc [m] .berrl; 
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nc[n] .berr - nc Ins] .berrl + nc In] .bref * nc [m] .ferr; 

nc[n] .Fswsqr = nc[m] .Fswsqr - Delta_sqr / nclm] .Bswsqrl; 
nctn] .Bswsqr * nclm] -Bswsqrl - Delta_sqr / nc [m] .Fswsqr; 

if( (nctn] .Fswsqr + nc [n] .Bswsqr) > 0.00001 | | (n < 5) ) { 

nc[n] .Gamma = nclm] .Gamma - nc [m] -berrl * nclm] .berrl / 

nc [m] .Bswsqrl; 

if (nc[n] .Gamma < 0.05) nc [n] .Gamma = 0.05; 
if (nctn] .Gamma > 1.00) nctn] .Gamma = 1.00; 

/* Joint Process Estimation Section */ 



15 



20 



nc tm] . Roh_a 
nc tm] . Roh_a 
nc tm] .)c_a 
nc [n] .err_a 

nc lm] . Roh_c 
nc tm] . Roh_c 
nc [m] . k_c 
nctn] .err_b 



*= LAMBDA; 

«■« nc[m] .berr * nc [m] .err_a / nc [m] .Gamma ; 
= nctm].Roh_a / nc [m] .Bswsqr ; 

« nclm] .err_a - nc tm] .k_a * nc tm] .berr; 

*= LAMDBA; 

nctm] .berr & nc tm] .err_b / nc tm] .Gamma ; 
= nc[m].Roh_c / nc [m] .Bswsqr ; 

m nctm] .err_b - nctm] .k_c * nc [m] .berr; 



) 

25 else { 

contraction = TRUE; 
for (i=n; i<NC_CELLS; i++) { 

nc ti] .err_a 

nc ti] .Roh_a 
30 ncti] .err_b 

nc[i] .Roh_c 

nc ti] .Delta 

ncti] .Fswsqr 

ncti] .Bswsqr 
35 nc [i] .Bswsqrl 



0.0; 
» 0.0; 
= 0.0; 
- 0.0, 
= 0.0, 
0.00001; 
0.00001; 
0.00001; 



♦out a++ - (int) ( (-exp(nctiil.err_a) ♦ 1.0) / scale) 
♦outlet = (int) ( <-exp(nctii] .err_b) +1.0) / scale) 



40 



Least Square Lattic 
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APPENDIX C 



/* Normalized Adaptive Noise Canceler */ 

# i nc lude < s t dl ib . h > 
# include <math.h> 



U define TRUE 
tfdefine FALSE 

#defirie FLOAT32 
#define INT32 
#define VOID 

#define FABSF 
#define SQRTF 

tfdefine MAX(a,b) 
#define MIN(a,b) 



1 
0 

float 
long int 
void 

f abs 
sqrt 



(a) 
(a) 



(b) 
(b) 



(a) 
(a) 



(b) 
(b) 



#def ine 
#define 
#def ine 
#define 
#def ine 
#def ine 



MIN_VAL 
MAX_DEL 
MINJDEL 
MAX_RHO 
MIN_RHO 
MIN BSERR 



0.01 

0,999999 
-0 . 999999 

2.0 
-2.0 

IE-15 



typedef struct { 



} 



FLOAT32 


berr; 


FLOAT32 


berr_l; 


FLOAT32 


delta; 


FLOAT32 


err; 


FLOAT32 


ferr; 


FLOAT 3 2 


gamma ; 


FLOAT32 


gamma_l ; 


FLOAT32 


rho; 


FLOAT32 


delta_l; 


FLOAT32 


Bserr; 


FLOAT 3 2 


Bserr^l ; 


LANC_CELLS; 



typedef struct { 
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} 



INT32 cc; 
FLOAT 3 2 lambda ; 
FLOAT32 min_error ; 
LANC_CELLS * ce 1 1 s ; 

LANC Context; 



extern LANC_Context * 
LANC_Init ( 

INT32 

FLOAT 3 2 

FLOAT 3 2 



num_cells, 
lambda, 
rain error) ; 



/* number of cells 
/* put in value for lambda 
/* parameter 

/* point to array of ANC_ CELLS 



/* number of cells 
/* lambda par am 

/* min error 



15 



extern VOID 
LANC_Done ( 

LANC Context 



^C; 



20 



25 



30 



35 



40 



extern VOID 
LANC_Reset ( 

LANC_Context 

extern FLOAT32 
LANC_Calc ( 

LANC_Con t ex t 

FLOAT32 

FLOAT32 



*anc) ; 



*anc, 

nps, 

noise) 



/* input, context handle 
/* input, noise plus signal 
/* input, noise reference 



/* The following macros provide efficient access to the lattice 



#def ine 


ANC_CELL_CIZE 


#define 


xBERR 


0 




#def ine 


xBERR_l 




1 


# define 


xDELTA 


2 




#define 


xDELTA_ 


1 


3 


#define 


xGAMHA 


4 




#define 


xGAMMA_ 


1 


5 


tdefine 


XBSERR 


6 




#define 


xBSERR_ 


1 


7 


#define 


xERR 


8 




#define 


xFERR 


9 




# define 


xRho 


10 





11 



SUBSTITUTE SHEET (RULE 26) 



WO 96/12435 



PCT/US95/13469 



-110- 



tfdefine berr 

tfdefine P_berr_l 

# define P_berr 

ttdefine berr_l 

5 

tfdefine Bserr 

#def ine Bserr_l 

#define P_Bserr_l 

10 fldefine P_delta 
#define delta 
#define delta_l 
^define P_delta_l 

1 5 #def ine err 

#define N_err 

#define Pjerr 
♦♦define ferr 

20 

♦♦define gamma 
# define P_gamma 
♦♦define N_gamma 
#define P_gamma_l 
25 #define gamma_l 

# define rho 



<Mp 

<Mp 
<Mp 

<Mp 



(*(p 



(*(p 

(*(p 

(Mp 

(*(p 



+ xBERR) ) 

(♦(p + xBERR_l - ANC_CELL_SIZE) ) 
+ xBERR - ANC_CELL_SIZE) ) 
+ xBERR_l) ) 

+ xBSERR) ) 
(* (p + xBSERR_l) ) 

<+(p + xBSERR_l - ANC_CELL_SIZE) ) 

(Mp + xDELTA - ANC_CELL_SIZE) ) 
+ xDELTA) ) 
(Mp + xDELTA_l) ) 

(Mp + xDELTA_l - ANC_CELL_SIZE) 

{* (p + xERR) ) 
+ xERR + ANC_CELIi_SIZE) ) 

+ xFERR - ANC_CELL_SIZE) ) 
+ xFERR) ) 

+ xGAMMA) ) 

(Mp + xGAMMA - ANC_CELL_SIZE) ) 
(Mp + xGAMMA + ANC_CELL_SIZE) ) 
(Mp ♦ xGAMMA_l - ANC_CELL_SIZE) ) 
(Mp + xGAMMA^l) ) 

(* (p + xRho) ) 



30 



Name: LANC_Init 

Abstract: Create an ANC context 



35 



40 



extern LANC_Context * 

LANC_Init( number of cells 

INT32 num^cells, / uu 

FLOAT32 lambda, /* lambda param 

FLOAT32 min_error) { /* min error 
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LANC^Context *anc; /* context 

anc - (LANC_Context * ) malloc (sizeof (LANC_Context) ) ; 
assert (anc != NULL) ; 
anc->cc = num_cells ; 
anc- > lambda = lambda; 
anc->min error = min^error; 

anc->cells = (LANC CELLS * ) malloc (sizeof ( LANC__CELLS ) 

+ (num_cells ♦ 2) ) ; 
assert <anc->cells != NULL) ; 
return (anc) ; 



Name : LANC_Reset 

Abstract: Reset an ANC context 



extern VOID 
LANC_Re set ( 

LANC_Context *anc) { 



FLO AT 3 2 
INT32 



m; 



p = (FLOAT 3 2 *)anc->cells; 

for (m = 0; m ,= ancocc; m+ + ) { 



rho 

err 

f err 

berr 

berr_l 

delta 

delta_l 

Bserr 

Bserr 1 



0.0 
0.0 
0.0 
0.0 
0.0 

= 0.0; 
0.0; 

= anc->min_error; 
ancO >min_ error ; 
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garnma « MIN_VAL; 

gamma_l = MIN_VAL; 

p += ANC_CELL_SIZE; 

p « (FLOAT 3 2 *) anc->cells; 
gamma =1.0; 
gamma_l =1.0;. 



/* Cell # 0 special case 



10 



Name : LANC_Done 

Abstract: Delete an ANC context 



15 



20 



extern VOID 
LANC_Done ( 

LANC Context *anc) { 



free (anc->cells) ; 
free (anc) ; 



25 



30 



Name : LANC_Calc 
Abstract : Calculate 



35 



FLOAT32 

LANC_Calc ( 

LANC_Context *anc, 
FLOAT32 nps, 
FLOAT 3 2 noise) 



/* input, context handle 
/* input, noise plus signal 
/*" input, noise reference 



40 



INT32 
FL0AT32 
FLOAT32 
FLOAT32 



m; 



B,F,B2, F2; 
qd2 , qd3 ; 
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INT32 output_cell; 

/* Update time delay elements in cell structure 

5 p = (FLOAT 3 2 *)anc->cells; 

for (m = 0; m <= anc-oc; m++) { 
gamma_l « gamma; 
berr_l = berr; 
Bserr_l - Bserr; 
10 delta_l - delta; 

p + = ANC_CELL_SIZE; 



15 



/* Handle Cell # 0 

p = (FLOAT32 *)anc->cells; 

Bserr = anc->lambda * Bserr_l + noise * noise; 
MAX (Bserr , MIN_BSERR) ; 



Bserr = 



f err noise / SQRTF (Bserr) ; 
ferr = MAX(ferr, MIN_DEL) ; 
20 ferr = MINlferr, MAX_DEL) ; 

berr = ferr; 

rho = anc->lambda * SQRTF (Bserr_l / Bserr) * rho + berr * 

25 nps; 

N err - nps - rho * berr; 

^r~<- . i. /* Assume last cell for 

output_cell = ancocc - 1, ' 

30 starter 

for (m = 1; m < anc->cc; m++) { 
p += ANC_CELL_SIZE ; 

35 B SQRTF (1.0 - P_berr_l * P_berr_l) ; B2 

1.0/B; 

P = SQRTF (1.0 - P.ferr * P_ferr ); F2 - 1.0/F, 

40 P delta = P.deltal * F * B _ P berr.l * P_ferr.- 
P~delta = MAX(P_delta, MIN_DEL) ; 
P~delta = MIN(P_delta, MAX_DEL) ; 
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qd3 = 1.0 - P_delta * P_delta; 
qd2 - 1.0 / SQRTF ( qd3 ) ; 

ferr = <P_ferr -P_delta * P_berr_l) * qd2 * B2; 
5 ferr = MAX (ferr, MIN_DEL) ; 
ferr = MIN(ferr, MAX_DEL) ; 

berr = (P_berr_l -P_delta * P.ferr ) * ^2 * F2; 
berr = MAX (berr, MIN_DEL> ; 
10 berr = MIN(berr, MAX_DEL) ; 

P_gamma * (1.0 - P_berr * P_berr) ; 



15 



20 



gamma 

gamma = MAX (gamma, MINJVAL) ; 

gamma = MIN( gamma, MAXJDEL) ; 

Bserr = P_Bserr_l * qd3 ; 

Bserr - MAX (Bserr, MIN_BSERR) ; 

rho *= anolambda * SQRTF ( (Bserr_l / Bserr) * (gamma / gamma.l) 
rho += berr * err; 
rho = MAX (rho, MIN_RH0) ; 
rho = MINtrho, MAX_RHO) ; 

25 N_err . - err - rho * berr; 

, , , ^ _ _n /* *ANC CELL SIZE */] ) ; 
p = (FLOAT 3 2 *)&(anc->cells[output_cell / ANt_^ _ 

30 return (N_err) ; 
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APPENDIX D 



/* QRD 

#include <stdlib . h> 
^include <math.h> 

#define TRUE 1 
#define FALSE 0 

# define FLOAT32 float 
#define INT32 long int 
#define VOID void 



15 #define FABSF fabs 
#define SQRTF sqrt 

typedef struct { 

INT 3 2 dummy; 
20 } LQRDJPEF — CONTEXT ; 



typedef LQRD J P E F — CONTEXT Sc LQRDJPEF_Handle; 



extern LQRDJPEF_Handle 

25 LQRDJPEF_Init ( 

INT32 NumCells, 

FLOAT3 2 Lambda , 

FLOAT 3 2 SumErrlnit, 

FLOAT3 2 Gams Init , 

30 FLO AT 3 2 MinSumErr) ; 



extern VOID 
LQRDJPEF_Done ( 

LQRDJPEF_Handle hJPE) 

extern VOID 
#define Reset ( 

LQRDJPEF — Handle hJPE) 



40 extern FLOAT32 
LQRDJPEF Calc( 

LQRDJPEF — Handle hJPE, /* handle 
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10 



15 



20 



25 



30 



35 



40 



FL0AT3 2 
FLO AT 3 2 



nps 

noise) ; 



/* noise plus signal 
/* noise reference 



tfdefine MAX(a,b) (a) > (b) ? (a) 
ttdefine MIN(a,b) (a) < (b) ? (a) 

#define LQRDJPEF_CELL_SIZE 



(b) 
(b) 



27 



typedef struct { 
FLOAT32 sinf , 
FLOAT 3 2 sinb , 

FL0AT32 cosf 
FLO AT 3 2 cosb 

FLOAT 3 2 epsf ; 
FLO AT 3 2 epsb 



sinf_l 
sinb_l 

cosf_l 
cosb 1 



epsb_ 1 

pief_l 
pieb_l 



FLOAT 3 2 pief 
FLO AT 3 2 pieb 

FLOAT 3 2 Bserr_l, BserrJ, SQRTF_Bserr_x , ow^x _ _ 

FLOAT 3 2 p_lf P_ 2 ; 
FLO AT 3 2 gams_l; 
FLOAT32 epsi_l; 



LQRD JPEF_CELL ; 

struct { 
INT32 
FLOAT 3 2 
FLOAT 3 2 



FL0AT3 2 
FL0AT3 2 

Bserr 

FLO AT 3 2 
FLOAT 3 2 
SumErrlnit 



NuraCells; 
Lambda; 
SumErrlnit; 



Gams I nit ? 
MinSumErr; 



/* number of cells 

/* Lambda 

/♦ initial value for 

Fserr, 
Bserr 

/* Initial value for gams 

/* Minimum for Fserr, 



SQRTF Lambda; /* square root of Larnbda _ 
SWRTFlsumErrlnit; /* square root of 
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LQRDJPEF_CELL *cells; /* point to array of JPE_CELLS 
LQRDJPEF_Context ; 

The following macros provide efficient access to the lattice 
Define variables offsets within structure 





#def ine 


xSINF 


0 






tfdefine 


xSINF_l 




1 


10 


#def ine 


xSINB 


2 






#define 


xSINB_l 




3 




#def ine 


xCOSF 


4 






#def ine 


xCOSF_l 




5 


15 


#define 


xCOSB 


6 






#def ine 


xCOSB_l 




7 




#def ine 


xEPSF 


8 






#def ine 


xEPSB 


9 




20 


# define 


xEPSB_l 




10 




#def ine 


xPIEF 


11 






#define 


xPIEF_l 




12 




#define 


xPIEB 


13 




25 


# define 


xPIEB_l 




14 




#def ine 


xFSERR 


IS 






#def ine 


xFSERR_l 




16 




#define 


xSQRTF_FSERR 




17 


30 


#def ine 


xSQRTF^FSERR^ 


_1 


18 




# define 


xBSERR_l 




19 




# define 


xBSERR_2 




20 




#def ine 


xSQRTF_BSERR 


_1 


21 


35 


ttdef ine 


xSQRTF_BSERR_ 


~2 


22 




ttdefine 






23 




#def ine 


xp_2 




24 



40 FLOAT32 epsb , epsb_l 

FLOAT32 pief , pief_l ; 



SUBSTITUTE SHEET (RULE 26) 



WO 96/12435 



PCT/DS95/13469 



-118- 



15 



FLOAT 3 2 pieb , pieb_l 

FLOAT 3 2 Fserr, Fserr.l, SQRTF.Fserr, S f T ^" r ^; r 2 . 
FLOAT 3 2 Bserrj, Bserr„2, SQRTF_Bserr_l , SQRTF_Bserr_2 , 

FLO AT 3 2 p_l# P_2; 
FLOAT32 gams_l; 
FLOAT32 epsi_l; 



1 0 LQRD JPEF_CELL : 

typedef struct { 
INT32 
FLOAT 3 2 



FLOAT3 2 
Bserr 



NumCells; /* number of cells 

Lambda; /* Lambda 



SumErrlnit ; 



/* Initial value for Fserr, 



20 



25 



30 



35 



FLOAT 3 2 
FLOAT 3 2 
FLOAT32 
FLOAT 3 2 
LQRDJPEF CELL 



Gamslnit; /* Initial value for gams . _ 

MinSumErr; /* Minimum for Fserr, Bserr 

SQRTF Lambda; /* square root of Lambda . 

— _ „ -r • ^ /* enuare roof of SumErrlnit 

SQRTF_SumErrInit; /* square root 

♦cells; /* Point to array of JPE_CELLS 



LQRD JPEF_Con text ; 

The following macros provide efficient access to the lattice 

Define variable offsets within structure 

#define xSINF 0 
#define xSINF_l 
^define xSINB 2 
tfdefine xSINB_l 



40 



#define xCOSF 
Sdefine xC0SF_l 
#define xCOSB 
#def ine xC0SB_l 



1 
3 

5 
7 



#define xEPSF 



8 
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#def ine 


xEPSB 




9 




#def ine 


xEPSB_l 








#def ine 


xPIEF 




11 


5 


#def ine 


xPIEF_l 








# define 


xPIEB 




13 




if define 


xPIEB_l 








#def ine 


xFSERR 




15 


10 


ft vlC 1- XilC 


xFSERR 


1 






#def ine 


xSQRTF_ 


FSERR 




#def ine 


xSQRTF_ 


FSERR_1 




#def ine 


xBSERR_ 


1 




15 


#def ine 


xBSERR_ 


2 






#def ine 


xSQRTF_ 


BSERR_1 




#def ine 


xSQRTF_ 


BSERR_2 




#def ine 


xp 1 




23 


20 


# define 






24 




#def ine 


xGAMSQ_ 


1 






ft UC X. X11C 


xEPSI 1 






25 


# define 


sin 




(* 




#def ine 


sinf_l 


(* 


(ptr 




#def ine 


P_sinf 


(* 


(ptr - 




true 1. ±iic 


P sinf 


.1 


(* 


30 


#def ine 


sinb 


(* 


(ptr 




#def ine 


sinb_l 


(* 


(ptr 




#def ine 


P_sinb 


(* 


(ptr 




#def ine 


P — sinb_ 


1 


(* 


35 


# define 


cosf 


(* 


(ptr 




#def ine 


cosf_l 


(* 


(ptr 




#def ine 


P_cos f 


(* 


(ptr 




#def ine 


P_cos f _ 


1 


(* 


40 


#def ine 


cosb 


(* 


(ptr 




#def ine 


cosb_l 


{* 


(ptr 




#def ine 


P cosb 


(* 


(ptr 



10 

12 
14 



16 
17 
18 

19 
20 
21 
22 



25 
26 



(ptr + xSINF) ) 
► xSINF_l) ) 

xSINF - LQRDJPEF_CELL_SIZE) ) 
(ptr + xSINF_l - LQRD JPEF_CELL_S I ZE ) ) 



+ xSINB - LQRDJPEF_CELL_SIZE) ) 



¥ xCOSF - LQRDJPEF_CELL_SIZE) ) 

(ptr + xC0SF_l - LQRDJPEF_CELL_SIZE) ) 



xCOSB 



LQRDJPEF_CELL_SIZE) ) 
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10 



15 



20 



25 



30 



35 



ftdefirie P_cosb_l (Mptr ♦ xCOSB_l - LQRDJPEF_CELL_SIZE) ) 



fldefine epsf 
ttdefine P_epsf 



(* (ptr + xEPSF) ) 

(Mptr + xEPSF - LQRDJPEF_CELL_SIZE) ) 



ftdefine epsb (Mptr ♦ xEPSB)) 
#define epsb_l (Mptr + xEPSB_l) ) 

ttdefine P epsb (Mptr + xEPSB - LQRDJPEF_CELL_SIZE) ) 
#define fcepsb.l ( Mptr ♦ xEPSBJL - LQRDJPEF_CELL_SIZE) ) 

#define pief (Mptr + xPIEF) ) 
fldefine pief 1 (Mptr + xPIEF_l) ) 

Idlfine P Pief (Mptr + xPIEF - LQRDJPEF_CELL_SIZE) I 
^define Pji.M (Mptr + xPIEF_l - LQRDJPEF_CELL_SIZE) ) 

#define pieb (Mptr + xPIEB) ) 

#define pieb 1 (Mptr + xPIEB_l) ) % 
Idefine P Pieb (Mptr + xPIEB - LQRDJPEF_CELL_SIZE) ) 
define P^ieb.l (Mptr ♦ xPXEBl - LQRDJPEF.CEI^SIZE) > 



#def ine 
#def ine 
#def ine 
#def ine 
ttdef ine 
ttdef ine 
#define 

#def ine 



Fserr (Mptr + xFSERR)) 

Fserr 1 (*(ptr + xFSERR_l) ) 

P Fserr (Mptr + xFSERR - LQRDJPEF_CELL_SIZE> ) 
P~Fserr_l (Mptr * xFSERR_l - LQRDJPEF_CEL»I* — SIZE) ) 
SQRTF Fserr (* (ptr >. xSQRTF.FSERR) ) 
SQRTF~Fserr 1 (* <ptr + xSQRTF_FSERR 1) > 
SQRTF'P Fserr (Mptr ♦ xSQRTF FSERR 
LQRDJPEF_CELL_SIZE) ) 

_ „ _ /* + xSQRTF FSERR 1 
SQRTF P Fserr_l (Mptr + CZ~ X % ~ ~ 
LQRDJPEF_CELL_SIZE) ) 



#define 
#define 
#def ine 
#define 
#def ine 
#define 
#def ine 



Bserr_l ( 
Bserr_2 ( 
p_Bserr_l ( 
P_Bserr_2 ( 
SQRTF_B s er r_l 
SQRTF_Bserr_2 
SQRTF_ P_Bs e rr 



* (ptr +■ xBSERR_l) ) 
* (ptr + xBSERR_2)) 

Mptr + xBSERR_l LQRDJPEF_CELL_SIZE) ) 
Mptr + xBSERR_2 - LQRDJPEF_CEI,L_SIZE) ) 
(Mptr + xSQRTF_BSERR_D ) 
(* (ptr + xSQRTF_BSERR_2) ) 
1 (*(ptr + xSQRTF.BSERR.l 
LQRDJPEF_CELL_SIZE) ) 

xSQRTF_BSERR_2 



40 # de f i ne SQRTF_P_B s e rr_2 (Mptr * 

LQRD JPEF_CELL_SIZE) ) 
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#define p_l (*<ptr + xp_l) ) 

ftdefine p 2 (Mptr + xp_2)) 

#define P~p 1 (Mptr + xp 1 - LQRDJPEF_CELL_SIZE) ) 

Jdefine P±2 (Mptr + xp^ - LQRDJPEF_CELL_SIZE) ) 

5 

#define gams 1 (* (ptr + xGAMSQ_l) ) 

#define P_gams_l (Mptr ♦ xGAMSQ.l - LQRDJPEF__CELL_SIZE) ) 

ttdefine epsi 1 (* (ptr + xEPSI_D) 
10 Udefine p!_epsi_l (Mptr ♦ xEPSI_l - LQRDJPEF_CELL_SIZE) ) 

static FLOAT 3 2 
RSQRTF ( 

FLOAT32 x) { 



15 



return l.OF / SQRTF(x) ; 



20 

Name : LQRDJPEFJEnit 

Abstract : Create a JPE context 



25 extern LQRDJPEF_Handle 
LQRDJPEF_Init ( 

INT32 NumCells, 

FLOAT 3 2 Lambda , 

FLOAT3 2 SumEr rlni t , 

30 FLOAT32 Gamslnit, 

FLOAT32 MinSumErr) { 

LQRDJPEF — Context * j pe ; 

35 j pe = malloc (sizeof (LQRDJPEF__ Context) ) ; 

assert (jpe != NULL) ; 

jpe->NumCells = NumCells; 
jpe- > Lambda = Lambda ; 
40 jpe->SumErrInit = SumErrlnit; 

jpe->GamsInit = Gamslnit; 
jpe-> MinSumErr = MinSumErr; 
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-jpe->SQRTF Lambda = SQRTF ( jpe->Lambda) r 
5pe->SQRTF:sumErrInit - SQRTF ( jpe->SumErrInxt) . 

jpe ->cells = malice (sizeof (LQRDJPEF.CELL) * (NumCells 

assert ( jpe- >cells !=NULL); 

LQRDJPEF_Reset ( (LQRDJPEF_Handle) jpe) ; 
return ( (LQRDJPEF_handle) jpe) ; 



Name : LQRDJPEF_Init 

Abstract : Create a JPE context 



extern VOID 
LQRD JPEF_Done ( 



20 LQRDJPEF_Handle hJPE) { 

LQRDJPEF.Context *jpe = (LQRDJPEF.Context *>hJPE; 



free(jpe->cells) 
25 free (jpe) ; 



Name ' LQRDJPEF_Reset 

Abstract : Reset a JPE context 



extern VOID 



LQRD JPEF_Re set ( . 
35 LQRDJPEF_Handle hJPE» \ 

LQRDJPEF.Context * jpe ( LQRD JPEF_Con text * ) hJPE ; 

FLOAT 3 2 *ptr; 
40 INT32 • m? 

per = (FLOAT32 *) jpe->cells; 
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for (m * 0; m <= jpe- >NumCells ; m++) { 

sinf = O.OF 
sinb = O.OF 
cosf * O.OF 
cosb = O.OF 
epsf = O.OF 
epsb = O.OF 
pief = O.OF 
pieb « O.OF 
p_l = O.OF. 

Fserr = jpe->SuraErrInit ; 

Bserr_l = jpe- >SumErrInit ; 
gams_l = jpe->GamsInit ; 

SQRTF_Fserr • jpe->SQRTF_SumErrInit ? 
SQRTF_Bfserr_l = jpe- >SQRTF_SumErrInit ; 

+= LQRDJPEF_CELL_SIZE; 



ptr 

ptr = (FLOAT32 * ) jpe- >cells ; /.# Cell #0 special case 
gams_l = 1 - OF; 



25 Name : LQRDJPEF_Calc 

Abstract : 



extern FLOAT32 
30 LQRDJPEF_Calc ( 

LQRDJPEF_Handle 
FLOAT32 nps, 
FLOAT 3 2 noise) 

35 LQRDJPEF — Context 
INT32 « 

FLOAT32 *ptr; 
40 FLOAT32 trap; 

Time update section 



hJPE, 

/* noise plus signal i 
/* noise reference */ 

*jpe = (LQRDJPEF — Context *)h<JPE; 
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ptr - (FLOAT 3 2 *) jpe->cells; 

for (m - 0; m <= jpe- >NumCells ; m++) { 

J* of following delay elements are not needed V 



sinf_l - sinf; 

sinb_l = sinb; 

cosf_l - cosf; 

cosb_l = cosb; 

10 epsb_l =epsb; 

pief JL - P ief ' - 

pieb_l - pieb; 

Fserr_l - Fserr; 

15 Bserr_2 = Bserr_l; 

p_2 * P— 1 ; 



SQRTF_Bserr_2 = SQRTF_Bserr_l ; 
SQRTF_Fserr_l - SQRTF_Fserr ; 

ptr += LQRDJPEF_CELL_SIZE; 
} 

/* Order update section 



* 



/ 



/* Handle Cell #0 ' 

ptr = (FLOAT32 *) (jpe->cells) ; /* point to cell # 0 •/. 
epsf = noise; 

30 epsb » noise; 

epsi_l = ^P s ? 



/* rest of cells */ 

for (m - 1; » < jpe->NumCells; nw+> { . 

/+ access next cell */ 
ptr +« LQRDJPEF_CELL_SIZE; 7 

/* Prediction section */ 

p'fiserrll = MAX ( P_Bserr_l . ]P e->MinSumErr) , 
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SQRTF P_Bserr_l = SQRTF (P_Bserr_l) ; 

tmp = RSQRTF ( P_Bserr_l) ; /* this comes free on DSP */ 
P_cosb_l = jpe->SQRTF_Lambda * SQRTF_P_Bserr_2 * tmp; 
P_sinb_l = P_epsb_l * tmp; 

tmp = jpe->SQRTF_Lambda * P_pief_l; 

epsf = P_cosb_l * P_epsf - P_sinb_l * tmp; 

P_j>ief = P_cosb_l * tmp + P_sinb_l * P_epsf; 

gams_l = P_cosb_l * P_gams_l; 

P_Fserr = jpe->Lambda * P_Fserr_l + P_epsf * P_epsf; 

p Fserr = MAX(P_Fserr, jpe- >MinSumErr) ; 

SQRTF_P+Fserr = SQRTF (P_Fserr) ; 

tmp - RSQRTF (P_Fserr) ; /* this comes free on DSP */ 
P_cosf = jpe->SQRTF_Lambda * SQRTF_P_Fserr_l * tmp; 

P sinf = P_epsf * tmp; 



25 tmp = jpe->SQRTF_Lambda * P__pieb_l? 

epsb = P_cosf * P_epsb_l - P_sinf * tmp; 

p_pieb = P_cosf • tmp + P_sinf * P_epsb_l; 

/* Joint Process Estimation section */ 



tmp « jpe- >SQRTF Lambda * P_p_2; 

trap J * ~- - - - • - p sinb 1 * tmp; 



epsi_l = P_cosb_l * P_epsi_l 
35 Pjp_l = P_cosb_l * tmp + P_sinb_l * P_epsi_l; 

} 

ptr + = L QRDJPEF_CELL_SIZE; /* access next cell */ 

40 , n 

/♦ Do minimum work for JPE of very last cell 

only four equations are required for prediction section / 
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P Bserr.l . jP e->Lambda * P_Bserr_2 ♦ P_epsb_l * P.epsb.l , 
P~Bserr_l = MAX (P_Bserr_l , jpe->MinSumErr) ; 

SORTF P Bserr 1 = SQRTF ( P Bserr_l ) ; 

tmp =~ RSQRTF(P_Bserr_l) ; /* this comes free on DSP V 
P_cosb_l jpe->SQRTF_Lambda * SQRTF_P_Bserr_2 tmp; 
P_sinb_l = P_epsb_l * tm P ; 

/* Joint Process Estimation for last cell */ 
tmp = jpe->SQRTF_Lambda * P_P_2; 

ep5i_l = P_cosb_l * P_epsi_l - P_sinb_l * tmp; 
p _ p _ 1 = P _cosb_l * tmp ♦ P_sinb_l * P_epsi_l; 

gams_l -. P_cosb_l * P_gams_l; 
return (gams_l * epsi_l) ; 
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is n aimed IS: 

1. In a signal processor for processing at least two measured signals S, and S 2 each 
containing a primary signal portion s and a secondary signal portion n, said signals S, and S 2 being in 
accordance with the following relationship: 

S, = s, + n A 
S 2 = % + n 2 

where s, and s 2 , and n 1 and n 2 are related by: 

s 1 = r a S 2 and n A =r v n 2 

and where r 8 and r ¥ are coefficients, 

a method comprising the steps of: 

determining a values for the coefficient r, which minimizes correlation between s, and 

n,; 

calculating the blood oxygen saturation from said value of r K and 

displaying the blood oxygen saturation on a display. 
2. In a signal processor for processing at least two measured signals S t and S 2 each containing a 
primary signal portion s and a secondary signal portion n. said signals S, and S 2 being in accordance with 
the following relationship: 

S 2 = % + n 2 

where s 1 and s?, and n 1 and n 2 are related by: 

«| = r a S 2 and n \ =r ^ 

and where r a and r ¥ are coefficients. 

a method comprising the steps of: 

determining a value the coefficients r„ which minimize correlation between s, and a,; and 



10 



WO 96/12435 PCT/US95/13469 

-128- 

processing at least one of the first and sefcM agtjils using the determined value for r, 
to significantly reduce n from at least one of the first 1 op us, nd measured signal to form a clean 
signal. 

3. The method of Cleim 2. further comprising the step of displaying the resulfmg clean signal 
on a display. 

4. The method of Claim 2, wherein said first and second signals are physiological signals, 
further comprising the step of processing said clean signal to determine a physiological parameter from said 

first and second measured signals. 

5. The method of Claim 4. wherein said physiological parameter is arterial oxygen saturation. 

6. The method of Claim 4, wherein said physiological parameter is an ECG signal. 

7. The method of Claim 2, wherein the first portion of said measured signals is indicative 
of a heart plethysmograph, further comprising the step of calculating the pulse rate. 

8. A physiological monitor comprising: 

a first input configured to receive a first measured signal S t having a primary portion. s v 

15 and a secondary portion ny, 

a second input configured to received a second measured signal S 2 having a primary 
portion s 2 and a secondary portion n 2 . said first and said second measured signals Sj and S 2 
being in accordance with the following relationship: 

S, = s 1 + n % 
S 2 = s 2 + n 2 

where s t and s 2 , and nj and n 2 are related by: 

s 1 = r a s 2 and n^Vfe 

20 and where r a and r, are coefficients; 

a seen reference processor, said scan reference processor responsive to a plurality of 
possible values for r,. to multiply said second measured signal by each of said possible values for 
r a and for each of the resulting values, to subtract the resulting values from the first measured 
signal to provide a plurality of output signals: 

25 a correfc tion canceler having a first input configured to receive said first measured signal 

and having a second input configured to receive the plurality of output signals from said saturation 
scan reference processor, said correlation canceler providing a plurality of output vectors 
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corresponding to the correlation conception between the plurality of output signals and the first 
measured signal; 

an integrator having an input configured to receive the plurality of output vectors from the 
correlation canceler, the integrator responsive to the plurality of output vectors to determine a 
corresponding power for each output vectors; and 

a extremum detector coupled at its input to the output of the integrator, said extremum 
detector responsive to said corresponding power for each output vector to detect a selected 
power. 

9. The physiological monitor of Claim 8, wherein said plurality of possible values correspond 
to s plurality of possible values for a selected blood constituent 

10. The physiological monitor of Claim 9. wherein said selected blood constituent is erterial 

blood oxygen saturation. 

11. The physiological monitor of Claim 9, wherein said selected blood constituent is venous 

blood oxygen saturation. 

12. The physiological monitor of Claim 9. wherein seid selected blood constituent is CBrbon 

monoxide. 

13. The physiological monitor of Claim 8. wherein said plurality of possible values correspond 
to a physiological concentration. 

14. A physiological monitor comprising: 

a first input configured to receive a first measured signal S t having a primary portion, s v 

and a secondary portion n Y ; 

a second input configured to received a second measured signal S 2 having a primary 
portion s 2 and a secondery portion n 2 . said first end said second measured signals S t and S 2 
being in accordance with the following relationship: 

S, = + n A 
S 2 = + n 2 

where Sj and s 2 , and n, and n 2 are related by: 

«1 = r a S 2 



and n 1 =r/7 2 



and where r, and r v are coefficients; 
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8 transform module, said saturation transfonn module responsive to said first and said 
second measured signals and responsive to a plurality of possible values for r. to provide at least 

one power curve as an output; 

an extremum calculation module, said extremum calculation module responsive to said at 
least one power curve to select a value for r a which minimizes the correlation between s and n. 
and to calculate from said value for r, a corresponding saturation value as an output; and 

a display module, said display module responsive to the output of said saturation 
calculation to display said saturation value. 

15 In signal processor for processing at least two measured signals, each containing a 
primary signal portion and a secondary signal portion, said first and second signals substantially adhering 
to a predefined signal model, a method comprising the steps oh 

sampling said first and second signals over a period to obtain a first plurality of data 

points representing said first signal over said period and a second series of data points representing 

said second signal over said period; 

transforming said first series of data points into a first transformed series of points having 
at least a frequency component end a magnitude component and transforming said second series 
of data points into a second transformed series of points having at least a frequency component 

and a magnitude component; 

comparing seid first end second transformed series of points to obtain a third series of 

comparison values having a magnitude component end at least a frequency component- 
selecting at least one of said comparison values that has a magnitude within a selected 

threshold; and # 

from said selected at least one comparison value, determining a resulting value cons.stent 

with the predefined signal model. 

16 The method of Claim 15. wherein said step of comparing comprises deterrmning a rat,o 
of the first transformed series of point to sad sec^d series of poinp ftp wherein said step of selected 
at least one of said comparison values comprise. ^ step of selecting the lower of the rafos. 

,7. The method of Claim 16. where*., said step of detent a residtmg value compnses 
calculating a blood oxygen saturation from the sorted ratios. 

18. The method of Claim 15. whe^i said resulting value is blood oxygen saturatmn. 

19. The method of Claim 15. wherein said resulting value is pulse rate. 
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20. In a signal processor for processing at least two measured signals, each containing • 
primary signal portion and a secondary signal portion, said first and second signals substantially adhering 
to a signal model for blood constituent saturation, a method comprising the steps of: 

sampling said first and second signals over a period to tbtain a first plurafity of data 
points representing said first signal over said period and a second se ;es of data points representing 
said second signal over said period; 

transforming said first series of data points from time i main to frequency domain to 
obtain first transformed series of points and a second transformed series of points, said first and 
second transformed series of points having a magnitude component and at least a frequency 
component; 

determining a series of retios of magnitudes with respect to frequency of ones of said 
first transformed series of points to ones of said second transformed series of points; 

selecting at least one of the ratios from said series of ratios that has a magnitude within 
a selected threshold; and 

from said selected at least one of the ratios, determining a resulting value consistent with 

the signal modeL 

21. The method of Claim 20, wherein said ratios correspond to blood oxygen saturation, said 
step of selecting at least one of said ratios comprising selecting at least one of the ratios corresponding 
to the higher values of blood oxygen saturation. 

22. The method of Claim 21, wherein said step of determining a resulting value comprises 

calculating the blood oxygen saturation 
from the selected at least one of the ratios. 

23. The method of Claim 22, further comprising the steps of: 

combining with a window function at least one of said first transformed series of points 
or said second transformed series of ponjjt' rfith said resulting value; 

performing a spectrum analysis on the combination to obtain the pulse rate. 
25. The method of Claim 22, further comprising the steps of: 

using a window function, combining at least one of said first transformed series of points 
or said second transformed series of points with seid resulting value; 

performing en inverse window function to obtein a plethysmograph. 

24. The method of Claim 22. wherein said resulting value is 
blood oxygen saturation. 
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26. In a signal processor for processing at least two measured signals, each containing a 
primary signal portion and a secondary signal portion, said first and second signals substantially adhering 
to a signal model a method comprising the steps of: 

sampling said first and second signals over e period to obtain a First plurality of data 
points representing said first signal over said period and a second series of data points] » presenting 
said second signal over said period; 

performing a fast saturation transform with said first and second series of data points to 
obtain a series of transformed data points in said frequency domain; 

determining a selected saturation value from said plurality of saturation values. 

27. The method of Claim 26, wherein said selected seturation value is arterial blood oxygen 
saturation. 

28. The method of Claim 26, wherein ssid selected saturation value is venous blood oxygen 
saturation. 

29. The method of Claim 26, wherein said step of performing said fast saturation transform 
comprises calculating first and second plurality of intermediery transformed points from said first and 
second series of data points, said method further comprising the step of determining a pulse rate from said 
selected saturation value and from said first plurality of intermediary transformed points. 

30. A low noise emitter driver comprising: 
an emitter; 

en emitter current source having an a control input and an output, said output providing 

current drive for seid emitter; 

an emitter current source input switch coupled at an output to the control input of said 
emitter current source, sail' input switch coupled at its input to one of at least two inputs; 

an emitter current cource output switch coupled at an output to said emitter and coupled 
at its input to one of at least two inputs; and 

an emitter switch control latch, said control latch coupled to seid emitter current source 
input switch and to said emitter current source output switch, said emitter switch control latch 
controlling which of the inputs to said switches is selected. 

31. The low noise emitter driver of Cits* 30, wherein said at least two inputs of said emitter 
current source input switch are a power source fag said current source and a ground. 

32. m combination: 

a detector responsive to a first signal which travels along a first propagation path and a 
second signal which travels along a second propegation path, to provide a representation of said 
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first and said second signals on an output, a portion of said first and second propagation paths 
being located in the same propagation medium, wherein said representation of said first signal on 
said output has a primary signal portion and a secondary signal portion, said primary signal portion 
of said first signal being subject to attenuation along substantially the entire first propagation path 

5 and wherein said representation of said second signal on said output has a primary signal portion 

and a secondary signal portion, said primary signal portion of said second signal being subject to 
attenuation along substantially the entire second propagation path; and 

a first signal processor having inputs coupled to said detector, said first signal processor 
responsive to seid representetions of said representation of said first and second signals from said 

10 detector to combine said first and second signals to generate either a primary or secondary 

reference signal which is a function significantly of ehher of. respectively, said primary or said 
secondary signal portions of said first and second signals. 

33. The combinetion recited in Claim 32. further comprising a second signal processor 
responsive to the secondary reference signal and to said representation of said first signal to derive 

15 therefrom an output signal which is a function of significantly said primary signal portion of said first 
signal 

34. The combination recited in Claim 33. wherein said second signal processor comprises a 
correlation canceler. 

35. The combination recited in Claim 33, wherein said second signal processor comprises an 

20 adaptive noise canceter. 

38. The combination recited in Claim 35. wherein said adaptive noise canceler comprises a joint 

process estimator. 

37. The combination recited in Claim 36, wherein said joint process estimator comprises a 
least-squares lattice predictor and a regression filter. 
25 38. The combination recited in Claim 32, further comprising a second signal processor 

responato to said primary reference signal and to said representation of said first signal to derive therefrom 
an out| Jt signal which is a function sionificantly of said secondary signal portion of said first signal. 

39. The combination recited in Claim 32, further comprising a second signal processor 
responsive to said secondary reference signal and to said representation of said first signal to derive 

30 therefrom an output signal which is 8 function significantly of said primary signal portion of said second 
signal. 

4a The combination recited in Claim 32. further comprising a second signal processor 
responsive to said secondary reference signal and to said representation of said first signal to derive 
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therefrtwn an output signal havin B a significant component which is a function of said secondary signal 

portion of said second signal. 

41. The combination recited in Claim 32. wherein said detector is configured to detect a 
physiological function represented by said first and second signals. 

42. The combination recited in Claim 41. wherein said detector is adapted to measure a blood 

constituent. 

43. The combination recited in Claim 42. wherein the blood constituent measured by said detector 
is blood gas. . 

44. The combination recited in Claim 41. wherein said detector comprises a sensor that is 

responsive to electromagnetic energy. 

45. The combination recited in Claim 32. further comprising electromagnetic means connected to 
said detector means for measuring a plethysmogrephic waveform depending upon said first and second 
signals received by said detector means through said propagation medium, said propagation medium including 
living tissue. 

46. The combination recited in Claim 32. further comprising a pulse oximeter connected to said 
detector said pulse extorter monitoring a physiological condition depending upon said first and second 
signals received by said detector through said propagation medium, said propagation medium including Dving 
tissue. 

47. The combination recited in Claim 32, f urther comprising a blood pressure monitor connected 
to said detector end configured to derive a physiological condition depending upon said first and second 
signals received by said detector through said propagation medium, said propagation medium including living 
tissuo. 

48. The combination recited in Claim 32. further comprising an electrocardiograph connected to 
said detector said electrocardiograph adapted to determine an electrocardiogram condition depending upon 
said first and second signals received by sljld detector means through said propagation medium, said 
propagation medium including living tissue. 

49. The combination recited in Cliilm 48. wherein said electrocardiograph includes a tripolar 

electrode sensor having three concentrically tfrri}nged electrodes. 

50. Apparatus for indicating the constituency of a material, said apparatus comprising; 

first signal processor means for recent first end second inputs, the first input comprising one 
of a plurality of reference signals, each of saitf ^ality of reference signals being related to the second 

input; 
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means for applying said second input and each of said plurality of reference signals to said first 

&!„nal processor means; and 

means for detecting the output signals of said first signal processor means for each of said 
plurality of reference signals applied to said signal processor means, wherein said output signals are 
5 indicative of the constituency of said material. 

51. The epparatus recited in Claim 50, wherein said first signal processor means is a correlation 

canceler. 

52. The apparatus received in Claim 50, wherein said first signal processor means is an adaptive 
noise canceler. 

10 53. The apparatus recited in Claim 52, wherein said adaptive noise canceler comprises a joint 

process estimator. 

54. The apparatus recited in Claim 53, wherein said joint process estimator comprises a least- 
squares lattice predictor and a regression filter. 

55. The apparatus recited in Claim 50. further comprising second signal processor means receiving 
15 each of the output signals of said first signal processor means, said second signal processor means 

integrating each of said output signals for forming a cumulative output signal which is indicative of said 

integrated output signals. 

56. The apparatus recited in Claim 55. wherein said detecting means is responsive to the 
extremum of said cumulative output signal of said second signal processor means for providing an indication 

20 of the constituency of the material. 

57. The apparatus recited in Claim 55. wherein said detecting means is responsive to the 
inflections of said cumulative output signal of said second signal processor means for providing an indication 

of the constituency of the material 

58. The apparatus recited in Claim 50, wherem said second input applied to said first signal 
25 processor means is one of a first or a second signal each of said first and second signals having an arterial 

signal portion and another signal portion that is indicative, of venous blood, and further comprising second 
signal processor means for receiving said first and second signals and a plurality of signal coefficients, 
wherein the output of said second signal processor means forms said first input applied to said first signal 
processor means, the output of said second signal processor means including a first component which is 
30 related to the arterial signal portions of said first and second signals and a second component which is 
related to the other signal portions of said first and second signals. 

59. The apparatus ched in Claim 58, wherein the other signal portion of each of t ; d first and 
second signals includes an indication of human respiration. 
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60. The apparatus recited in Claim 58, wherein at least one of said plurality of signal coefficients 
received by said second signal processor means relates to the arterial signal portions of said first and 
second signals of said second input to said first signal processor. 

61. The apparatus recited in Claim 58, wherein at least one of said plurality of signal coefficients 
received by said second signal processor means relates to the other signal portions of said first and second 
signals of said second input to said first signal processor means. 

62. The apparatus recited in Claim 50, wherein said second input applied to said first signal 
processor means is one of a first or a second signal, each of said first and second signals having an arterial 
signal portion and another signal portion that is indicative of venous blood, and further comprising second 
signal processor means for receiving said first and second signals and a plurality of signal coefficients, 
wherein the output of said second signal processor means forms said first input applied to said first signal 
processor means, the output of said second signal processor means including a component which is related 
to the arterial signal portions of said first and second signals or a component which is related to the other 
signal portions of said first and second signals. 

15 63. The apparatus recited in Claim 62. wherein the other signal portion of each of said first and 

second signals includes an indication of human respiration. 

64. The apparatus recited in Claim 50. wherein said material is human tissue. 

65. The apparatus recited in Claim 50. wherein said detecting means is responsive to the power 
of the output signals of said first signal processor means for indicating the constituency of said material. 

20 66. An apparatus for computing arterial and venous signals in living tissue, said apparatus 

comprising: , 

a detector configured to receive a first signal which travels along a first propagat.cn path 

and a second signal which travels along a second propagation path, a portion of said first and 

second propagation paths being located in a propagation medium, wherein said first signal has an 

25 arterial signal portion that is indicative of arterial blood and another signal portion that is indicative 

of venous blood, and said second signal has an arterial signal portion that is indicative of arterial 
blood and another signal portion that is indicative of venous blood; and 

signal processor means having an input coupled to said detector and responsive to said 
first and second signals and to combine said first and second signals to generate a signal having 

30 a significant component which is a function of either of said arterial or said other signal portions 

of said first and second signal 

67. The apparatus recited in Claim 66, wherein the $4$** portion of each of said first and 
second signals includes an indication of human respiration. 
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68. Apparatus for computing arterial and venous blood constituent values in living tissue, said 
apparatus comprising: 

si in al processor means for receiving first and second inputs, a first of said inputs 
comprising ne of a plurality of reference signals, eacif .aid plurality of reference signals being 

related to the second output; 

means for applying said second input and each of said plurality of reference signals of said 

first input to said signal processor' means; and 

means for detecting the power of the output signals of said signal processor means for 
each of said plurality of reference signals of said first input wherein said power is indicative of 
said arterial and venous blood constituent values. 

69. The apparatus recited in Claim 68. wherein said arterial and venous blood constituent values 
are the oxygen saturation of arterial and venous blood, respectively. 

70. A method minimiiing noise in a physiological monitor, said method comprising the steps of: 

detecting a first measured signal having a portion that is indicative of a frequency of the 

hear beat of a living human; 

detecting at least one second measured signal having a portion that is at least, in part. 

indicative of the oxygen saturation of blood in a living human; 
isolating said frequency from said first measured signal; 
applying said at least one second measured signal to a tunable fiter; 
applying said frequency to tune said tunable filter; and 
filtering said at least one second measured signal with said tunable filter. 
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